


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1988-12 


Parameter plane analysis of automatic control 
systems using an IBM compatible microcomputer 


Kranz, Richard John, III 


http://ndl.handle.net/10945/23341 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 

























































































































































































































































































































































































































































































































a a ae e a. | ‘ ‘ 4, LS i oe . AB El ‘ a i UT | lle ee ee ene, ee | DEMS NS Bo Sy we Oa OFS: Cag & 
‘ emron 1 f x t % ty bh ee eee ee ae =) tre 4 446 sen rgied hip phipiten ind daetie lan ceni 
Day Rute #8 Oa a s oe 4 Cr a ree ee ei AO be Uh BAM a Slle ty CiOE Re Me & naan Rp Ao Gn: ty & 9.0) BA ID Be Oi Oumrernem One. 8 
. oe | om 6 9 40 ham A 1 Shah Ba aee® © bt afsadnctia Raat a arB® Ostet Vetere O7Aae Reareri am Gs t-4 Cee eee et ee 
ej 1 ats e % 1 ia ' a ¢ ia,4, re en wr) PiQatia In tg mays OCR. 9.8, O70 ed aed, ae Re Rete ® Po BS “WAaytutiies 
. ‘ or * Oar rt ae ma 4 e se 4 48 Lom ® 94A5B PB Goma Bhat ba Pathe SIMI AM 8B 0-8 Bed. Reed 18 mse BS. Rao 
on Sb | . ‘ ” tae O,Rcm 824, 49 # me * UN at Me Ger sme pePnag cv) op te Se BAe hy UM Af Pm Mags Cole 
i . . ‘ eae 1 See Mu. * or a nr 2 aoe 1 Be MMe Bo! Rhee DOR sae rsa & Bend yn choy Pheis A: 81a MeV Ashes eRe reekeh Pen 
; . 1 ’ A a e a a a % $ @r { at? 8 #iese? Aa raging & Be @® Botan. deem W OFqd: (mt Bites Be Ress @ OO Ue OO Dla & Wh 
1 wa a e > ae 4 ve ee er a | 0 Ata a trp uite eR — Coe Br ABh "AR Pty — BIA OrAtnn mas Lapeam ss anti ieigey 
i 1 4 4 Po No = nt so MEM? Or Rae ’ as Soma! Boh © tee giade he mgddt Gm —Cole'te® tems te Bote tomas 
. Ad e ' etn 81 QBs a a PP ee Pe er ee ee ee oo ae # 4d. 90908f8, ‘Ook ® er ee ee a! 
. a ae 1) 4 Pw a) e > © A Ret & beds © DG b me LApn aeuds Dene MeteanAc® 826, de %* 1 Be ¢ SAAe iy hg rhe Mes Ay Ny 1AMy) mote AB 
‘ o ? » rs be rs ite,h & ty Re 2&2 RM Vee d & Be to Nye a4 @ ee Bet AR, HH dem bs NMEA eee pep: se0ehe 
¢ @ t 4 Febas.. Ba: & te%imy ite dt. tt ® da ah On We or etter DS Ben MH 8. Re We Cs Be & 4D, OF +94 dg tows &: Fa a el Aomky 21s Me Behe Oy Ot ae 
. Cia ity 1 E> @a@wryar ee ee ee on ie Pe ee eo Aince 2 Oh Oh eat. Ge a & hth iw 04100, 9141 8. Ba HOG Oded. Re keiner @abe 
' ie) % 4 i i a.ts? & ee rOS 1 Oe Ses, m mee 2 O ASMA LARA MRO Me OF Oy Bem a A ae Get ae ee a 
. ‘ + 2 mht «Ot gage Veh cw NUL ® mew B82 tea tsb CO Aa Mey Slew Dh Gao obeg t Op, B TPL Oe nh ban Raed: MResewe 
' a er) a a 1A. ms. Pe ee eGR ‘s Mts at ss Med ane ‘ Bole tet MWe? Bese ' RAs Ae Re Re Bode GR ete ity, oe dr ads oo 0 My he 
. 1 @ eres ‘ * ese a at rah REA Meneg ee 4% + 6 Ae hoy ae #2ede™ @.%2at Om BORER, os € A0ps Porth Aine 
i . i ee a | ae#e ute au \Ay 4 oe "sh 279,50 age + 0 Se 0D DOR AD @ Bot: to B, Refihe tp Os to Be ae COO Ee, & Py hs Gr mepeves item 
: A . . 1 @tw Sw me ' “ ‘ it a) eh ee] ‘es was * 4 Vaedissree DA Vd gears Made WR NO Be te Rg atte e ly Camm 
i . 8 ry x 1 es @ \ Se tet eee oF 4 Mahe is a Meal AFgR ta — RAD 05) haw, Medel a tat Be % so AB A MAR MBee kaw 
. ee ee en 7 a we } we og Nee Ct rs | Les eRsg emsase @ 8 Bum 2AM ABO DR e os. 8 BO | ohare Sel geem & MRR. MIDE OOo aD Ow Reteme (EGR BH O6te 
5 8 ria a Og seta e Q*' H's as He Me Oy, Fe GLANS Dire hy GMAT Ad aA gg ha Al Mefee A cd! eee Rats % Mel 9.405 Poe Q 995 80 00 RO Rhy Meo”. IEEE RO 
py hen . D0 14 » Seow 0 ee i FP ahe yae Bethee Fea MORe Re ee aneh 8 04 Se Seaetnn ntehn me eth me. 
2 ez . LP 2 ’ aes seo oa + Bare Rime aekewi gy 84 FO Nee OTR OM Seti RE oS 4 My Ans Gy Asn Cee Mom, Oy Se aseltemer 
: eer hat es es eon ite th 4 th, @ © MONS ew, BoM G19 @ oh te OQ at eats tA. a6 Ty fee MAD OM Ne Ap ehNG My Opa ae te BHR 
A 7 ae 4ihs oGityts ,MoAq UpBate Mablt. 6.0) ArKag: & -_ & Fah 19418 Bog. cel Da gla tol Mehr Or Mp hy Manse Mens Goat e 
. «4 = a . hehe Aes ue bs ab hg Ge 4a 3 a am a sg th Wp Mowl tnt tgeg ase & tatb als op 0 Am ee ee Pe ee rw Yo 
4 5 ‘ aia s how Oo B.A Um 8 that Agd lec taSi dot la ty ® WENGE: &: Qaige ir, S04 2 Rhee Hg Potten ROMs Me My Re Ms we AAI FEMI Oymeg, 
. ae i 5 . as Hg Me Bete OB BF OR Uae OR Beets pte me UA. Be Me DG Retg At OD ee Dem meds WTR LANE Ay ee Qc 2 ROMER OEE Wh ns 
Fi eee 1 2 UE Rew OF OO Rey Fete Rg tee Hien My Gn Mr Ry O44: Ry Awd cate Bp Oe tee he an® 1 A te: & © tere O, OE Oe Oomimade™ 
A OF . wieet 7 6 Ropers et a yir ra * fe 
Pa A wRetat Bal 2 By ter Be QAP Dh, Wht s WR teeta, # 
‘ . . " . 4 mw tot,) © ee ee % 9 tat Seer egret 
7 ae ‘ ~| toga “a =e * Re yh ,@n & i 4OP ES «O'R MTR WRC AIPrtgln © MOD mw soMea “205% 
eo “ 1 en wr i Cs ee Oe a ee ee ew min § aig ms be et CoM Bom F IA01d ree Mire Get MELO TOs — ote BoA SOR -M8 Roman ten: 
- 5 a ' , Pa | yoo4 1 oo ee th trl |ay Fe Badge acts Myeet tA PR tu d Me ae Re BL Uy A Meche hh BeBe AAR 09 o, Be, LOG Rs Bam he We On vse Ge dome tsa, mS 
ae a , wan is ‘uw @ 1 Ay 6 1 ae eo? Ci Le hd She ee 1B Gee pws, | 2 © fe On AR e cone NR st Get e Rede hamsUMReemimtinaD & 
1 i bas . aa oye aii A ‘teh 8 ware. 0 Pp SEM Dotre 8, 8 pBORy! M.LoPeardos Rm Ap VePate Rj Orm MiGs ALD B Wy Seite a he soph Os Me Agee Uw, Mela ms 0s ed 
ts one af 4 faiai s a8 a @té9 que 418 ® ®e DsRa GPa bie hehe O49 aides OAS Qed ets Oe Gem! OM. nssOGh Bones Oeeaan. a a 
' . * gar in 6 ola. Xe Re deh 8 tm Py 2 tue Fa% PoHgrae MR, Red trek M low Pesotewem Aad My mgidi My Marlioshs Sygate 
: . Par . ' * r sat ry & é Ay 5 ot sem bete § ptary’? “8 eo Somes trade tem oft Ai Yoke ara. tors aga t one we Semel AAES 4 eet Dem 
aaa: a ‘ 7 \ 110d a6 a Sry « Pa be Ce ih ee Se ee ee er ty Lay oe eel 4A. Sey Wi MB Me & Sime ity ree 5A my hy wa 
t a fo ry yr? Aw Bre Rol bedt> 2h BV IGM MR 20m Diltace ys O60rO 1 8 oy LO WA kgeee emerge Fede sad, CMs tek Ades my «nag, bs ea eee 
a . . Epa sit a 8 Age 1 NS Pee OA Ue er aes 8 Rye ath Dy Pe Site Mat 1B POR Re 1 Be ems AM hy mnie 9 0 Os Det Be Wy My Mee Bey Ra Spite ten de Pos Maite Sade eather des ine 
soa . 4 "Roe @ue e a oF aw Ag * eae WeAse dots tet DIM & 4 8 hor QQ aur Age, @: Mater tas ae Wan eotedatne: eas as -« 
ae . . . oy 1 | W ‘aad WA ere Cae Key Be Md Oe FOES, BeRe M BItereRe te tes Meri apited 19s dee b4 et 2 Sg be @ WOO edn eran. 
. - .: . : wa ned Ghee és 1 <M ow? tots fase ve oA Mey 4 one Be *400R et Be eddy ds Aik Ie Ns lints ede. getty Bin @ M4 rOpem acne 7 
. a ae OP: ie ee Po Ca CM a a ' a0 es * + A NRA Bee Ue may & teem ty & OY Bt ley ae By dete Oi Oe nash Py, Bots Op He Mh IDEAS 2am 1M, Ry Patn > 
‘ é 2 t oma , 4@4 " Le] 1*4 Aaron O' tees ete Bm at Rt ie fF § Cet Feat ate @ stn Aer acest apy 6 Soe Cec Rom Ge crt ewes 
7 rf ' . Care ir eee ant oe '« Ce ek ed a4°R oe, & ee Wienghe, 5 ot 5D. hg ans Ogee 25P oh Bi Ree AS Ap Pte My eee amine Ry RAE Ry gt % 
P eee ‘ “8 7 , aut ey ° Mae Ds Sty 1 ho aave orite «te to Bay eae Binge te EL Yims Behe 84 sem We Rae ee. 
A A é ' aa "8 s eu 8 . BR se2% & te mint “as ept a9 aoe ee ee ee 2 wipers GM It we lam gee Bad. MaKe 00a ee Oy ae ate ony REM 
A . 4 ' ane ’ . % se sas ewig 1 6 Diaeaen 2 2 avast fH © SPAM TRAWAS ote Dewi ktdag! | Br a @ Utes Metis,’ alee Yaad Gor 2 “0, was RRe! 
A ' . ‘ oe seeds e » ii ae HO ame tg oF Ee me RE Bo wR Fe FD PRs Rebel PEA MENS Wyiihe 616 Oy GEFH Hag oy Leh De MeMAD tee Ree 
. 1 me . an . a a of & aan” % ae ee) ' oat 80 top ee Bo te BAe eet Ast. esse Adee gip © pb IMs #2 O40 4 TERD she Coa nirstp had DP Oe RDA td 
. : oho a* ge uo Ase 8 4 8e@eue vu ee ENS BAY Ha te Bat Bates BCH PCOS 1g Mes Sug ES OMEMG 3 1OUG Od Dhews Ay td 28 Me My hs SLES 9, Bap oveta f-beRrTMedy om ty > 
. ' 7 4! Mai! * aoe OD e ao ‘t-? Toran Se 8a shee @ fm Ae Bike 8 He AALAM Oh arene €) 410i © Hy eee ee a ee 
i ‘ a a "4 er o Sma on “ 5 Ah ee hee ee Pee Ca Syfornvegu ® Nem aE ab tly? OTRAGCE TempAL La FAteMs® 6 Bes Ae on we (MRO We 
‘ 1 . to. oe | ' 4% ts . oe «he Re he slams w pw, 8 Om abe eed © 0 © 2601p Oe or IH tebe APs: the PUM ETom CARSON se We tamed. AHO eaaen 
. on ‘ an t ' wee one | NS 1 tems, 1 & Me ” Ra ke AL Ra RE eRe he tes fate Be oe he heey At aye gns Cle O* Bom 
A = ' 1 ‘ " 5 4 t a ree | . ta & 36h a ee | SOS A a a ee a Retotis nl im We Otay Oop Peas Pye O Laem rm 0, Pte fr 0 heme ein D1 enhye mtr sap intel be 
a Mpeg b 7 Aly Sanaa Pri | nw toe 8g a eet EM a 8 Tee tg Metab t ae 96990 F COP Ree 4 mom Bais Oe ot RE op tity Meeks wo Of elegy Sa W157 Bd lipPu Re wimAr' dg toMy mE 
Ps ry 1 Cons . ,» G Per) ce | 8 Dene a 08 a ” TY were ' ' 44 o tA Rs ae BOOZ rAw I GRaeratiMe —— Fy Witt mgs ere weg nen hej 18) Or BA wee! 
e . J wor a et a7 4 ot “gy gt Mo Mpa te Fat Rede QAR © Wek t Oy OR OF Bet we fare tee Gat deme NH A hhegeeg Beat ye ene Date tay | Aptian 
: . 7 oe ones - COC ack en Tey wae Te ee . 1 a atone? : os oR Fangs 408 5e MBs Om hr ane fae 044d te BEN e @ hard dintanes B! Samm Ree fee Notes Gere 
. x ' . a . . « a a fe & em earemyphets kV “times of ee ke oe ey 
Par ; F F P ’ 2 aaa ” A, ooae. F beara ox 9 4 abamsd Meoeesi* 1 Bertiag® C8 ‘Ome FH Alay D S.A A, fh wlige eves WAGs s Mie es Byte Melee 
" oth ; le mee Meco ee a % . SA ato et net BN are As AP RAG 4 ow we fag Gai ty Mead reag tied em Pl 1 Ag Da De ATR ge ds Res ene 
< : ‘ F i i a on oe “A tt ‘ feca 8 Fa.s bee mw 8 8 tke Cee Oe ee ee 
. c Py % -8 4 . © gh tem He ee te Tg Meas Rene 8 getty BIA rave oF COR baa antag. MS Vem etm ene: tet Me Rect ae ng 
: ‘ r tuo . ‘ u von * a Vig ee LAM el Da ORO Bhan ote ~ temuted iMeea Avvo et MOD Ot Re aN ym wees Me Bite tele ee! 
; - ee oe « . e a ue ' aie Pare’ Se a ee Vy ee ee, er PY 1 emetghw On 0 Oltetag Beet Od ame 6 a amatan Ra tates We Bs tend © noene re 
. Py a rf aos ‘if is re oa 4808 Me) ee, ee oe, ed fees 8 IM Et fe eee ewe 8 aa ms cht ap bahens © er Me CR Se Seihe 8 Wee @ Ogee: ete ee 
« . a ~% e . . of . ‘ ‘ ara. aa eo. r 1m, toh ates 6 8 ee ee my 1 Dame Petee gfe Metra Ty ae KH PUD os te Mt 60 Ae AWARE he Bete RO ae ME 
. a. 4 wots aa 2 @ a Ripa + oe fF 4 va “ae Metdrenet sop mcg ee tate srteaeen OO te We On ees Re Re Oe Ne Aa nee gn eV ie int ocd 
' ot "@asn dak awe SAFael « m ¢ 54 teh fo 2, We ROW ote Ty seg Med Rae eee RGD OR ay De Saye TREES Megas any Wi ATs at 
. . 18? so "w a a . tgs ee Pet, ee ee, ee Se OP eh eee eT 
. . . . 1 ‘ae ts . as 1 6 ot of ads MS we L] Wl UK A ee 
ks . ails a F ' 8 s . ye 1 ‘ = Pre Tee ey ee Oe ee To ek a | wot tod Wave mere. sit wee vee war a 
. 1 1 * ate ce eet Se Scr a Yuen * mite et © be WA te mam Sgt g ope np legit & eS My Moke FR A ohh tee See Dee ah 
E ry i is F; 7 5 1 “4 ay . a . tae orm « ~ a ee 4a > Wein Ge Anean har tant 4 Saga" freed sma aleemy, gree os 
(i A a . . ’ 2 i 1 2 woes ‘ 4p 2 RPM A ERA oF Agd ROL Reety ee ne Ane ty Dee s we Rete ets me ee 
. . Reriet ine es " ' at, . 2: 5 0) 20 he tema er ae. Rvetes, Sfuvueaes 14 wit . tS SRE ee cape 
. ae aa ‘ roms " os © ire ~ wv he Ane ARst SS ade O88, O50%, Me we 8 store ON ed 
re a ry 8 at | ae Sd % a ag? we mest WW ket 6 gate adam tonne mi boos redeemed eRe why tote ted = 2ay? 
a ‘ . ’ oe 2 8 Chee fous i! . 4 a @ "ty » 4 Pe | mee Ne Og 8a ko te eres te he ene ua eu ey tee Se By EE Oe Ree ete Bate ae 
. ae 4 rn) ames, Sed the  “« i ie Oe “fe <M Rete Be temo 
- Ae * = is hey rd Aye 4 . ' wef g there wee fgte een “'atnhots pron wae w fee Se ee Oe Re oy 
, . sh he 8 as Pr oy "*e * ' atin 1670S oh ‘* « Weidens Fekel ow game thd spss hden %-1- 
. Par 2 a po. 1 ’ Pier) a . = en ye Reds mf . tofiteaa nt Steet te teh ep! 2 Tm le oe 
e « . ‘ 4 ' so 48 . * 6 vet the om Te ee i eT 1 wef Rte sene WE oli 
. ‘ = A enya aa wacdatg Of els + gti pane Tete SoM ls Set of mip ee @ we ty Ne tei 
é i ‘ ° Peer ' i i} a oN at been cae 2a ts 
p 7 = pose a ' wg Xe segs Whe pies be wee i nets See te Toe ety OESe 
Fe s as ' s F A oat mae e« ‘ ow dea ae Cn wy Comin geor ohn s WARS | 
' ne ie fig] « i . e Bod eaesast Op ‘ SR ee Ate e tale te aed 
oN i , ' r) ' n = id Perth as Oi tiwicehs” qttetide a &." BT amt SE tt 
: . es ’ Te ee ore ier! wee Pern a rk i eo te baa 08 Ke ee 
"—,° . Oe ie r *,¢0' <, a Adne stn dededi ened . fos 16 ft wint AI ate ogo, ihe! 
% ' » Guar ‘ rr] en . wat hs bahia | t= abt ate %s os 
1 . 1 4 . 1 ay . . 4 ae e918. resem, oO! FZ patti ane nat white oaduse 2 tm nth S" ASD 6” say. 
. a 8 ' a Cer teat 1 “4% ery Td ee oe ee ee ee | cry oP ee ty ee “hy sites whe Sane 
Pe a ' . aos ar, . - . whe « 2 t adaergtcasBn mm tyne 2 240 ta cneg erg hath Meme OF UA Bae $e Pi’. who le 
ne 5 a ‘ a & Pete - hae 2k eee, eh ake b ge gas digaan ns eet mw haett We le ong wee SS Siem Bei 
ave : t 8 oe Bots? tity etl tae ° dearare bbhee ber Boney 
' “4 Sy rie . . Pf Ce 6s if te Ae dee = Wee kee 
8 : ‘ 
. ae U) ‘ ’ ' s 
: oa ’ cts A 2 ‘ 1 , 
. Stns =f Acs are aha = ‘8 
SY gee 8 ‘ ae eyes . . 
. fa ‘ i . ‘ aa 
8 . ' 7 1 a a fae oe 
soe ae ' : . kel ts 
7 i La | - 7 ad Ay ieee: os 
. . 1 of « 4 % : feo: 
. ' ae .< Palos 1 
at ° eo. te base Cha Pat) 
a o@e “ae . 
. 4 " ‘ 1 on of , 7 
- 8 1 te c*% * r 
' ' ’ #8 a4 r A a bane he wo mun vied ei 
. . . fe tame 1 - . a ' bv ney ; “tee we é ofa" Aw 
. « a ' Car ve eM ee a ° Py ne we gene bee ww ae 
‘1 ‘ . + . ,* ' sear ‘ S Set ere Fas” af: um «oF 
. . . a 6. haa, . e4 ya 8 an Shem twet Shaw’ 's 
1 Pa « . eunete ie eee re Ie — fia oy, 5 ff 
. ia - a 1 ' . wk Fe wide 024% WS 
. oe « I tee t : ae ee 
. . i < Cera te ORC ae a ar i ee 
ae . . iar Care oa 
: . A ; mas . 8 . a 
. . . ow . ae 
« oe ' : . 4 He i Soe eee 
a on i , a et wm Wah. ct 
. ' . 2 « ’ ‘cs fat mee oe at my LAR the ‘a 
. ira. +6 . - , Pee ate al PRS 
1 | » ee é i] 1 CF 
° ‘ io0 . 1 aoe oon 
"I . . ef ra 1 ages % 
i] om . ' . 1 oa OR . . . 
. mrt ' eur a x. ‘ . A” 
' . . . ' . . : . o «4 1 = te 
. . a ny ' te ts ut ares ‘ rit aae 
ae an a: ra FS a 5, ' ary 
' ae eee Pa en 1 A 
. eu ' oo: slbge, te A ome mea? a 
. ‘a ® é : fac . 
‘ ‘ 
« ‘ 
oe. 8 ss 
: : " fg tet" 
' Ya (27 
z San © “ety 
. tes “2s 
. Cena 
‘ ry ' . of #8 
. wae ate 
eas 
. 
id 1 
; . 
ne. s 
1 . 4 
. tua 
me a a 
. ten 
1 . 
on . 
8 . 
. ‘ . ‘ " Be eae PAs 
1 ‘a ' ‘“e ge etytecget we 8 . weet nb eer sae eW ap ge se ny rere 2s. Pee cd 
‘a Cy Coy ee | of 1 ‘eo "sy ' ane re : vue besedéat? ase pte ase gt ae@4 2004,35* . eo cee 78 ee 
ee? ? F = wast ng 4 prec te Fy tt ates =e es Fou rer Fg IR aqotre ge, Sat gt 7 ot ote lady 
Pa . ace a Pie 1g mo 8 tage gsy 27eTe hg ace fe 0, A At Boy OH turer 1 
‘ a aga Of pee r bey aries fete Boogie wate wae Fe ee el 73 “-* 
. eae etowe Ra R EC Oe CLC oie! 2 5 lie gaG ave oO ets? ws LS alee 
Pon) ope tr etfe Ly ee « ” a gt geen, acm te 
' * ia toa . at 4 ote cer 6 ad Fg 8 gost © porter 
1 on un ater oul Mapp oa este we ae 
‘ ‘ 1 sare Ot Speer eater * Aare gee 
1 . an at ate . eer A «fs te 
on oer ee ft gta wt eh aemrpt sae gt 
4 . tae . . to oo try ot 2 ea ne eee are Or 
ine i a a) ; ‘ ' Ow ae te Dg aie airs wee ta 8 rere g LA oe 
A . a ee et ee 1 ta rier | toa . CU? Wb pms Her Gre CM gr | ee ee | arpen 4 ted MARISOL 
eres ‘ 1 Pe ee te 0 a = 4 eae ar Cary bo opatary Ci od ee ee Getty aeea Pee F bt 280 P7SLPI Ne Rie Oe, Br gt eaet etna. mw mye ee 
bal ‘ r 1 . eo. ot . a eee ol nt ' 1 on ot 86 ate st goatee at ~wPeten aoe ' he #8 6 fa fF, eet ree Rabe Bcataene otal af re ee 
1 . en oe . Ped eC "papa 8a Pets p*qge gae oo MO oe ery CR ee ed Ce PRTC Me gh Pele gra ah aol Fao g ee eEE say 
ee aapeeun ae. Sane thee dee) wal ene hs ap staan ' , I rch lhe Rpieeother ne fem evelieet ams ce mer a 5s 
. oo . 1 Pe ee | ieee “p48 ooeTae ° 03a et * 4 siceer as z 
1 Boe . ore . ‘aiiwe err dea “ 1 . ‘ = nts ¥ or QR OT Mey TY 
' a 1 re ar ‘cc oeaiiece ab apie rae a Pa ata ‘ ee ee ee eu Boge r unevate a oatee ae 6 eet e a rat sg eavat ede rhOte 
‘ ‘ ra. as Pi 1 ‘ ie Le . Ar) ' oie ea 8 4d yh ee 44 sfFfs ¥ © pasties © peer Webi Ds o a ae Ate aed 
F; ~ A en 7 cw a eles ered] Armen aa "4 e a eypteon ns te es ono othe o abe p 2 tn ee sdveve wld al et ew # Wide pe eat peat 
1 tne te 7 . he ate la pay fae | Paton “gene sae 4 aye Pte aa ANSE ow ere ot ge wets get See ee eae 
7 - aa tan i St- eaea a ‘ ‘ Ceci tan rnge at anata Mme Spree se arntdeae $ 4h Rerelpr esewes peagee Hw ae aR! oe an ye te 
on oc Bie eer) aoe s e@-8 ae . " . e = ¢ re a | eases al wns fe see, "gt fon 2229 8? wage 0.08 gvdet 
. rae 1 ens . aC a ry rear | 1 Fareed a Je eege sy pee ere tee e ry Bev ahiam So- twene® Memon ge pave emar 
. bee i : a . » # ' : tee oa ee Pe tae ? a 
7 . 1 a 4 . . A 68s s ~~ er 8 te Ce Le oe ie i] othe Soph Of 4° Fee atk we ee aed y ve <9 of phy! OF Gr 88 Og oer aden oe, 
1. a a 1 ‘ ‘ Py é Kieagia ieaees oe Pw pe 0) ae A eine aa ay oe Fe tte 6 FO Ot wm at pe titers are Le all f Si gtevat state wm we 6 org ft He Pt 
7 ane wg as aNnece ae le ude as me A ne Wed «sagt Pe i eer ee Pe 60 2 mg inet Vipmatgtes ae 6am FF oa* Menara Peer re aterae 
aa . Paar tr ee Pere 1 fe ‘ aoe ane oa it ee eat eteta ea we ty "eon a prate of ay “p.2 ete © tated a omy eters toes Woe Gere Obs adr peat Wane 8 ee te Chel yt ane 1 Many 
. 48 ta ‘1 tera oe . son eeett ern Oy eta portage a ete tthe ft ep ato et etlépes @ Freee we ehyen Cmesiarer ae eG pede aa gt A tan wlan}? hale eS 8%, 
an o 8 84 a er ee ean Cie ane Jer War’ Sue Teer a a Pe 4 del ea afe fate @ wt Fae? iP gm s TT a ete * DP FEST Ft reat gant 
. ae « at oe tor ar] ' stp * oe Galan Sa ee te one "4 ewe 48 a AO GNSS Rete wrhete tee pentne weeas ema fe ee ene GE ae 
fe usa es gee ee joa Be es rea et at fee Gi cies tie an Ena fe Dare eee pee ge SR Dd Tt eae 8 Fat g 48 be AON ACL sae gu sed ee Mt 
1 © 4 dion hui 2 apte © aahe a ane Oat to ate 2 On tee 4 ta af fe gaare ae cele Sie a tithetatmiat ating © sP eenwemal reget 1 ate eWnen dane me age mn 
a a = ' ae Geos se eu. Cr ee ot eon ' jf water a oe Ome be alan sen brated pa riar stase Od hk Did Pf MEI ae AERO EO Om 
' . my . aa ' wena 0 oe aede Fp ee re ee hel ee oe ee el Lee ee Oe hes) Tenae Ot gee eR ee ea! ee een © ape lew ee 18 
7 ‘aoa : ie ‘ Pea ar Oe Cae | 1 Oe de Soe Det Le Le BT Ie slay ey a ptatee. ora erate 0 PS Palen af Ce ate ed 
a 1 ee er on ems sO Biya oe e faneea Oe. Bo rays lage EF! we? Meade gs Lael iF 
i; ua co te ear fi ‘gaye 8 Oph ee ee tae os ee 4 Mase pete eretipdy be alam FATTY PET a OE 
e ‘ ae owe = \e ot dnote nd atavarre eo rye Mergent ae pagead tpt oad ae awe A fh eB 28U Moen, 
a Pier | Fi a4 . « etacue Cr, Pesca eek da at eee Sat ee oy Bec let hae aD LU Gort asptoutre ae SADE eee she Ose 8 at san Oe ge py oF 18 wate ’ 
ae ae ea ba a)'s,.0 6 ate gcho os tw be te gel wide ae etanng Peg tals da a7 Ag Me at weet eee am ake do. al ere og ae seers wa % 
aoe ' ) . a * 8 oe, 8 yaad "an COU arte ROP ere ieee ee iri te ee eee ar re 
. es ta oni s, 545 Gaee raee 4a twa coun 2 MORAN 6 eg Etna, AR te ot phate? BT Oe axe aire 
’ ' s tt peua 1 F wm mt Ui WENA ME PONE TEt pe Sw apah a cuee gang 
a a 7 Es .e aye “°F tea ace ery ee er 9 aha gee ce enti ye te OAL OTN eg ap 
0 tee aa ' tw Mie ore 18d eta ee dena arp atetgey tee Ce ea er ee org gee GPO ee Gta UFO pie ge FOR etree eharmame® a 
i oa 7 Fi F a ; > Res + 30 is peoen wa ate Bee hate sadben 4 a jae a Oats fee et eee ait Pie ed tlt ee Lee er, eee 
1 F 5 Pu Tet] 2 oa Mao ies tise, i oe 1 yr etheae ys aet@ ee abe gp ga tneps ete ta 2 eee Peete epee te UP “2 Spe B08 8) 1 ot of got etre rh pet ata mee gm hat. 4* 
a a a rs F pa Qe 8 4 aie a aes 8 ae ae the wg ee ee A a er a bl FE arate omer ata eat FW eek ge gt ATE Lite a gt arp eat af 
F ba é ‘ toa et Caer en hs eats Per ha | ry 2 . ph pt atated ar ee FOB ee Pt POON 8 1 AEE om mt ga Paty Fge 4 
ae . ‘ Ps r - a y e abt 09 wae rz thar Pyta ara Nett fF 
ea ae ‘i a i rae err » %e ted tne ae att oe on ore oan Sat baal oh A dak npr axiele ts ah cal ial 
; nen : P rergare haat Toe Qiale efattg . a Mk) Sweet toa 8 ei Oeaeniant 0 Fett ore Fa Bare a ehe wkend yt hige yeaah aed she VERE ERiWe rae Aw. Peign 
7 ay a ee oF P ‘ 12g we ute at te kg tlhy atta ate 0th wth sg. re ee a S918 eT OO gta het at etme gp ba ae O° EPMA 6r Fog cnrkey 





PEt Cat et ge EEN © eet ale! OS ye ge tears tlh agval oF prierg abt S! OTRO 1d PRO! Adel geage,w 

BOER TA We RP A ae pt HOO Vi ak sigh got bd pa 08? Wp Eat it od 

PAVE LAGU INR Oba wae BET s ee Sak pha usd ae AS FUE etd ee Vo waal Roe 
’ 






Sgt eve Pateeate ace gs Mee tee a te Beem WEG Bateman sey terete 4 
paced Tyas 8 Ot ys ate of i 08) ge ange gan Paeeets a at ee thee Phang 
I> & ase pe Pests e Pal eVaenee a” 





















ee . . 4 ° 10 “s gare i PY a ae are Lo lat al ed ok ee te 
‘ ” ae ort @ ‘tev ie Lie ay att GAT e te ghat CX aw al went p88 s OS ee Prt wat fh pare Mgt 
‘ tote Pe 8 tte sh pee ahteat) 6 MR 1 ste tg pat LEV eI EEN OOF A ON FD 





Fee ge ara agente apa nre gata 
cor Ny feats bth atten @ Log aéugtr G0 8 gh ae Aid- Ohio ah IF ath tee Vite ttle seal a eile 49 Spi, Fae? 
er ie re ir wel Petey Stree ei Lr re eet ee teil.) se ee re ee aa de 





© oad4pel 
163% 4 0 otafevisee’ 
btae Rests ote o 





1 
4ha note ow ahd eth ety satg OF 
a on at a) oho ie Ay '’ ' 





6 






































waa shy Ota 08 Mgt yg MONT OsOT HAs red SE ® Bathe taser ane Oe abl bt pe Ut ge PON Pe gee alo oO 

1 of sea rabelia 4 by eo ee ee De ee OP ed ead TROL LTT i ae ee bl ad tel RO ee 
aoa'tt ‘ be nnitia 84 oy SL eh Meh Se he eee ee OO OF stew sande At. pale hake Bats Cg ad T97 

. af ar arres Gow dy bee cr i er | Ane! BOP 0 Ih Gabel pee egilsSei® 2" ep oF stout ames 

sa ht tae s guetta: & met seit ar % ’ mee even ate ten Te ae UN re watts ae fae Pe ie a eT ool 
te seen 8 Serta cee Hote Bh ge gee a Geeks 6 te bere eg s is ee ie ed ee ee oe ee ea a Se 
af wee LG RUE E ET VR ae Ete Geet ONE oy ate? BLT ee AAR FOREN OKAMOTO” ry hae gg Ter Pet AO 1d Si Ose se 


> eb aas s afh 2 ot 1 tot.4 6 betued bere 
arse. cee DB be eet 4B 8 sew wIel for 8 g beeen By yu aed 
tara ge tiar e UQSRNT sagt Wh 





See oe see bytes 


REO bar Feld oC in gouge rhs gs wt pei pay etl ie sQt GES Cpt yevenns 
Vers tarpelyperates a 


4 
ae ae ee ee es Pe ee Ot ee a et od dal ea 









































































































































ol Fj “yy ate eth ba mee ae rahe hehe te is oe eat ee emule en wire thie Let tak te 
i . Ine . a A at emt sav aan eee sd 6 ot era "gery CARS “henry ah 8 tps wo we are Soutres sat Oy eee er oll 
gan 1o#eg tot et & ho a wo igeat ie oe ae oe 2 PT eh ee ee wi 
aly YL seta ths aye DBs S watt gaps pretal | 190520 Vente ar & PAD GPT DP ieitg ad APT tana LGW Ma ge, 
1 “~ se fet. “ae 'k al yree te neater oven otta er Poe ee ee ee itl 2s eee 
n't fri : ue og te Ok ee OTR ara ie ia heed ee ETT eee EP ee oe eee 
‘ ree Been herd ea ere ar an ee ee es Oe hd eo eT i oe eee) 
' oR) eg Aye PaO ea he Wee Te Caw OES hs pate: rT reas tte Ft oe eka yd 
1 err “4 4 oe eee Pas iimanm yea tater meen agree 4 he ee ee ee ee oe ed De a 
. owas «Weare att aie 4 a] da sete Ree ef $ Li ase Teak Ui era h 
pelan Ulett Dade Fe Ra heat F Pith ee iy Oe ikem Tere Pre w dhs e Sy Ye 
e 7 Reet wleoeriaee Ape ‘et face § teases? % Om penta or oe eyo ewe “pGad ge eS ONES 4 0s tel aig ials 
oe 4 iar ee ee tL ee er Cn oe ed oe ee | FBS Ger 18 14 NEUE ETN EH Dee “he CALS TE oe betas 4 
. f mis a = (na) mgd NCD Wy ew ete gry oe ee ee pipe ensey see Vat ore Pim trs ty Oty Sel Me eee ee tw Pee ioe ST se ee Tare he ee ha ee el 
af 0 oe 9 OO aT ae Oe BT ene ar ihe eeteteye Sek ATO ay EP ANG Gerais epey rina hae aye? OTe Paya ee sabe 'eset seven 
sa 4 be otsee ot sathats $4 Of WE aatd se aye opps te Red Fe Seana 81 OTA Le Syl g EU SN A EO 889° Fame Bhss 
ana Sa ‘ Pu Bh eahatned 2egare ay! 2a be wade Cat es PTA a & ae) Voted es ypatee Les othe 0094 a! et era 
a a aees eee Meds 84 yw Gere oe } a (eens @8ea Sem Vag weep iles ph eR Ae Meee site 1B AIMAER SS HET GOUR 
8 He 8 tee 8 on ngham yee FE te PENNE RET 8 BT BON TREE TH er pnzede a Te i0h ade POPP U TE MONY fa Meg BEU YMe Oe ow ata pees pone 08 URS re 
uo > offer a we Pye = OM Cad ote ote POE ce ere Mates CWS Oy Bem ge ental coe gigs raed 12s Oe Teeth be! pan ty Te? dye Ex eee i ee <ewaeeee 
serora “sue 0 maa) ose 44th baad Ge rr ht eck oe ee Oe et Peele Vinee eke Or tod Cee Ca ae PL ad a ll ll 
. als Fier Pe Co es Ty ee i ou Pape were pts oo te ata be es AT byte sil Lelie 1 Cee ele wig Rg le a 9 renin 
. ive ae ’ 4? oO "yee © Ohad oe at a ee A wep aM ae'p Hag CHRP Deer Ras Wl et Re Poe ee a rs 
' ' a & Eager! te a at oy yee Were oF as ge rou .£¥GE ee < te gods. Bae eden aT Lay sl scomen pra ee boc ™ 
; fee fh oe - © ieetade Pe rere Pee Py Oe ee pa 289 argh Ph yah U4 oho tit Oth ya ate ee oe setter eh ae 
t a ee Ce re Oe See atone TE Rie er Cire cs Pech a oc TU ee cea ol SUC Ee RG Lh es eee Ca Peer torr et : 
‘ ee ee kn ZANE a pee ey ba pret ahace boo everest as aren, 14 2898 Re Ra gy EBA rebefey Ce nemaain gt MoNreNe 
sia ew Pa re ee Cee UL Ce eae Me MC A es De ee ee Ue dl pik ae Santee e Ae oe 
ho wae ws ‘ Py POT nes eer I a a Te ee Pi iat aie tak EM th hee hon ed ue ae « DITA ERMA! oF rar pt gt FOVAHUGS 
, - 8 th ¢ oS ne a arate & 4 Cee ee Me LY Diy Da ae PMID ed oar rales Ge This Site LIE E aia et ht Lee a 
i < Bet Wana m4 1 fe ae¥? best a8 ” Oe late (UR ANY OWE ILEY anh iT iE teode Lees Oava qe brinaintestesiee 
. 1 ‘4 Sa ae Q4e%s 2) paths tense see "qn 3 aP oS ho oF, thet Pit n graricgletsuaa ir are 
ry or ae | Py |e « £ ees 08° Ecce we wee we ae OED LAUT LION * po aA E Lagcan ete FM cy Beaty ye Bet ULE OG) Ve egret elpE eke 




















NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





THESIS 


YAR IE 
= ? Nd Lad 


PARAMETER PLANE ANALYSIS WITH AN 
IBM COMPATIBLE MICROCOMPUTER 


by 


Richard John Kranz Lil 


December 1988 


Thesis Advisor: George J. Thaler 





Approved for public release; distribution is unlimited 


T242018 





UNCLASSIFIED 
ECURITY CLASSIFICATION OF THIS PAGE 





REPORT DOCUMENTATION PAGE 


a REPORT SECURITY CLASSIFICATION Ib RESTRICTIVE MARKINGS 
Unclassified 


a. SECURITY CLASSIFICATION AUTHORITY 3. DISTRIBUTION /AVAILABILITY OF REPORT 
Approved for public release; Distribution is 
unlimited 






b DECLASSIFICATION ‘DOWNGRADING SCHEDULE 






. PERFORMING ORGANIZATION REPORT NUMBER(S) 5 MONITORING ORGANIZATION REPORT NUMBER(S) 


a. NAME OF PERFORMING ORGANIZATION 6b OFFICE SYMBOL 


(if applicable) 
62 


7a. NAME OF MONITORING ORGANIZATION 






Vaval Postgraduate School Naval Postgraduate School 





x. ADDRESS (City, State, and ZIP Code) 
fonterey, California 93943-5000 


7b. ADDRESS (City, State, and ZIP Code) 
Monterey, California 93943-5000 






la. NAME OF FUNDING, SPONSORING So OrriGessyIViBOL 9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 
ORGANIZATION (if applicable) 








ic. ADDRESS (City, State, and ZIP Code) 10 SOURCE OF FUNDING NUMBERS 


PROGRAM PROJECT TASK WORK UNIT 
ELEMENT NO. NO NO ACCESSION NO 


1 THLE (include Security Classification) DARAMETER PLANE ANALYSIS WITH AN IBM COMPATIBLE 
4ICROCOMPUTER 





2. PERSONAL AUTHOR(S) 
See, Richard J. IIil 


3a. TYPE OF REPORT 13b TIME COVERED 14. DATE OF REPORT (Year, Month, Day) 115 PAGE COUNT 
(jeter's Thesis FROM TO 1988 December 14] 


6. SUPPLEMENTARY NOTATION The views expressed in this thesis are those of the author 


and do not reflect _ the official policy or position of the Department of 
Jefense or the U. S Government. 


7 COSATI CODES 18 SUBJECT TERMS (Continue on reverse if necessary and identify by block number) 
FIELD SUB-GROUP Parameter Plane, Automatic Feedback Control, 


Td SSCS~*S Mttrovic's Technique 
a 


‘9. ABSTRACT (Continue on reverse if necessary and identify by block number) 

A group of lesser used analog control system design techniques, generally termed pa- 
rameter plane methods, is examined through the use of an IBM compatible microcomputer 
program developed as part of this thesis. 

The coefficients of a system's characteristic polynomial are determined by the plant 
and any added compensators. As these coefficients are varied, so too are the roots of the 
characteristic equation and therefore the system response in terms of bandwidth, settling 
time, etc. 

In the parameter plane method, a designer selects two parameters of a system's 
compensator(s). The parameters commonly represent such attributes as a compensator gain, 
pole or zero but can be any linear system function. One or more system characteristics 
dictating desired system performance, suchas relative damping or undamped natural frequency, 
are computer model inputs. The associated parameter values to achieve the input 
characteristics are output in graphical and/or tabular form. 


20 DISTRIBUTION / AVAILABILITY OF ABSTRACT 
By) UNCLASSIFIED/UNLIMITED [ SAME AS RPT. {} DTIC USERS UNCLASSIFIED 
22a NAME OF RESPONSIBLE INDIVIDUAL 22b TELEPHONE (Include Area Code) | 22c OFFICE SYMBOL 







All other editions are obsolete 2 U.S. Government Printing Office. 1986—606-243 


i UNCLASSIFIED 


Approved for public release; distribution unlimited. 


Parameter Plane Analysis of Automatic Control 
Systems Using an IBM Compatible Microcomputer 


by 


Richard J. Kranz III 
Lieutenant Commander,’ United States Navy 
B.S., United States Naval Academy, 1974 
M.S., University of Southern California, 1982 


Submitted in partial fulfillment of the 
requirements for the degree of 
MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 
from the 


NAVAL POSTGRADUATE SCHOOL 
December 1988 


ABSTRACT 


A group of lesser used analog control system design 
techniques, generally termed parameter plane methods, is 
examined through the use of an IBM compatible microcomputer 
program developed as part of this thesis. 

The coefficients of a system’s characteristic 
polynomial are determined by the plant and any added 
compensators. As these coefficients are varied, so too are 
the roots of the characteristic equation and therefore the 
system response in terms of bandwidth, settling time, etc. 

In the parameter plane method, a designer selects two 
parameters of a system’s compensator(s). The parameters 
commonly represent such attributeS aS a compensator gain, 
pole, or zero but can be any linear system function. One or 
more system characteristics dictating desired system 
performance, such as relative damping or undamped natural 
frequency, are computer model inputs. The associated 
parameter values to achieve the input characteristics are 


output in graphical and/or tabular form. 
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I. INTRODUCTION 


A. BACKGROUND 

Over the past decade, a resurgence in interest in 
analog design techniques has been triggered most notably 
with the advent of the switched capacitor networks. These 
networks combine analog and digital technology on one chip. 
Proven analog design techniques which had been eclipsed by 
digital design methodology have had new life breathed into 
them. It is widely recognized that they are still valid 
tools in the design of modern control systems. 

This thesis will examine a group of lesser used analog 
techniques termed parameter plane methods, parameter space 
methods or Mitrovic's method. This method enables the user 
to observe changes in the dynamic behavior (e.g., natural 
frequency, damping ratio, settling time) of the system when 
certain parameters, such as poles, zeros or gain, are 
adjusted. 

Through the transformation of the differential 
equations that describe a linear system into algebraic 
equations in the s-domain, a characteristic polynomial for 
that system can be obtained. In a feedback control systen, 
the coefficients of this polynomial are determined by the 
plant and any added compensators. The roots of the 


characteristic polynomial theoretically determine the 


system response in terms of bandwidth, settling time, 
steady-state accuracy, overshoot, etc. As the coefficients 
of the characteristic polynomial are varied, so too are the 
roots and the associated system response characteristics. 
Two parameters of the characteristic equation @ and f, are 
varied within user defined ranges to produce a plot with 
these two parameters representing the abscissa and 
ordinate. In essence, the characteristic polynomial acts as 
a mapping function whereby s-plane contours are mapped onto 
the a-—f plane. One can plot a family of parameter plane 
curves for various constant values of ¢, W,, ¢W, and/or o.” By 
selecting a desired operating point using these curves, the 
associated a and fB values can be graphically determined. 
The parameter plane plot provides the user with a visual 
means of deducing how the dominant roots of the 
characteristic equation move about in the s-plane as the 


values of aand Pf are varied. 


*In a controls’ systen, ¢ is the relative damping 
coefficient, Ww, is the undamped natural frequency and o is 
the real part of a root of the system's characteristic 
equation. The product ¢w, provides an indication of the 
settling time of a system, where the settling time is 
commonly equated to 4/C€w,. All four of these parameters can 
be precisely determined given a specific characteristic 
equation root value. Chapter 2 discusses these parameters 
and Figure 2-2 depicts a root on the s-plane. 


B. PURPOSE 

There currently exists at the Naval Postgraduate School 
a parameter plane subprogram of a major controls systems 
program residing on the school's IBM 370 mainframe. While 
offering excellent graphics resolution coupled with the 
availability of other analysis tools, such as the root 
locus method, in the same controls analysis package, it 
lacks the portability and ready accessibility inherent with 
personal computer compatible programs. Current 
microprocessor capabilities justify the development of a 
parameter plane program for analysis on a= personal 
computer. 

The parameter plane routine developed for this thesis 
incorporates algorithms which were originally proposed by 
D. Mitrovic [Ref. 1] and expanded upon and made more 
versatile by D.D. Siljak [Ref. 2] and G.J. Thaler [Refs. 3 
and 4]. These basic algorithms were used by R.M. Nutting in 
the development of his thesis [Ref. 5] which included a 
parameter plane program supported by a mainframe computer. 
Over two decades had passed when D.M. Potter, taking 
advantage of improved computer codes and vastly superior 
processors and computer architecture, updated Nutting’s 
program as part of his thesis [Ref. 6] by coding in Fortran 
77, simplifying user/machine interactions and utilizing the 


DISSPLA graphics package available on the school’s 


mainframe. 


In addition to the portability and versatility gained 
through the use of microcomputer compatible programs, the 
parameter plane program developed in conjunction with this 
thesis improves user friendliness over existing similar 
programs. Improvements include the incorporation of menu- 
driven prompts which provide a user with multiple available 
options on a Single screen. This presentation permits rapid 
selection of desired functions and plots, including the 
selection of many commonly desired curve groups with one 
key stroke. It also allows a user to bypass options not 
needed on a particular analysis. In addition, the existing 
mainframe parameter plane routine requires an exceedingly 
long time to generate and plot the user selected constant 
curves on the a—f-plane. This problem is exacerbated when 
mainframe computer usage load is high. The parameter plane 
routine presented in this thesis generates and displays 
plots in significantly less time and with much greater 


flexibility than was previously possible. 


II. PARAMETER PLANE DEVELOPMENT 


A. HISTORY 

A linear system can be formulated as a single ordinary 
linear differential equation with constant coefficients. 
Through the application of the Laplace transform, this 
differential equation is conveyed from the time domain to 
the frequency or s-domain and consequently can . be 
manipulated using algebraic techniques. The solution to the 
algebraic problem of synthesizing a control system was 
published in a paper by Dusan Mitrovic in 1959 [Ref 1]. 

Mitrovic's method operates in terms of both the 
frequency and time domains. It is an analysis and design 
technique of linear feedback control systems which applies 
graphical methods to algebraic equations. Mitrovic's 
procedure is based upon conformal mapping from the s-plane 
to the coefficient (a—f) plane through the characteristic 
equation. This methodology transforms an s-plane 
presentation into the real domain. The real domain in this 
case is defined by a coefficient plane whose coordinate 
axeS are two parameters («a and f) of the characteristic 
equation. These parameters appear in the coefficients of 
the characteristic polynomial. 

In the masters theses of H.H. Chon [Ref. 7} and C.H. 


Hyon [Ref. 8], Mitrovic's method was applied to a variety 


of linear feedback control systems thus illustrating the 
many useful applications of this technique. However, a 
major limitation of Mitrovic's method exists in that only 
two coefficients of the characteristic polynomial can be 
considered as variables. In real world systems, the desired 
adjustable system parameters, such as system gain and a 
pole of a compensator, are often located in more than two 
coefficients fone the characteristic polynomial. 
Straightforward application of Mitrovic's technique could 
not be accomplished in such cases. 

Recognizing this limitation, D.D. Siljak expanded upon 
Mitrovic's work by introducing Chebyshev functions in the 
graphical procedure of the method in 1964 [(Ref. 2:pp. 451- 
453]. This addition simplified Mitrovic's procedure and 


made it more convenient for computer simulation. 


B. BASIC ALGEBRAIC RELATIONSHIPS 

The characteristic equation of a feedback control 
system is simply the denominator of the closed loop 
transfer function of that system. Given a simple unity 
feedback system as presented in Figure (2-1), the closed 
loop transfer function is defined as: 


C(s) N(s) p(s) 


Ris) Ns) + Dis) sls) (a 











Figure 2-1 
Simple Unity Feedback System 


p(s) and /fls) are polynomials in the s-domain. The 
characteristic equation of this system is f(s) and is of the 


form: 


2 


fis) = s"ta_ys™' +4,557 7 + .. + ays? + ays + a ian 


This form occurs after dividing all terms of the equation 
by the highest order coefficient (coefficient of = 

Mitrovic illustrates the transformation of a 
characteristic polynomial in the s-domain into one which is 
a function of w,(the undamped natural frequency) and ((the 
relative damping coefficient) by assigning: 


ets 


A e) . 


—w, sind + jo,cosd = —w,f + jo,/l — C0? (2-3) 


where 


IA 
a 
IA 


0O< 0< = and 0 


The original lowest order coefficients of the s-domain 


characteristic equation are then defined as: 


a, = a,$2(0)O, + a363(C)O,° +... + a,6,(6)o, n—| 


(2=4)) 

ae if wee) 
a = —w, La2?,(C) + a362,(C)w, +... + a,6,_,(C)o, (2-5) 
Functions ?,(¢) are fixed values regardless of the 


degree oor coefficient values of the characteristic 


equation. %,(S) is calculated by: 
O,(0) = —L26b,_1(0) + ,_2(0)J (226) 


with 


oo(f) = 0 and $,(¢) = —I 


A plot can be constructed, with parameters @g and a, as 
coordinate axes, by calculating the values of ap, and a, as 
Cand w, are varied over a user defined range. This method 
obviously lends itself well to digital computer solution. 

Analysis of stability and system compensation are 
important applications with this technique. Through the 
selection of particular values for ag and a, corresponding 
to desired ¢ and/or w, values, the transient response of 
the system can be molded. 

The ¢ equals zero curve corresponds to the imaginary 
axis of the s-plane. This is the critical line determining 


system stability. By examining the ¢ equals zero curve on 


the parameter plane defined by a, and aj, values can be 
determined that yield a stable systen. 

As previously stated, Siljak introduced Chebyshev 
functions in the graphical procedure of Mitrovic's method. 
As did Mitrovic, he expressed s in terms of ¢ and wo, as 
depicted in Equation (2-3). The relationship of 6 and ®, in 


the s-plane is illustrated in Figure 2-2. 


Real 5 





Figure 2-2 
S-Plane 


s= -o,C + jw, /1 - Ee 


He then applied Chebyshev functions T,{¢) and U,{¢) to 


transform Equation (2-3) into: 


* = 0, [h(-0) + NICU (-9) (2-7) 

where 
Ty(—0) = (-1)°T, ©) (2-8) 
Uf —o) = (-1)°7'U,{0) (2-9) 


The argument ¢ of the Chebyshev functions is restricted to 
0<C€<1 for stable systems. 
The functions 7,(¢€) and U,(¢) can be calculated through 


application of the recurrence formulas: 


je 9 9 ae er (2-10) 
Ugay(S) — 2CU,(S) + Op (S) = 0 (25 iy 
where 
ion =a 
iiG)e =sec 
Ulf) = 0 
Uf) = 1 


Alternately, the Chebyshev function values can be 


computed using trigonometric functions as follows: 


T,(¢) = cos(K arccos C) (2-12) 


sin(k arccos ¢) 


O,(G) = sin(arccos ¢) (2-13) 
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Siljak then substitutes the value of s derived in 
Equation (2-7) into the equation for a characteristic 
polynomial as shown in Equation (2-2). By applying the 
condition that the summation of the real terms and the 
imaginary terms must go to zero independently, the 
characteristic polynomial can be rewritten as two 


simultaneous equations: 


~ k 0 
ay ads T,< a C) = (2-14) 
Xan Ui(—0) = 0 (2-15) 


The function 7,({) may be expressed in terms of U,(¢€) as 


follows: 


PG) = 60, (6) = O,.3(0) (Be) 


Redefining 7;(—¢) and U,(—¢) in terms of Equations (2-8) 


and (2-9) yields: 


~.( -1)'a,0,°Uz(-0) = 0 (2-17) 
r ( —1)"a,w, “U; aC) = 0 ae) 
A= 


The coefficients of the characteristic polynomial can 


be divided into «, and constant terms in this manner: 


a, = ba + eB + dy (2-19) 


ab 


Substituting the preceding equation into Equations (2- 


17) and (2-18) results in the following simultaneous 


equations: 
aBi(w,,6) + BC,(w,,6) + D,(w,,6) = 0 (2-20) 
2B(Wp,l) + BC2(w,,6) + D2(w,,6) = 0 (eae 
where 
a k k 
B, aa x (-1) 5, , Ui, 
k=0 
m . ° 
By = > (-1)"bw, “U; 
k=0 
m 
C= Eiger vm 
k=0 
m - r 
Ca ¥ (-1)* cy, U; 
A=0 
— \k k 
D, = X(-1)'4,w, U;_ 
k=0 
D, = (=I ene 
2 = =)! dw, U;, (2-22) 


Cramer's rule can be applied to the solution of the two 


Simultaneous equations defined in Equations (2-20) and 


C2 — 210 )e: 
ae C,D, ad C,D, (2-23) 
B,C, — B,C, 
B = B,D, — B,D» 
BC, — B,C, (2-24) 


Holding ¢, w, or (w, constant and solving the preceding 


equations while varying the other parameters yields a loci 
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of points corresponding to the roots of the characteristic 
equation with constant relative damping, undamped natural 
frequency or settling time. 

However, if adjustments of a and fB to achieve a 
settling time associated with a particular (w, value are of 


primary concern, Equation (2-3) should be rewritten as: 


sk = Pi(lw_, On7) + jog/l — 67 O,(la_, O,”) (2-25) 


The functions P, and Q, are related to the Chebyshev 


functions 7, and U,; as follows: 


Py (CWp, Wp 5) 


nT (-0) = (- 1)", “7,(0) (2-26) 
Oe oma a lu 0) = (1), Ur) — (2-27) 


The recurrence formulae associated with P, and QO, 
are: 
: 2 = 
Pras + 20w,P, + Wy Mate 0 (22S) 
4 
On) + 20w,Q;, + Wp hay = 0 (2-29) 
where 


iil 
= 


Po(Cap On) 
Pi(CWq On’) 
Oy(S0n, Wn) 
Q,(Lo_, On”) 


a Co, 


IN 
° 


HT 
= 
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As with 7, and U,, P, can be expressed in terms of Q,: 


Py = Ce — On Q-1 (2-30) 


In the same manner as Equations (2-17) and (2-18) were 


derived, Siljak produced two simultaneous equations in 


terms of Q;: 


5 4,0, =u (2 See 
=I) 
2 a,.Q, = 0 (2-32) 


Once again, as in Equation (2-19), a, is made up of z,f 
and constant terms. Therefore, the solutions to a and f in 
this case are identical to Equations (2-23) and (2-24). 


However, the expressions for Bj, Boj, Cy, C5, Dy and D5 now 


become: 

Wr 

B, = pa Ore 
k=0 
m 

B, = » b,x 
Ak=0 
m 

Co = LeyQy 
k=0 
mm 

CG = DC, Ox 
k=0 
7 

D, = Dd Qy_ 
k=0 
m 

D, = » 4,0, (2-337 
A=0 


As with functions 94,(¢), both Chebyshev functions are 
fixed values, not affected by differences in coefficient 
values or the degree of the characteristic polynomial. Once 
again, solution Of Mitrovic's method incorporating 
Chebyshev functions lends itself nicely to digital computer 


techniques. 
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III. PARAMETER PLANE PROGRAM DESCRIPTION 


A. BACKGROUND 
1. Introduction 

The parameter plane program is menu driven whenever 
possible. Many menus have an option to select other menus 
and these, in turn, may have that same feature. If the 
program user does not wish to examine particular curves, 
printer or labeling options, roots, etc., he is not subject 
to stepping through unwanted menus or responses. 

2. Software/Hardware 

The programming language used in coding this 
program is Microsoft FORTRAN77 V4.01. The Plotworks PLOT88 
graphics library is used to generate output plots. 

All simulations conducted for this thesis were 
performed on IBM-AT or IBM compatible 80286 machines. Due 
to the universality of FORTRAN coding, this program could 
be implemented on any machine capable of being programmed 
in FORTRAN. The source code for the parameter plane 
program is listed in Appendix. 

A note before eaea ae the program. ANSI.SYS should 
be incorporated in the personal computer's CONFIG.SYS file 
prior to running the parameter plane program. Through the 
ANSI.SYS device driver, system calls to clear the screen 


and position the cursor are enabled. If not, screen 
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readability suffers although the program is still fully 
functional. More will be said on this in the ANSI Module 
section of this chapter. 

3. Parameter Plane Program 

The parameter plane program itself consists of less 
than a page of code. Other than assigning default settings 
to certain variables used throughout the program, its sole 
function is to serve as a central switchboard to route 
calls to the eleven major subroutines which make up the 
parameter plane package. 

There is a minor subroutine which is only called 
when the program is initiated. This subroutine presents an 
introductory menu entitled 'LOAD/INSERT MENU' (Figure 3-1). 
It provides options to load a problem from an existing 
file, input the characteristic polynomial of a system the 
user wishes to examine, view an example of how to input a 
characteristic polynomial or quit the program. A _ flag 
corresponding to the selected option is set, control is 
returned to the main program and the selected subroutine is 


called and executed. 
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1 LOAD Problem from File | 
2 | INPUT Characteristic Equation | 
3 EXAMPLE Characteristic Equation Input | 
9 | EXIT to Main Menu | 


Enter integer number for selection ===> 


Figure 3-1 
LOAD/INSERT Menu 


Upon completion oof loading a characteristic 
polynomial either manually or through a previously created 
file, a subroutine providing the user with and titled MAIN 
Menu is called (Figure 3-2). From here, subroutines 
encompassing all user available options can be selected. 


These subroutines are grouped into four broad categories. 


af | CURVE Selection Menu | 
| 2 | PRINTER Selection Menu | 
3 | REVIEW/CHANGE Selections l 
4 LOAD Problem from File | 
l 5 l SAVE This Problem | 
6 | INPUT Characteristic Equation | 


| 7 | PLOT Curves 
8 ROOT Finder 

RE EERE EEE EER EEE FREER EEE REEF EEE EEE EE FE FETE PEt ttt ett 
| 9 | EXIT Program | 


Enter integer number for selection ===> 


Figure 3-2 
MAIN Menu 


There are four primary modules associated with the 


parameter plane routine: a user utilities module, a curve 
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selection module, a plotting module and a root finding 
module. In addition, there is a small module consisting of 
two subroutines; one to clear the screen and the other to 
position the cursor on the monitor screen. Figure 3-3 is a 
schematic illustrating the organization and 
interrelationships of the modules and associated 


subroutines of the parameter plane package. 






PARAPLAN 
alll re 
INTRMN MAINMN 





USER CURVE 


ROOT 
UTILITIES SELECTION 


FINDING 


EDEILE SAVFIL REVCHG CRVSEL RUNPRO ROOTS 


MONPRT |} | CHAREQ |} SAMPLE ZETA FEOT ROOTS] 
OMEGAN ZETAOM NORMAL 


Figure 3-3 
Parameter Plane Schematic 


B. USER UTILITIES MODULE 
1. Overview 
The user utilities module contains six subroutines 
which provide primary call routing procedures, data 
Manipulation and review and output setup. In addition, a 
sample problem subroutine is incorporated within this 


module. 


ig 


2. Subroutine CHAREO 


The provision eee manually inputting a 
characteristic polynomial is provided through Subroutine 
CHAREQ. The initial prompt displayed upon calling this 
subroutine requests entry of the order of the 
characteristic equation. A ninth-order polynomial is the 
maximum size limitation. 

Chapter II goes into great detail explaining the 
development of the parameter plane method. Equation 2-19 
defines the algebraic form of each coefficient of the 
characteristic polynomial. One can see that each 
coefficient term is composed of three parts, in this case 
defined as an alpha, a beta and a constant part. Alpha and 
beta are the two user defined parameters which the designer 
wishes to set to obtain certain desired system response 
characteristics. A coefficient containing only a constant 
part is not affected when system characteristics such as 
relative damping or undamped natural frequency are altered. 
These alterations will of course change the values for 
alpha and/or beta however. Subsection 6 in this section 
titled Subroutine SAMPLE provides a simple example on the 
mechanics of assigning values to the constant, alpha and 
beta coefficient terms. Chapter IV illustrates solving more 
complex problems with the parameter plane method. 

Following input of system order, requests are made 


to separately enter the three parts of each characteristic 
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polynomial coefficient. Constant coefficient terms are 
entered first. Input is requested from highest order term 
to lowest. For example, in the case of a third order 
polynomial, the first prompt would be ‘CONSTANT Coefficient 
of S ** 3 = ' with the cursor positioned just after the 
equals sign, awaiting user input. In many instances, a 
polynomial will not contain constant terms for certain 
coefficients. In these cases, the user must enter a zero 
for the applicable term(s). When insertion of constant 
coefficients is complete, they are echoed back and an offer 
to change incorrectly entered values is provided. 

After satisfactory entry of the constant 
coefficients, the user receives a prompt to enter the alpha 
coefficient values of the characteristic equation. Once 
again, these terms apply to the first of two user selected 
system parameters corresponding to a system or compensator 
gain, pole, etc. which the user desires to fix in order to 
achieve a particular set of system characteristics. 
Chapters I and II discuss the background and development of 
this method while Chapter IV contains a number of problems 
which illustrate its application. The same procedures used 
to enter and verify constant coefficients are used in the 
alpha value inputs. The user is then prompted to enter the 
coefficients for the second of two parameters, the beta 


terms. 


Pat 


When entry of all constant, alpha and beta 
coefficients is complete, a prompt is made for input of the 
minimum and maximum undamped natural frequency values over 
which to calculate all requested relative damping 
coefficients. Subsection 6 provides an example on 
determining these values in discussing Subroutine SAMPLE. 

A final function of Subroutine CHAREQ, in 
conjunction with the root finding module, allows one to 
enter specific values for alpha and beta and then calculate 
the roots of the existing characteristic equation. 

3. Subroutine LDFILE 

Subroutine LDFILE permits program input of a 
preexisting data file. The user is queried as to file name 
(not to exceed eight characters) and file extension (three 
characters or less). The code checks for existence of the 
file within the working subdirectory. If it does not exist, 
a message is returned stating such and offering the option 
of entering another file name or returning to the MAIN 
Menu. Entry of file extension is optional. In fact, a file 
can be labeled with as little as one character or number. 

The data file contains graphics output options (to 
the monitor or a particular printer type), minimum and 
maximum frequencies over which to span constant damping 
coefficient selections, the order of the characteristic 
polynomial and associated coefficient values. When loaded, 


all parameter values are displayed on the screen for 
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verification. It is most easily created using MAIN Menu 
option 5, ‘Save this problem', after an analysis has been 
run. This option calls Subroutine SAVFIL to be discussed 
next. A data file could be created using a text editor and 
aligning the previously stated parameter values in the 
format identified by the Subroutine LDFILE source code 
listed in Appendix A. Of course the first method is quicker 
and less subject to errors. 
4. Subroutine SAVFIL 

After entering in Subroutine CHAREQ the previously 
discussed parameters which define a system and the 
particular range of frequencies of interest, it would be 
convenient and time efficient to save the system definitive 
numbers, if one is to look at that same system in the 
future. Subroutine SAVFIL provides this service. As with 
Subroutine LDFILE, a prompt is displayed requesting input 
of file name and file extension. In addition, a check is 
made to determine if the entered file name plus optional 
extension already exists. If it does, the option is 
provided to overwrite the existing file or enter a new file 
name and extension. The current parameters as outlined in 
the Subroutine LDFILE subsection are then written to the 
indicated file. All parameters are also echoed to the 
screen as a final check for accuracy. 

Even high order polynomials require only 300 to 


400 bytes of memory, therefore storage of multiple systems 
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will have little impact on disk space used. Approximately 
10,000 high order systems could be stored on a double 
sided, double density diskette. 
5. Subroutine REVCHG 
Occasionally, a need arises to review the system 
coefficient values and frequency range under consideration. 
This can be done by loading or saving a file or accessing 
Subroutine FREQ in the Curve Selection Module. In most 
cases, this procedure can be accomplished more quickly by 
selecting option 3, 'REVIEW/CHANGE Selection', in the MAIN 
Menu. This selection calls Subroutine REVCHG. As indicated 
by the option title, changes to system parameters can be 
made directly from this subroutine. 
6. Subroutine SAMPLE 
Subroutine SAMPLE provides a quick tutorial on the 
derivation and entry of a system's characteristic 
polynomial. As an example, a universal third order system 
is presented. 
A universal third order system is represented by 


the characteristic equation s? + As? 


+ BS + 1= 0. A and B 
are the variables defining the two system/compensator 
parameters to be determined so as to achieve desired system 
response. This characteristic equation is obtained with a 
*) 


l/s? plant incorporating both acceleration (s and 


velocity (s) feedback. A and B (the alpha and beta terms 
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identified in Equation 2-19) are the respective gains for 
the two compensators. 

To load this characteristic polynomial into the 
parameter plane model, one would first select option 2 in 
the introductory LOAD/INSET Menu or option 6 in the MAIN 
Menu. Each of these option titles is ‘INPUT Characteristic 
Equation'. 

A simple examination of the characteristic 


2 


polynomial shows that the s“ and s terms have no constant 


3 and constant (1) terms have no 


coefficients while the s 
alpha (A) or beta (B) terms. Therefore, the response to 
‘Enter the CONSTANT Coefficient Values of the 
Characteristic Equation' would be: 
"CONSTANT Coefficient of s**3 =! 
"CONSTANT Coefficient of s**2 =! 


‘CONSTANT Coefficient of s**1 =! 
‘CONSTANT Coefficient of s**0O =! 


IR |O IO |e 


The response to the alpha and beta coefficient request is: 


"ALPHA Coefficient of s**3 =! 
'ALPHA Coefficient of s**2 =! 
‘ALPHA Coefficient of s**1 =! 
"ALPHA Coefficient of s**0 =! 


lo fo lr lo 


"BETA Coefficient of s**3 =! 
'BETA Coefficient of s**2 =! 
'BETA Coefficient of s**1 =! 
'BETA Coefficient of s**0O =' 


[SO |r fo fo 


Entry of minimum and maximum values of undamped 
natural frequency would be dictated by the desired system 


Meeponse Characteristics. For example, a settling time (TQ) 


Zo 


under 2 seconds is required. System settling time is 
approximated by four divided by the product of relative 
damping coefficient (zeta) and natural undamped frequency 


(omega Since zeta is restricted to values between zero 


n)- 
and one, omega, ranges can be calculated (using minimum 
value of .1 for zeta). The result is»that omega jwicanmary 
between 2 and 20, therefore enter these values as minimum 
and maximum frequency. 

7. Subroutine MONPRT 

The largest of the subroutines within the User 
Utilities Module is Subroutine MONPRT. A number of smaller 
subroutines are associated with it. Subroutine MONPRT 
allows selection of the graphics output device, be it the 
monitor or a wide variety of printers. Graphics output 
defaults to the monitor when the program is first executed. 
Graphics output can always be dumped from the screen to a 
printer, but there may be times when direct output to a 
printer is desired. All NPS printers are included in the 
printer output options. In addition, one option permits the 
user to directly enter the values for IOPORT and MODEL, as 
outlined in the PLOT88 Manual [Ref. 10]. 

The first menu displayed when Subroutine MONPRT is 
called is the PRINTER/OUTPUT Menu (Figure 3-4). Selection 
of a printer in this menu automatically directs graphics 
output to the most commonly associated output port, either 


parallel or serial, for that particular device. However, 
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usually there are multiple parallel or serial ports to 
which a printer can be attached. The default printer output 
setting may not direct plots to the appropriate port. To 
provide maximum output flexibility, the PRINTER/OUTPUT Menu 
allows selection of a specific output port through access 


to another menu. 


| PRINTER NO. | PRINTER 
| al | Epson FX-80, All | 
| 2 | Epson FX-100, All | 
| 3 | Epson MxX-100, All | 
| 4 | Epson RX-80, All | 
| 5 | Epson MX-80 & IBM Printer | 
| 6 | HP 7470A Graphics Plotter | 
| 7 | HP 7475A Graphics Plotter | 
| 8 | HP 758xB Series Plotters | 
| 9 HP 2686A Laser Jet | 
TEPER ETHERS FRPP RE REAP HEAP RELA APEEFEEEE PH E PPE EHH Ett tte 
| 10 Graphics Monitor (default) | 
| algk HARDWARE Interface Menu | 
| | 
| | 


: 
a2 | Input PLOT88 Values for IOPORT and MODEL 
| EXIT to Main Menu 


eee nee Sse SOS SS SOP SP See Se SSP SS SP SP SP SP SP SP SP SP SF SF S88 SF SF SF SSF SF 8 SF SVS SF SF Ss SO €VB Ss Ss Ss Ss Ss £€S SF SF S88 SF SS S&S &S &S = = 


Enter integer number for selection ===> 


Figure 3-4 
PRINTER/OUTPUT Menu 
The output port selection option is called the 
HARDWARE Interface Menu (Figure 3-5). This menu is part of 
a subordinate MONPRT subroutine titled PORT. The options 
allow selection of one of three parallel ports (LPT1 - 


LPT3) or either of two serial ports (COM1 - COM2). 
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| | LPT1 Printer Port | 
| 2 | LPT2 Printer Port | 
| LPT3 Printer Port | 
| COM1 Serial Port | 
| 5 COM2 Serial Port 
+ ~ 
| | 


FEEAEEEE EE EEEEEFEFEEHEEFHEFEEFEFEE FEE FFEF+ FEE + 4444+ 4+ 44444 4+ 4444444 
99 EXIT to Main Menu 


Enter integer number for selection ===> 
Figure 3-5 
HARDWARE INTERFACE Menu 
If a serial port is manually selected in the 
HARDWARE INTERFACE Menu, an associated data transfer (baud) 
rate must also be assigned. Immediately after selection of 
serial port options 4 or 5, the BAUD (data transfer) RATE 
Menu is displayed (Figure 3-6). Transfer of graphics data 
over a Single line is usually time consuming, therefore the 
highest printer capable transfer rate available (9600 baud) 
is normally the best selection. A slower transfer rate 
would only be indicated if improper graphics output is 


generated. 
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ll ‘SSELECGRENO.. | BAUD RATE | 
| 1 | 300 l 
| 2 | 1200 | 
| 3 | 4800 | 
| 4 | 9600 | 


Enter integer number for selection ===> 
Figure 3-6 
BAUD (data transfer) RATE Menu 
Bit by bit transfer of data can also have a check 
Sum (parity) associated with it to provide an internal 
check of correct data transfer. Following selection of a 
data transfer rate, the PARITY Menu is displayed. The user 


has the option of choosing odd, even or no parity. 


fn SE LE Clh NO: | PARITY | 
| “ | NO Parity | 
| 2 | EVEN Parity | 
| 3 | ODD Parity | 


Enter integer number for selection ===> 
Figure 3-7 
PARITY Menu 
The PLOT88 Software Library Reference Manual [Ref. 
10} contains a number of tables which provide common output 
settings for many more printers than can be selected in the 


PRINTER/OUTPUT Menu. 
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C. CURVE SELECTION MODULE 
1. Overview 

The Curve Selection Module contains only one 
subroutine. However, Subroutine CRVSEL is the single 
largest subroutine in the entire program and contains all 
of the coding that makes this program a parameter plane 
analysis tool. Attempts were made to subdivide this large 
subroutine into many function specific smaller ones. 
Apparent Fortran compiler limitations due to the overall 
Size of the program frustrated this effort. 

Numerous menus in Subroutine CRVSEL permit the user 
to rapidly select a variety of data display and file 
options as well as commonly desired constant parameter 
curves such as the zeta equals zero curve, defining the 
system stability limit. These menus also allow one to 
easily determine the curve selection process level and 
provide the ability to backtrack if necessary. 

Three constant parameter choices are offered for 
solution of alpha and beta values associated with a 
particular root position and subsequent plotting on the 
alpha/beta plane. The first of these choices is the 
constant zeta contour. Zeta 1s the variable normally 
assigned to the relative damping coefficient. It is 
represented on the s-plane aS a radial vector extending 
outward from the origin and is directly related to system 


steady state overshoot. The zeta equals zero curve is the 
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imaginary axis of the s-plane while the zeta equals one 
curve falls on the real axis. The second contour available 
for plotting is the constant omega, curve representing 
system undamped natural frequency. Its s-plane 
representation is a curve of constant magnitude about the 
Origin. The intersection of these two contours exactly 
defines a system root location. The final curve is the 
product of the first two contours and indicates the 
transient response of the system. This curve plots as a 
straight line on the s-plane parallel to the imaginary 
axis. System settling time is usually defined as four 
divided by the zeta-omega, product. The interrelationship 
of these terms with the s-plane is depicted in Figure 2-2. 
2. Subroutine CRVSEL 

Subroutine CRVSEL is called from the main program 
through selection of option one in the MAIN Menu, labeled 
CURVE Selection Menu. Upon selection of this option, the 
CURVE DATA POINT DISPLAY Menu appears on the screen (Figure 
3-8). Option selection from this menu determines if the 
user is provided with a screen display of every tenth 
alpha/beta pair and/or all alpha/beta solutions are routed 
to a data file for future reference. Another available 
option provides neither a screen display of alpha/beta 
pairs nor a dump of these data points to file. Selection of 
this option speeds the curve construction process slightly 


by avoiding I/O and may be the preferred option if only a 
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quick graphical overview is desired. The CURVE DATA POINT 
DISPLAY Menu also provides the option of computing system 
closed loop root values associated with each alpha/beta 
pair. These roots can be displayed on the screen and/or 


written to a file. 


| OPTION NO. | OUTPUT SELECTION | 


| | NO output DISPLAY or save to FILE | 
| 2 | DISPLAY every 10th alpha/beta value | 
| 3 | DISPLAY alpha/beta values and roots | 
| 4 | DISPLAY & FILE alpha/betas and roots | 
| 5 | FILE all alpha/beta values | 
| 6 | FILE all alpha/beta values and roots | 


woe weew ew me ew www ww ew www ww wm ewe wee ene ewe rrr ererrerer**r>>r>>"—"EeQ———~ 
see = 


Enter integer number for selection ===> 
Figure 3-8 
CURVE DATA POINT DISPLAY Menu 

When options 4,5 or 6 are made in the CURVE DATA 
POINT DISPLAY Menu, indicating that calculated data is to 
be filed, the user is queried as to name and extension of 
the data file to be opened. A check is made to determine if 
the entered data file name and extension already exists. If 
it does, the user is given a choice to overwrite that file 
with newly generated data or open a new file under another 
file name. 

When the data display and storage procedures just 
Giscussed are completed, the primary menu of this 
subroutine is displayed. It is called the CURVE SELECTION 


Menu (Figure 3-9) for obvious reasons. In addition to 
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options allowing selection of any one of the three constant 
parameter curves previously outlined, options to change the 
frequency range over which constant zeta contours are 
calculated, plot the selected curves or return to the MAIN 


Menu are offered. 


m_en an oe a a S62 a oe SP SP oe SS SS ca ae ae Se ee ee ie ee a i i es es Se ee SE GS es eas OS OSC SE EE 


| OPTION NO. | CURVE SELECTION | 
i | Constant ZETA Curves | 
| 2 | Constant OMEGA Curves | 
| 3 | Constant ZETA*OMEGA Curves 
| 4 | Change Frequency Range | 
FEAF AEEEPEFEPFAFEFEFRAFAFEFHEHFEFEFE FEAF EFEFEHETHHE HHH +H tet tee tet 
5 | PLOT Selected Curves | 
FAFA EEEEPEHALEEEFEEEEFEEEEEEEFEFEFEEE EEF TEE EE EE E+E HE +H +++ $4 44444 
9 | EXIT to Main Menu | 


Enter integer number for selection ===> 


Figure 3-9 
CURVE SELECTION Menu 

Each of the three curve options permits’ the 
selection and plotting of up to ten constant contours. In 
most instances, a designer is interested in viewing the 
stability limits of a system, as represented by the zeta 
equals zero curve. Due to this zeta curve attribute and the 
fact that for a stable system, zeta is restricted to the 
range zero to one, a CONSTANT ZETA Menu (Figure 3-10) is 
provided. The user is given the option of entering 
particular constant zeta curve values or selecting with a 
Single key stroke the zeta equals zero curve or two other 
options which equally subdivide the range zero to one into 


three or five parts. Following identification of constant 
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zeta curves, the program returns to the CURVE SELECTION 
Menu to allow selection of other constant curves or plot 


the current selections. 


| OPTION NO. | ZETA CURVE SELECTION | 
| ak | Select particular constant ZETA curves | 
| 2 | ZETA = 0 curve | 
| E | | ZETA = 0,0.5 and 1.0 curves 
| 4 | ZETA =f0,0.25,0.5,0.75,1.0 Curves | 
a oe a a ee eo Se oe a 
9 | EXIT to Curve Selection Menu | 


Enter integer number for selection ===> 


Figure 3-10 
CONSTANT ZETA Menu 

Fach constant curve is calculated using 100 data 
exepbghess,  Uble\ale value was selected as a compromise between 
speed of curve generation and precision of plot. A possible 
update to this program would permit a user entered data 
point number. This would enable a designer to more closely 
tailor the parameter plane plots to his or her needs. 

Successive data points are linearly incremented in 
the case of constant omega, and zeta-omega,, curves. In the 
generation of successive constant zeta points, a 
logarithmic increment is used due to the large possible 
variation between minimum and maximum omega, values. (In 
calculating alpha/beta pairs associated with zeta curves, 
omega, is varied between minimum and maximum while zeta is 


held constant). 
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High order systems often exhibit large variations 
between successive alpha/beta points. Unfortunately, the 
PLOT88 graphics package occasionally fails when these 
variations are extreme (i.e., in excess of five orders of 
magnitude). To prevent inadvertently halting program 
execution, two checks are made. The first prevents 
calculation of an alpha or beta value when the denominator 
of the defining equation is less than 1.e-12. The second 
compares successive alpha/beta points for a difference 
exceeding 1.e5. If either situation occurs, a warning is 
sent to the monitor and, if applicable, to a data file. The 
offending alpha/beta pair is then assigned the value of the 
previously calculated point. 

A final procedure available in Subroutine CRVSEL is 
the selection of minimum and maximum frequencies over which 
constant zeta curves are calculated. Should this option be 
exercised after selection of constant zeta curves, the 
program automatically recalculates the associated 


alpha/beta points. 


D. PLOTTING MODULE 
1. Overview 
The Plotting Module consists of three subroutines. 
The primary function of this module is to call PLOT88 
Graphics Library Routines which generate the constant curve 


plots on the alpha-beta plane. The alpha/beta arrays are 
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passed to this module via a common block. This method was 
adopted when attempts to pass large arrays as arguments of 
subroutines resulted in stack overflow errors. 

As discussed early in this chapter, the latest 
released Fortran version was used in compiling and linking 
this program. Microsoft FORTRAN77 v4.01 contains an 
optimizing compiler that does not directly translate the 
source code into an object file but instead alters the code 
to achieve various user defined objectives. Two common 
optimizations are increased speed and reduced code size. 
The Plotting Module holds the distinction as the only 
module that could not take advantage of this optimization. 
One or more calls to PLOT88 Routines were distorted to the 
point that program execution was halted unpredictably. 
Therefore, the optimization functions of the Fortran 
compiler were bypassed. 

2. Subroutine RUNPRO 

Subroutine RUNPRO makes up the bulk of this module. 
When first called, the PLOTTING Menu (Figure 3-11) is 
displayed. This menu provides the user with a number of 


plot presentation and labeling options. 
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| 2 | TITLE output graph and PLOT data | 
| Z | PLOT data (no title) | 
| 3 | SIZE output graph | 
| 4 | SYMBOL to be plotted at each data point | 
| 9 | EXIT to Main Menu | 


Enter integer number for selection ===> 


Figure 3-11 
PLOTTING Menu 

Option 1 permits entry of a title up to 30 
characters in length. A routine computes the length of the 
entered title then centers it at the top of each graph. The 
third option permits adjustment of plot size to both the 
screen and the printer. This option was utilized to scale 
plot outputs for inclusion in this thesis. Another option 
allows the user to select a particular symbol to be applied 
to the plot at each data point. Figure 3-12 illustrates 


some of the symbol choices available. 


Type in an INTEGER number from 0 through 13 to place a 
symbol at each calculated data point. There are 100 data points 
for each curve. Examples of symbols with associated numbers: 


2 Triangle 
3 + 

4 xX 
8 Z 
9 Y 
1i * 
13 Vertical Line 


Enter INTEGER number (0 - 13) ==> 


Figure 3-12 
Symbol Selection Display 
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A final option permits immediate plotting of data 
arrays without enhancing the graphical presentation. As 
with most of the menus, the ability to return to the 
calling routine, in this case the MAIN Menu, is also 
available. 

Once again, three sets of system parameter curves 
can be plotted; constant zeta, constant omega, and constant 
zeta-omega,. Each set can contain up to ten separate 
contours. The program sequentially checks for the existence 
of any curves in each of the three categories. In addition 
to alpha/beta values, total number of curves and data 
points within each category are passed to Subroutine RUNPRO 
from Subroutine CRVSEL. When one or more curves are 
detected in a particular category, Subroutine PLTCRV is 
called to actually compute each curve. Requested curves in 
all categories are first calculated before the first plot 
is displayed. Sequencing to succeeding plots is 
accomplished by depressing the Enter key. All zeta contours 
are drawn as a series of straight line segments one tenth 
of an inch long. This procedure gives the appearance of a 
smooth, continuous contour. Omega, and zeta-omega,, 
contours are plotted similarly, but with dashed lines made 
up of different segment lengths. 

If both zeta and omega, curves are present, the 
last plot consists of both sets of curves. The minimum 


alpha and beta values and axes increments from the zeta 
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plot are used to scale this graph. If any omega, constant 
contour selection frequencies fall outside the minimum to 
maximum frequency range designated for constructing the 
constant zeta curves, these curves may not be displayed 
within the plotting box. The user can still view the 
offending contour lines by exiting then reentering the 
plotting routine and selecting the plot size reduction 
option in the PLOTTING Menu. 

After viewing the requested plots, the user is 
given the option of focusing on a particular area of the 
plot. This is accomplished through the EXPAND PLOT Menu 
(Figure 3-13). 


| 1 | Expand area defined by axes values | 
| 2 Expand around a selected point 

| ERE EAEEEAE ETHEL EEEALE EF FEFAEEFEFEPEFFE EEE FEEEA EE THAT++E+ +++ +++ | 
| 9 | EMEE plotting routine | 


Enter integer number for selection ===> 
Figure 3-13 
EXPAND PLOT Menu 
A choice is offered to expand a specified area 
Gefined by minimum and maximum alpha and beta axes values 
or expand the plot around a selected point relative to 
absolute alpha/beta axis positions (e.g., expanding about 
the center of the plot would be accomplished by entering a 


value of 3.5 for alpha and 2.5 for beta). The second option 
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provides the user with a quick look at a particular point 
however the axes are usually lost with any significant 
expansion. The first option rescales and displays the alpha 
and beta axes, therefore it is normally the preferred 
expansion method. 

A few remarks on the mechanics of producing hard 
copy plots, as they apply to Naval Postgraduate School 
specific hardware, are in order. Graphics output to a 
device other than a monitor requires an extended period of 
time as the data file is loaded into the printer buffer and 
formatted for hardcopy output. In addition, the user is 
unable to view the plot when the program has been 
configured for a non-monitor output device. To avoid these 
inconveniences, the ‘'Print Screen' key can be depressed 
while a desired plot is displayed on the monitor. Printer 
output is then available in less than a minute. The Laser 
Jet Printer located in the Controls Lab must be configured 
aS an Epson Printer for the 'Print Screen' function to 
operate correctly. 

3. Subroutine GRID 

Subroutine GRID is a short, simple routine called 
before each plot. It first defines the plot physical origin 
in relation to the output device. This is usually the 
monitor screen although it would be a printer page if so 
defined in Subroutine MONPRT. The origin is presently set 


at eight tenths of an inch above and to the right of the 
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output device lower left corner. It cannot be modified by 
the user. 

A box of unmodified dimensions 7" x 5" is drawn 
from the physical origin, defining the plotting area. In 
addition, a dashed line grid is inserted at integer 
intervals of the vertical and horizontal axes. The box, 
grid lines and contour plots are all size modified (both 
screen and printer output) through the sizing option of the 
PLOTTING Menu. 

4. Subroutine PLTCRV 

This subroutine accepts alpha/beta data arrays and 
graph title passed by subroutine argument. Common block 
inputs to size the output plot and place a user identified 
symbol at each data point are also resident. 

An initial call is made to a PLOT88 subroutine 
specifying axis annotation characteristics such as tic mark 
placement and axis title height. Next, axes scaling 
intervals are determined based on alpha and beta value 
ranges and a respective seven and five unit division of the 
axes. Additionally, the minimum alpha and beta data point 
values are identified at this time. The x and y axes 
minimum and major increment values are displayed on the 
screen with a message indicating that curve construction is 
in progress. Meanwhile, a PLOT88 subroutine is called which 
accepts the 100 data point pairs associated with each 


constant curve and creates a smoothed contour using the 
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scaling data previously calculated. Data point symbols are 
also inserted on each curve if selected by the user in 
Subroutine RUNPRO. When all contours have been created for 
a particular parameter (e.g., relative damping coefficient 


- zeta), control is passed back to Subroutine RUNPRO. 


E. ROOT FINDING MODULE 
1. Overview 
The Root Finding Module is unique among the three 
major modules of the parameter plane program in that it has 
no menus associated with it. It is called predominantly 
from the CRVSEL subroutine of the Main Module when a ‘'y' 
response is entered to the question 'Do you want to view 
the roots of the polynomial?' This question is presented to 
the user each time calculations for a set of constant 
curves is requested. 
The other option of entering the Root Finding 
Module is offered in the MAIN Menu as Root Finder. Should 
the user make this selection, a flag is set, Subroutine 
ROOTS is called and subsequently the subroutine for 
inputting a characteristic equation is activated. A choice 
is given to enter a characteristic equation or find the 
roots of the current characteristic equation. If the user 
wishes to enter a new or initial characteristic polynomial, 
the standard characteristic equation entry procedure is 


carried out. If not, this portion of Subroutine CHAREQ is 
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skipped over. In either case, entry of a characteristic 
equation via the root finding module requires that a user 
enter specific values for alpha and beta. All associated 
coefficients of these two variables are multiplied by the 
respective alpha or beta entries. Subroutine CHAREQ then 
returns the program to the root finding module for 
Calculation of the input polynomial with specified alpha 
and beta values. In this manner, one can select particular 
alpha and beta values and compute all roots of the 
polynomial with those values. 
2. Subroutine ROOTS 

ROOTS is the highest level subroutine of the root 
finding module. Terms defining the order and coefficient 
values of a polynomial are passed to ROOTS as parameters of 
the subroutine. Additional parameters passed include the 
particular alpha and beta values to be used in calculating 
the coefficients of the polynomial. Another parameter is 
the flag differentiating between calculation of roots 
associated with constant curves or the calculation of root 
values for specific alpha/beta user inputs. A call is made 
to Subroutine CHAREQ if this flag indicates that the roots 
of a polynomial associated with a user inputted set of 
alpha and beta values are to be calculated. 

Once the characteristic polynomial is fixed, 
Subroutine NORMAL is called. This subroutine normalizes the 


polynomial to ensure that the coefficient of the highest 
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order term is unity. Finally, Subroutine ROOTS2, which 
contains the root calculation procedures, is called. 
3. Subroutine NORMAL 

NORMAL is a subroutine that normalizes’ the 
characteristic polynomial. It is also the subroutine which 
multiplies the alpha and beta coefficient terms by the 
alpha and beta values passed via subroutine parameters from 
ROOTS. These two products represent a portion of the 
coefficient of a particular order term associated with the 
alpha and beta parameters. They are added to the constant 
coefficient of that same order term to form the total 
coefficient of that term. A do loop provides the means for 
combining the constant, alpha and beta coefficient values 
into a single coefficient for each term of the polynomial. 

When the roots of a polynomial are calculated based 
upon user input of specific alpha and beta values, 
Subroutine NORMAL provides another function. In addition to 
showing the user the roots of the polynomial, the routine 
also outputs the coefficient values of the polynomial based 
upon the specific alpha/beta inputs. 

4. Subroutine ROOTS1 

ROOTS1 is a short subroutine that models’ the 
quadratic equation. It calculates the roots of a second 
order polynomial. This subroutine is called only from the 
ROOTS2 subroutine after the initial polynomial has been 


reduced to one of second order. 
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5. Subroutine ROOTS2 

Subroutine ROOTS2 is the heart of the root finding 
module. The algorithm used in this subroutine iteratively 
searches for a quadratic factor of the input polynomial. 
When this factor is found, the original polynomial is 
reduced by an order of two. This process is repeated until 
the remainder is a first or second order polynomial. Each 
time a quadratic factor is found, Subroutine ROOTS1 is 
called to calculate the roots. 

There are two terms in the algorithm that affect 
both the speed of computation and the accuracy of the 
solution. The first, labeled EPSLON, sets an acceptable 
accuracy level for the computed roots. The numerical 
procedure applies ever decreasing correction factors to the 
constant and first order terms of the quadratic Aa Somer: 
is calculating. When the sum of these correction factors is 
less than an EPSLON value of .00005, the refinement is 
halted and the roots of the quadratic factor are solved 
through a call to the ROOTS1 subroutine. In the vast 
Majority of cases, the quadratic factor is determined after 
relatively few iterations. However, there are cases, 
especially when the characteristic polynomial contains 
multiple identical real roots, when the sum of the 
correction factors is not reduced below EPSLON even after 
1000 iterations. ithe this situation occurs during 


computation and the routine determines that more than 500 
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iterations have taken place in determining a_e single 
quadratic factor, the procedure is halted and control 
returned to the calling routine. Should this happen during 
the computation of roots for multiple alpha/beta values, 
the alpha and beta values will still be displayed without 
the associated roots and an attempt to determine roots for 
the succeeding discrete alpha/beta pair will commence. 

The methodology Of continually reducing a 
polynomial by a quadratic factor until all roots are 
determined is known as-_ Bairstow's method. R.L. Wood 
incorporated this method in his thesis when he coded a 
routine which included control systems design through root 


locus techniques [Ref. 9:pp. 42-43]. 


F. ANSI MODULE 

The ANSI Module contains two subroutines. Subroutine 
CLR clears the monitor screen. It is used primarily as a 
precursor to displaying a menu or positioning queries or 
informational text. Subroutine PSIT positions the cursor 
immediately after a request for response. This subroutine 
serves not only to give a pleasing screen display but also 
provide the user with an immediate reference as to where he 
is in the program and if an input is required of the user. 
Both routines can be found in the PLOT88 Reference Manual 


[Ref. 10:pp. 85-86]. 
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The two subroutines which comprise the ANSI Module make 
use of ANSI escape sequences. These sequences are a series 
of characters that can be used to define functions to MS- 
DOS. In addition to clearing the screen and positioning, 
moving and saving the cursor's position, ANSI escape 
sequences can control screen graphics and reassign key 
definitions. 

In order to utilize the escape sequences, the 
Microprocessor must first be able to interpret the 
commands. It does this through a systems file entitled 
ANSI.SYS. ANSI.SYS must be defined as a device in the 
CONFIG.SYS file to enable proper interpretation of ANSI 
escape sequences (e.g., device=\bin\ansi.sys where \bin\ is 
the path to ansi.sys). If ANSI.SYS is not resident in the 
CONFIG.SYS file, the parameter plane program will not clear 
the screen or position the cursor when so directed. 
Instead, the indicated escape code sequence will be printed 
on the monitor whenever a call is made to CLR or PSIT. 
Although program readability is degraded, parameter plane 


meee still function correctly. 
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IV. COMPUTER AIDED PARAMETER PLANE ANALYSIS 


A. OVERVIEW 

An algebraic solution of feedback compensated control 
systems is usually performed by selecting particular zeta 
and omegan values that provide desired system response 
characteristics. A pair of complex conjugate roots are thus 
placed at specific s-plane locations that will provide the 
designed for response, if these roots are dominant. A check 
must therefore be made to satisfy the dominance question. 
If the positioned roots are not dominant, the design scheme 
must be altered. 

Design procedures which select values of zeta and 
omegan to position a pair of complex roots ina position of 
dominance assume that the contributions of all additional 
roots can be neglected. Although this is certainly not 
always the case, applying this simplification permits 
modeling of complex systems with second order polynomials. 
As such, the following system response parameter equations 
can be utilized to derive zeta and omegan values’ which 
produce desired system behavior: 

Maximum Overshoot (step input) 


114 
M,, = ] pevi-2 (4-1) 
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Settling Time (4 time constants) 


4 
CWn (4-2) 





I; 


Time of Maximum Overshoot 


UT 


t = ee! 
‘ wrV1-6? ee 
Bandwidth (response magnitude > .707*input) 


Wny 1-207 + /2- 467 + 40 (4-4) 


Wp 


Resonance Peak 
W, = W,nV¥1-2C7 (4-5) 
Transient Oscillating Frequency 


W, _ Mn 1-167 (4-6) 


Phase Margin 


YD, = see. (a=) 
— 207 + /404+1 
A graphical solution using the parameter plane 


technique, with reference to the resident root finding 
module, provides a designer with a rapid means of 
determining parameter values satisfying design requirements 
while ensuring both stability and root dominance. When a 


family of curves is computed and displayed, the flexibility 
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of choosing a number of parameter value pairs which may 


force design compliance is offered. 


B. PROBLEM SOLUTIONS 
1. Example 4-1 

Figure 4-1 illustrates an ideal instrument servo 
which is to be damped by velocity feedback and a cascaded 
gain amplifier. The common single parameter root locus 
method could easily be employed in the solution of this 
problem by deleting the gain amplifier. However, this stage 
provides” the designer with increased flexibility and 
lllustrates the ease of parameter plane solution in low 


order systems, despite an additional parameter. 





Figure 4-1 
Ideal Instrument Servo with Velocity Feedback 


Design requirements specify: 

Maximum Overshoot Mpt (step input) < 20% 

Settling Time Ts < 2 sec 
Applying Equations (4-1) and (4-2) with design 


specifications to solve for zeta and omegan yields: 
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C > .46 


(4-8) 
W, > 4.35 rad/s (4-9) 
Let: 
a = Kk, (4-10) 
io = K (4-11) 
The characteristic equation for this system is then: 
s* + 500as + 5008 = 0 412) 


Run the parameter plane program by typing in the 
executable file name, PARAPLAN, when in the subdirectory 
under which it resides. The IBM microcomputers in the NPS 
Controls Lab are set up under the file management system 
1DIR. To execute the program under this system, enter the 
appropriate subdirectory by using the up/down arrow keys to 
position the cursor adjacent to the subdirectory name then 
depress the ENTER key. Use this same procedure to select 
PARAPLAN.EXE and push ENTER to run the program. 

The LOAD/INSERT Menu will appear. Select Option 2, 
INPUT Characteristic Equation. Enter coefficient array 


values after each prompt as follows: 


‘CONSTANT Coefficient of s**2 =’ ] 
‘CONSTANT Coefficient of s**1 =’ 6 
"CONSTANT Coefficient of s**0 =’ 0 
‘ALPHA Coefficient of s**2 =’ 0 
‘ALPHA Coefficient of s**1 =’ 500 
‘ALPHA Coefficient of s**0 =’ 0 
‘BETA Coefficient of s**2 =’ 0 
"BETA Coefficient of s**l =’ 0 
"BETA Coefficient of s**0 =’ 500 
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As this is the first example, let’s view the 
parameter plane over the entire range of zeta. After 
entering the coefficient array, the program sequences to 
MAIN Menu. Select Option 1, CURVE Selection Menu. Choose 
any of the curve data point display options as the next 
menu appears. Option 1 of the CURVE Selection Menu allows 
selection of constant zeta curves while Options 2 and 3 
provide this service for constant omegan and constant zeta- 
omegan curves. Input a range of zeta values from zero 
through one anda range of omegan values from one through 
ten. Input avalue of 2 for the zeta-omegan curve. Select 
Option 5 to plot the designated curves. The PLOTTING Menu 
now appears. Either Option 1 or 2 will start the constant 
curve plotting calculations and subsequently display the 
plots. The following graphs were obtained by selecting 


Option 1 and titling the curves. 
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Figure 4-2 
Constant Zeta Curves 
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Figure 4-3 
Constant Omegan Curves 
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Figure 4-4 
Constant Zeta-Omegan Curve 
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Figure 4-5 
Constant Combined Curves 

As expected due to the stable nature of this 
system, all curves are inthe first quadrant because all 
coefficients must be positive to provide roots in the left 
half of the s-plane. A stable system requires compliance 
with this condition, although the condition alone does not 
guarantee stability. 

At this point, we could select an option from the 
expansion menu to better determine the corresponding alpha 


and beta values satisfying our calculated zeta and omega 
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— Pa a a 


parameters. However, we can more precisely determine these 
values by plotting the actual desired zeta and omegan 
curves. 

Exit from the Plotting Module. An option to save 
the contour values is offered but not necessary at this 
point. Reselect the CURVE Selection Option and enter the 
desired zeta and omegan curves. Equation (4-2) indicates 
that the zeta-omegan product must be greater than two to 
keep the settling time below two seconds. Therefore, also 
select for plotting a zeta-omegan curve of two. 

Applying the previously discussed procedures’ for 


producing a plot provides the graph of Figure 4-6. 
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Figure 4-6 
s**2 + 500As + 500B 
All alpha/beta pairs falling to the left of the 
zeta-omegan curve satisfy the given settling time design 
requirement. The alpha/beta pair that just meets all 
requirements is the intersection of the three contours. 
This value can more accurately be determined by expanding 
around the position of this point. Using Option 1 of the 
expansion menu and inserting an alpha range from .006 to 


.010 and a beta range from .02 to .06 yields the following 
plot: 
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Figure 4-7 
s**2 + 500As + 500B 


The alpha and beta values determined from the preceding 


graph are: 


a = .008 (4-13) 
B = 038 (4-14) 
Therefore: K = .038 and Kl = .211 


A determination of associated closed loop root 
values for this alpha/beta pair can be accomplished by 


selecting Option 8, ROOT Finder, in MAIN Menu. At the 
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prompt, enter the alpha and beta values. Confirming system 


stability, the outputted system root values are: 


COEF( 3) = 1.00 

COEF( 2) = 4.00 

COEF( 1) = 19.0 
First Root = -2.00 + j} 3.87 
Second Root = =) 2 laa) = 363.87 


2. Example 4-2 


The system response of the previous plant could 
also be shaped with a cascade compensator. Figure 4-4 


illustrates the block diagram of such a system. 





Figure 4-8 
Ideal Instrument Servo with Cascade Compensator 


The same procedures apply for deriving the 
parameter values associated with a particular desired 


system response. 


Let: 
p 
a = = 
| - (4-15) 
Bp = p (4-16) 
The characteristic equation for this system is then: 
s> + Ps? + S00as + S00B = 0 (4-17) 
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Since this is our first look at a third order 
system, a selection of zeta curves ranging from zero to one 
and omegan curves ranging from one to ten was again 


selected. The following graph was the result. 
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Figure 4-9 
Constant Zeta Curves 
s**3 + Bs**2 + 500As + 5SOOB 
As can be seen, most alpha/beta values are negative 
over the selected frequency range from 1 to 10 rad/s. Both 
parameters must be positive to provide all closed _ loop 


roots in the left half of the s-plane. Therefore, an 


alteration in frequency range is indicated. A_ range 


61 


Spanning 10 


to 40 rad/s was entered with the following 





-329,00-80,.00 -UW0,.00 0,00 W0.90 80,00 120,00 180,00 


a) eamie 


Figure 4-10 
Constant Zeta Curves 
s**3 + Bs**2 + 5SOOAS + 500B 


It can now be seen that positive alpha/beta point 


pairs are 


designed. 


possible, indicating that a stable system can be 


However, specification limits cannot be met 


precisely due to the increased natural undamped frequency 


required. Figure 4-11 presents the zeta and omegan contours 


over a frequency range of 20 to 50 rad/s. 
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Figure 4-11 
Constant Zeta & Omegan Curves 
s**3 + Bs**2 + SOOAsS + 500B 
Expanding about an alpha range of zero to ten and 


beta range of zero to 200 yields a more easily interpreted 


graph as seen in Figure 4-12. 
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Figure 4-12 
Constant Zeta & Omegan Curves 
s**3 + Bs**2 + 500As + 500B 
It appears that a minimum omegan value of 25 rad/s 
is required to maintain an acceptable relative damping 
coefficient. Selecting an operating point where zeta equals 
-5 and omegan equals 25 rad/s provides the alpha/beta pair: 


a 6.3 (4-18) 


B = 125 (4-19) 
This corresponds to a compensator with a pole at 125 rad/s 
and a zero at 19.8 rad/s. The resulting closed loop roots 


with this set of parameters is: 
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COEF( 4) = 1.00 


COEF( 3) = 125. 

COEF( 2) = .315E+04 

COEF( 1) = .625E+05 
Root( 3) = -12.65 + 4 21.61 Root( 2) = -12.65 = ee. G1 
Root( 1) = -99.69 + 3} .0000 


Thus, a clear set of dominant roots is present and all 
roots reside in the left half of the s-plane. 


3. Example 4-3 
This example deals with a plant incorporating both 
velocity (tachometer) and acceleration feedback. The block 


diagram of such a system is depicted in Figure 4-13. 





Figure 4-13 
Velocity and Acceleration 
Feedback Compensated System 


The characteristic equation for this system is: 


s> + (3+50a)s? + (2+508)s + 50 = O (4-20) 

where: 
a = Kk, (4-21) 
Se = Cai (4-22) 
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Rather than beginning with a required set of system 
responses, system options will first be examined to 
determine possible zeta and omegan values. Figure 4-14 
provides a look at the combined zeta and omegan graph 


spanning a frequency range from one to 100 rad/s. 
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Figure 4-14 
Constant Zeta & Omegan Curves 
s**3 + (3+50A)s**2 + (2+50B)s + 50 
This system appears to be stable over this 
frequency range. Design specifications requiring a settling 


time under one tenth of a second with an overshoot not to 


exceed ten percent translate into: 
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Games So; een 


Wn = 70 rad/s (4-24) 

Reentering the CURVE Selection Menu and adjusting 

the frequency range to 60-80 rad/s, while specifying curves 

of zeta equals .55, .57 and .60 and omegan equals 65, 70 
and 75 rad/s, produces the plot of Figure 4-15. 
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Figure 4-15 


Constant Zeta & Omegan Curves 
s**3 + (34+50A)s**2 + (2+50B)s + 50 
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Expanding about the desired operating point yields: 
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Figure 4-16 


Constant Zeta & Omegan Curves 
s**3 + (3+50A)s**2 + (2+50B)s + 50 


From Figure 4-16, alpha and beta can be easily 


extracted. Their values are: 


a 1.635 (4-25) 


B = 98.0 (4-26) 
Thus the acceleration feedback coefficient, Kl, is 1.635 
and the corresponding velocity feedback coefficient, K2, is 


98.0. These values provide root locations of: 
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COEF( 4) = 1.00 


COEF( 3) = 84.8 
COEF( 2) = .490E+04 
COEF( 1) = 50.0 
Root( 3) = -42.37 + } 55.73 Root( 2) = 42.37 =) 95.735 
Root( 1) = <-.1020E-01 + j .0000 


Obviously, the dominant roots are not at the 
designed locations. The dominant real root and relatively 
distant complex roots indicate that the system is probably 
overdamped. The omegan value chosen in an attempt to force 
a short settling time is too large. A redesign is indicated 
if a longer settling time can not be accommodated. 

4, Example 4-4 

The final example will examine a forth order system 
to demonstrate the difficulty involved in choosing an 
alpha/beta pair from a set of convoluted contours 
associated with higher order systems. A third order, type 
one plant will be shaped with a cascade compensator. The 


system block diagram is illustrated in Figure 4-17. 





Figure 4-17 
Third Order, Type One System 
With Cascade Compensation 
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The characteristic equation for this system is: 


sf + (11+H)s> + (104+11f)s? + (150a+10f)s + 1508 = O (4-27) 
where: 
o> : (4-28) 
B = p (4-29) 
System specifications call for: 
Peak Overshoot < 50% 
Settling Time < 1 second 


The corresponding relative damping coefficient and undamped 
natural frequency are: 

ae 22 (4-30) 

Wn = 18.2 rad/s (4-31) 

Based on Equation (4-7), the phase margin of a 

system with a zeta value of .22 iS approximately 20 

degrees. Although the system may be stable, it will be very 

sensitive, especially considering that the resonant peak of 

this system occurs at about the same frequency as omegan. 

In any case, with higher order systems it is important to 

perform a time response simulation to more accurately 

determine system behavior. Figure 4-18 is a combined plot 


over a frequency range of 10 to 50 rad/s. 
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Figure 4-18 
Constant Zeta & Omegan Curves 
s**4 + (11+B)s**3 + (10+11B)5**2 
+ (150A+10B)s + 1508 
There appears to be alot of activity about the 


operating point located at zero. Figure 4-19 is a result of 


expanding about that point. 
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Figure 4-19 
Constant Zeta & Omegan Curves 
s**4 + (11+B)s**3 + (10+11B)s**2 
+ (150A+10B)s + 150B 

It is apparent that higher order plots are not as 
well behaved as the lower order systems. To reduce 
confusion and make the graph more readable, fewer contours 
must be selected. Reducing the range of frequencies to a 
span from 10 to 20 rad/s and selecting three zeta and two 


omega curves about the desired operating point produces the 


following graph. 


72 


} 


25). 


<3 8 5L) 


i 
1n.d 1Sob.00 


Cas 


$1] 





OH. 00 


1.00,’ 80,00 180,00 2Y0,00 920,00 400.00 We0,00 680,00 
wa Alpha 


Figure 4-20 
Constant Zeta & Omegan Curves 
s**4 4 (114+B)s**3 + (10+11B)s**2 
+ (150A+10B)s + 150B 
The alpha/beta pair associated with zeta equals .22 
and omegan equals 18.0 follows: 

a = 280 (4-32) 
fe) a 130 (4-33) 


The roots derived through application of these 


values are: 
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COEF (05) 3200 
COEF( 4) = lai. 
COEF( 3) = .144E+04 
COEF( 2) = .433E+05 
COEF( 1) = .195E+05 
Root( 4) = -3.975 + j 17.49 Root ( 3)5=) -—3. 975 


Root( 2) = -.4570 + j .0000 Root( (1) = s-13276 ~ 
Pause - Please enter a blank line (to continue) or a DOS command. 


j 17.49 
j  .0000 


Unfortunately, with that particular selection of 
alpha and beta values, a dominant real root is present. 
Examination of the step input time response indicates an 
overshoot of 50% and a settling time of 1.2 seconds. An 
increase in zeta to better dampen the system results in a 
decreased omegan term if the system settling time is to 
remain at approximately one second. A zeta value increase 
to .3 with an associated omegan value of 14 yields: 

| T= ee (4-34) 


p= 115 (4-35) 


COEF( 5) = 1.00 
COEF( 4 )@=— 126. 
COEF( 3) = .128E+04 
COEF (2) =e. 23 /7E+05 
COEF( 1) = .173E+05 
Root( 4) = -4.218 + j 13.30 Root( 3) = -4.218 13.30 


=a 
Root( 2) = -.7591 + j .0000 Root( 1) = <-116.7 - 3 .0000 
Pause —- Please enter a blank line (to continue) or a DOS command. 


The complex roots have made a relative shift in the 
s-plane toward the dominant real root. Once again, however, 


it can be difficult with high order systems to make a 
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definitive judgment on system behavior without viewing 
transient time response curves. The unit step response of 
this system provides a maximum overshoot of 40% with a one 
second settling time. Thus, a cascade compensator with a 
gain of 150, a zero at 0.76 and a pole at 115 provides a 


suitable transient response. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

As has been demonstrated, rapid solution of a wide 
range of linear automatic control systems problems can be 
achieved through use of the parameter plane program 
developed as part of this thesis. An IBM compatible 
microcomputer program offers the inherent advantages of 
portability and ready access to machines capable of 
Supporting the code. 

The graphical solution of the parameter plane technique 
provides the designer with a visual means of deducing how 
the dominant roots of the characteristic equation move 
about in the s-plane as two user defined parameters 
defining compensator attributes are varied. As a corollary, 
input of specific system response characteristics can be 
made to determine the parameter values providing that 
response. In addition, the parameter plane design procedure 
is independent of system configuration since it deals with 
the characteristic equation only. 

There are some drawbacks to this technique. Only two 
variables can presently be considered in the computer 
model. However, a three variable parameter space method 
could be developed. In place of a series of curves in the 


two dimensional case, the parameter space graphical output 


We 


would be represented by a surface in three dimensions. 
Another shortcoming associated with the parameter plane 
analysis method is that information on open loop poles and 


zeros is not available. 


B. RECOMMENDATIONS FOR FURTHER STUDY 

The most significant addition that could be made to the 
parameter plane program would be the incorporation of a 
graphical cost function analysis package. In a linear 
feedback control system, the difference between the input 
to the system and its output is termed the error. This 
error differs with differing compensators. In most cases 
involving automatic feedback control systems, minimization 
of the error over a certain period of time (e.g., twice the 
settling time) is deSired. Minimization of other system 
attributes is also common. One may wish to minimize the 
absolute value of the error, the square of the error or any 
one of an infinite number of cost functions, depending on 
the particular system application. Constant cost value 
contours can be plotted in much the same presentation as 
the parameter plane plot. Using standardized _ scales, 
desired system response characteristics can be compared 
with the error cost at that operating point and a suitable 
compromise between minimum cost and desired system response 


can be reached. 


a7 


An alternative presentation method is available with 
the plotting package currently used by the parameter plane 
program. The cost value associated with a particular pair 
of alpha/beta values could be plotted as the _ third 
dimension of a parameter space surface. Although precise 
determination of the three variable values would be 
difficult, a quick synopsis of system response trends would 
be available. This would permit ready localization of 
interesting operating areas for more refined analysis ona 


two dimensional plot. 
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APPENDIX 
PARAMETER PLANE PROGRAM 


$NOfloatcalls 
Program ParaPlan 
dimension CONST(10) , ALPHA(10) , BETA(10) ,A(1002) ,B(1002) , ZETA(10) , 


- WN (10) ,AOMEGA(1002) , BOMEGA (1002) , ZWN(10) , AZTAOM( 1002) , 
zs BZTAOM (1002) 

real R,WNMIN,U,U1,U2,B1,B2,C1,C2,D1,D2,CONST, ALPHA, BETA,A,B, ZETA, 
> WN , AOMEGA, BOMEGA , WNMAX , WNMIN1 , ALPHA1 , BETA] , ZWN, AZTAOM, 
* BZTAOM, INCZWN 

integer SEL, SEL3,ORDER,NUMCOF ,NUMZTA,1I,M,N,K,L,J,ROW2,COLM, 
z IOPORT , MODEL, START , STOP, FLAG, NUMWN , STOPO, FLAGRT ,CURV, 
* NUMZWN , STOPZO, OUT 


character*12 FILEIN, FILEOUT 

common A,B,AOMEGA, BOMEGA , AZTAOM, BZTAOM 
common/plottr/ IOPORT ,MODEL 

data FLAGRT/0O/ 


IOPORT = 99 
MODEL = 99 
OUT = 3 
mR = 10. 


10 call INTRMN (SEL) 

Vi W suleewed, 4) .eor. (SEL .é€q. 6) .or. (SEL .eq. 10)) go to 4 
1 call MAINMN (SEL) 
> FLAG = 0 


3 if (SEL .eq. 1) then 
call CRVSEL (CONST,ALPHA, BETA, ZETA, NUMZTA, STOP, NUMCOF, 
= WNMIN , WNMAX , FLAG, SEL, NUMWN, WN, 
- STOPO , ORDER , ZWN , NUMZWN , STOPZO) 


else if (SEL .eq. 2) then 

call MONPRT 
else if (SEL .eq. 3) then 

call REVCHG (ORDER, NUMCOF , CONST, ALPHA, BETA, WNMIN,WNMAX, SEL, FLAG) 
else if (SEL .eq. 4) then 

call LDFILE (ORDER,NUMCOF, CONST, ALPHA, BETA ,WNMIN , WNMAX) 
else if (SEL .eq. 5) then 

call SAVFIL (ORDER , NUMCOF ,CONST, ALPHA, BETA,WNMIN , WNMAX, FILEOUT) 
else if (SEL .eq. 6) then 

call CHAREQ (ORDER,CONST,ALPHA, BETA, NUMCOF ,WNMIN,WNMAX,FLAGRT, 
= ALPHAI1 , BETA] ) 
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Cae OO) QO 


QO 


else if (SEL .eq. 7) then 
call RUNPRO (NUMZTA, STOP, ZETA, CURV,NUMWN, WN, 


> STOPO, ZWN , NUMZWN , STOPZO) 
else if (SEL .eq. 8) then 
FLAGRT = 1 
call ROOTS (ORDER,CONST,ALPHA, BETA ,NUMCOF,WNMIN, WNMAX, 
* ALPHA1 , BETA1 ,FLAGRT, OUT, R) 
PAUSE 
FLAGRT = O 
else if (SEL .eq. 10) then 
call SAMPLE 
goto 10 
else 
go to 999 
2 end if 
if (FLAG .eq. 1) go to 4 
go to l 
999 stop 
end 
end 


9K OK Kk OK ok KK ok ok KK ok KK ok OK OK OK ok OK ok OK OK ok 2K OK OK OK OK OK 2K ok OK ok OK ok ok OK ok OK OK KK OK ok OK OK ok OK OK KK OK ok OK KOK OK KK OK ok OK OK OK OK OK OK OK 


USER UTILITIES MODULE 
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Subroutine INTRMN -- allows input of C. E. from file or keyboard 


Sm mm we i ee we 
ee SS SS SS SS SS i ae iS ee a ee a a SS SS SS 


Subroutine INTRMN (SEL) 
integer SEL 
78 call CLR 
write (*, 80) 
write(*, 81) 
write(*, 82) 
write(*, 81) 
write(*,83) 
write(*, 81) 
write(*,93) 
80 format(1x,///////,10x, ° -------------------------- ----------------- 





81 format (10x, ° --------------- - - - -- - 


82... format(10x,” |° , 2x, “OPTION.NO.~ , 2X eae ce Oar T LON) xen) 
83 format (lOsx. |e 2 ih ‘ 2x,°!1’,5x,°’LOAD Problem from File’, 








80 


KIO Re 

















co LOR, ~ | PR 2 *,2x,°]° ,5x,°INPUT Characteristic Equatio 

tm, 1s, | * 

* SMO KOA , Dxept 3 ‘,2x,°’|°,5x,°EXAMPLE Characteristic Equat 

*vemminput’ ,5x,° |” 

EwyMeOx ,- | * 5. 2x, * 9 7, 2x,°1°,5x,°EXIT to Main Menu’ ,25x,’|7) 
93 format(1x,///,15x,’Enter integer number for selection ===> ”) 


call PSIT(22,56) 
read(*,*,err=78) SEL 
if (SEL .eq. 1) then 


SEL = 4 

else if (SEL .eq. 2) then 
Sebo 

else if (SEL .eq. 3) then 
SEL = 10 

else 
continue 

end if 

return 

end 

end 

C Subroutine MAINMN -- allows selection of items from the main menu 


am rn ww em i ee 
SS SS eS eS i ee ee a a a SS SS i. a SS iS EE SE SS ES SS ES SS 


Subroutine MAINMN (SEL) 
integer SEL 
79 call CER 
write(*, 80) 
write(*,81) 
write(*, 82) 
write(*, 81) 
write(*, 83) 
write(*, 81) 
write(*, 84) 
Wwiete (* 484) 
write(*, 85) 
write(*, 87) 
write(*, 86) 
write(*, 81) 
write(*,93) 
80 PoOniawil 7 (Ox. = _ _ Se _____ ae 
Danaea a ye Ox, |) 25x, MALIN MENU (28x, |~*) 
81 HOI LOR Oe a ee 


8 | 



































87 format (LOX, ~ FA+FEFEFFHFEFE FEAF HEHEHE FF FAF4A FHF 4444444444 444444444444 
*+++44444+ ) 

82 format(lOx, | 7233 OF Tlon NO.” , 2x2 |" , 20x, “OPTEON@ = 2 ia 

83 format( 10x, | 42 1 *,2x,’|°’,5x,°CURVE Selection Menu’, 
*2 One ale 
a / ROX sa exe 2 ° 2x,’|’,5x,’ PRINTER Selection Menu’ 

A PAD) ae, | 
e/a Omen ee cx oe 3 7 ,2x,°|%,5x, “REVIEW/CHANGE Selections’ 
aS ecl |) 

84 format (1Oxaelh 2x) 4 ‘,2x,°|°,5x,°LOAD Problem from File’, 
* DOE | 4 
Oe eene Xe 5 7 2x,°’|%,5x,°SAVE This Problem’ , 
2S ee) 

AGO X57 ee ae 6 *,2x,°]°,5x,° INPUT Characteristic Equa 
trices es 17) 

85 forma til Ox. «|| 2 xue 7 *,2x,°|°,5x,° PLOT Curves (aue.ae 
¥ of VOXG | eee 8 *,2x,°|*,5x, “ROOT Findemag sae) 

86 forma t (Ox eay 2 xe 9 *,2x,°|*° ,5x, EXIT Program ,30% ae 

93 format(1x,///,15x,’Enter integer number for selection ===> ”) 

call PSIT(24, 56) 

read(*,*,err=79) SEL 

return 

end 
C 
C sseesesenn2s22s2225282222222282558888828552525525 525525888 8888555===== 
C Subroutine LDFILE -- allows loading of data from an existing file 
C ssssesas2s2222s222=522222552885288588555 55552585 S5e5e88 228888585 255===2== 
C 

Subroutine LDFILE (ORDER,NUMCOF ,CONST,ALPHA, BETA, WNMIN,WNMAX) 

real WNMIN,WNMAX,CONST(10) , ALPHA(10) , BETA(10) 

integer ORDER,NUMCOF ,MODEL, IOPORT, TEST, I 

character*12 FILEIN, ANS*¥1 

logical EXIST 

common/plottr/ IOPORT,MODEL 

1 call CLR 

write(*,100) 

100 format(1x,////,10x, * ---------------—-------------- --- -- --- --------- 
*—————— “,/,10x,’*** This is the routine to LOAD data from a file ** 
KK. 7 10x, * ----------------------- ----- - - - - - - + “ 
*7/,t0x, File name should net exceed 6 chanacvensa, « 

*10x,°File extension should not exceed 3 characters’ ) 
write(*,101) 
101 format(1x,////,10x,’What is the file name (fn) and extension (ext) 


*? °  //,10x,’Enter in form fn.ext (e.g. PARAPLAN.INP) ===> ” ) 
call PSIT(18,57) 
read(*,102) FILEIN 
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102 format(al2) 
inquire(file=FILEIN,exist=EXIST) 
if (EXIST) then 
open(7,F ILE=FILEIN , STATUS=’ OLD’ , ACCESS=’ SEQUENTIAL’ ) 


else 
write(*,103) 
103 format(///,10x,’That file does not exist.’ ,/,10x, 
* “Do you want to try again? (Y/N)’) 


read(*,°(Al)*) ANS 
tee CANS .@€q. "1 )  -Or. (ANo .eq.  y )) then 
goto 1 
else 
goto 2 
end if 
end if 
read(7,105) MODEL, IOPORT 
read(7,110) ORDER,NUMCOF 
read(7,115) WNMIN,WNMAX 
do 190 I = 1,NUMCOF 
read(7,120) CONST(I) ,ALPHA(I) , BETA(I) 
190 continue 
euese (7, STATUS=" KEEP”) 


moo ftormat(1x,12,2xX,12) 
moe tormat(1x,1I2,2X,12) 
imeoee tTormat(1x,2F15.7) 
ieee format(1x,3F15.7) 
write(*,3) WNMIN,WNMAX, ORDER 


3 fommatitx,///////7,0x%, WNMIN = ° ,f12.3,5x,°WNMAX = * ,f12.3,5x, 
* SORDERS= ~.,.12%/7) 
do 9 I = NUMCOF,1,-1 
J=I-1 
write(*,4) J,CONST(I),J,ALPHA(I),J,BETA(I) 
9 continue 
4 Peano CONST( ,12, je=" | flO.3,5x, ALPHA( ,12,° )™= °,f10.3,5x, 
* eBerAtC ,i2. ) =" ,f10.3) 


write(*,5) 
5 hermac(ix,//////7/) 


PAUSE 
C 
Z return 
end 
C 
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Subroutine SAVFIL -- allows entered data to be saved to a data file 


100 


FO 


102 


103 


i i SS Se ee 
eSelSCCececCcCcClfCQcQSa SS SS SS SS SS Le ee ee ee ee ee ee ee ee ee ee ee ee en ee ee ee ee ee ee i SS SS SS ES LS SS ae ae ee Ee ie cee ee ee ee ee eee eee ee ee 


Subroutine SAVFIL (ORDER, NUMCOF ,CONST, ALPHA, BETA,WNMIN, WNMAX, 
FILEOUT) 

real WNMIN,WNMAX,CONST(10) , ALPHA(10) , BETA(10) 

integer ORDER ,NUMCOF ,MODEL,IOPORT, TEST,1I,J,SEL,SEL3 
character*12 FILEOUT,ANS*1 

logical EXIST 

common/plottr/ IOPORT,MODEL 


Calle CeER 

write(*,100) 

format(1x,////,10%,° -—-----=--=---=—— === eee 
*x*---—’ /,10x,°*** This is the routine to SAVE data to a file ***’ 
*,/,10K, =S- >= => == ---- en ne ee . 


* //,10x,’File name should not exceed 8 characters’ ,/, 
*10x,’File extension should not exceed 3 characters’ ) 
write(*,101) 
format(1x,////,10x,’What is the file name (fn) and extension (ext) 
*? “,//,10x,’Enter in form fn.ext (e.g. PARAPLAN.OUT) ===> 7) 
Call SeSITCs 357) 
read(*,102) FILEOQUT 
format(al2) 
eall CLR 
inquire(file=FILEOUT, exist=EXIST) 
if (EXIST) then 
write(*,103) 
format(///,10x,’That file already exists.’ ,/,10x, 
* “Do you want to overwrite? (Y/N)”) 
read(*,*(Al)*’) ANS 
1f((ANS .eq. “Y’) .or. (ANS .eq. ’y’)) then 
open(7, file=FILEOUT, status=’ OLD’ ) 
rewind 7 
goto 1 
else 
goto 2 
end if 
end if 
open(7,file=FILEOUT, status=’ NEW’ ) 
write(7,105) MODEL, IOPORT 
write(7,110) ORDER,NUMCOF 
write(7,115) WNMIN,WNMAX 
do 190 I = 1,NUMCOF 
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write(7,120) CONST(1I),ALPHA(I) ,BETA(TI) 
190 continue 
close(7,STATUS=’ KEEP’ ) 


mos) format (1x,12,2xX,12) 
imo format (|1x,12,2xX,12) 
ies fommat(1x, 2F15.%9 
mzO format(1x,3F15.7) 
write(*,3) FILEOUT 


3 format(1x,////////,10x,’File name and extension = ” ,al2,///) 
write(*,4) WNMIN,WNMAX, ORDER 
4 format(5x,°WNMIN = “” ,g11.4,5x,”“WNMAX = “ ,g11.4,5x, 
* SOR DEHe= “| 12777) 
do 9 I = NUMCOF,1,-1 
J=tI-l1 
wet eG 5 oot CONST (1) seh, ALPHA(1) , JgBETA(I1) 
9 continue 
5 memiaci@rn, COMoia( 12°) = ~ £10.3,.5x, “ALPHA(” ,12,°) = ° ,#10.3,5x, 
* OBmbAt ele.) = ~f10.3) 


write(*,6) 
6 BOrMakaulx. // 77 //) 


PAUSE 

return 

end 
C 
Ee ee eee eee 
C Subroutine CHAREQ -—- allows input of characteristic equation 
Neen eee eee ene eee ee ea eH a oe eS eS = = SS = SS == === 
C 


Subroutine CHAREQ (ORDER,CONST, ALPHA, BETA, NUMCOF ,WNMIN,WNMAX, 
* FLAGRT , ALPHA1 , BETA1 ) 
real CONST(10) ,ALPHA(10) , BETA(10) ,WNMIN,WNMAX , ALPHA] , BETA1 
integer ORDER,I,NUMCOF,IMIN1,ROW1 , FLAGRT 
character*1 CHG 
200 call CLR 


if (FLAGRT .eq. 1) then 
write(*,190) 
190 format(1x,//,9x,°Do you want to find the roots of the existing Cha 
*xracteristic’ ,/,9x, “Equation, entering values for ALPHA & BETA?’ , 
a eos Enter “y" or “"n" ==> *) 
call PSIT(7,30) 
read(*,° (A)* ,err=200) CHG 


end if 
iwiehiGeseq. ~“y ) .Or. (CHG .eq. ~Y )) then 
go to 293 
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201 


210 


PALA 


220 


PAALE| 
els 


230 


P| 


235 


239 


240 


end if 
call CLR 


*** Enter ORDER of Characteristic Equation *** 
write(*,201) 
format(1x,////,10x, “What is the order of the Characteristic Equati 


*on? *,//,10x,’Enter integer value less than or equal to 9 ===> ”) 


call PSIT(8,59) 
read(*,*,err=200) ORDER 
NUMCOF = ORDER + 1 


*** Enter CONSTANT Coefficients *** 

call CLR 

ROW1 = 10 

write(*,211) 

format(1x,////,5x,’Enter the CONSTANT Coefficient Values of the Ch 


do 217 I = NUMCOF,1,-1 
IMIN] =1I-1 
write(*,220) IMIN1 
format(1x,/,5x.° CONSTANT Coefficient of S **° ,12, == 
call PSIT(ROW1,40) 
read(*.*err—210) “CONST (1) 
ROW1 = ROW] + 2° 
continue 
call CLR 
write(*,230) 
format(1x,////,20x, ***** CONSTANT COEFFICIENTS, ***22 78 
do 235 I = NUMCOF,1,-1 
IMIN] = [1-1 
write(*,237) IMIN1,CONST(TI) 
format (26x, “CONST @ = !2- j=——. 197377) 
continue 
ROW1l = ROW1 + 3 
write(*,239) 
format(1x,//,5x,°Do you want to make any CHANGES to the CONSTANT C 


*oefficient Values? .// 5x. Enter "y)'s0re'n  ——-a. 


call PSIT(ROW1, 26) 

read(*,°(A)’ ,err=218) CHG 

if ((CHG .eq. “y’) .or. (CHG .eq. ’Y’)) then 
Po to 210 

end if 


*** Enter ALPHA Coefficients *** 
call CLR 
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ROW] = 10 
write(*,241) 
241 =format(1x,////,5x,’Enter the ALPHA Coefficient Values of the Chara 


do 247 I = NUMCOF,1,-1 
IMIN1] = I-11 
write(*,250) IMIN1 
250 format (lx,/,5x. ALPHA Coefficient of S **’ ,12,° =) 
call PSIT(ROW1, 39) 
read(*,*,err=240) ALPHA(I) 
ROW1 = ROW1 + 2 
247 continue 


call CLR 
write(*,260) 
Peomeronmat¢ix,////,20x,° ***** ALPHA COEFFICIENTS *** 8" ,//) 
do 265 I = NUMCOF,1,-1 
IMINI = I-11 
write(*,267) IMIN1,ALPHA(TI) 
267 format (26x, ALPHA(” ,12,°) = °,£9.3,/) 
265 continue 
ROWL = ROW1 + 3 
write(*,269) 
269 format(1x,//,5x,°Do you want to make any CHANGES to the ALPHA Coef 
Prrelrent Values?  ,//,5x, Enter "“y" or "n" ==> ”) 
call PSIT(ROW1, 26) 
read(*,° (A)” ,err=240) CHG 
mi ¢(CHG 2eq. ~y ) .or. (CHG .eq. °Y’)) then 
go to 240 
end if 


*** Enter BETA Coefficients *** 
e270), call CLR 
ROW1 = 10 
write(*,271) 
271 format(1x,////.5x,’°Enter the BETA Coefficient Values of the Charac 
*teristic Equation’ ,/,5x, °“---~------- ==a=——-——-—-—--—--—------—----------- 


do 277 I = NUMCOF,1,-1 
IMIN1 = I-11 
write(*,280) IMIN1 
280 hormmat¢(lx,/.ox, BETA Coefficient’ ors ** Kr,” = ”) 
call PSIT(ROW1, 38) 
read(*,*,err=270) BETA(I) 
ROW1 = ROW1 + 2 


87 


PEEL 


290 


294 
292 


299 


293 


370 
371 


372 
373 


QO 


2907 


298 


296 
299 


continue 
call CLR 
write(*,290) 


format (1x,//// (20K, ***** BETA COBPFICIENRSS 22". 

do 292 I = NUMCOF,1,-1 
IMINI = I-11 
write(*,294) IMIN1,BETA(I) 
format(z26— BETAG 42. je=" - {9.3.79 

continue 

ROW1 = ROW1 + 3 

write(*,295) 

format(1x,//,5x,’Do you want to make any CHANGES to the BETA Coeff 

*¥icient, Values? “.//.ox%- eater “vy 'sor "ns )——-a®) 

call PSIT(ROW1, 26) 

read(*,’° (A)”’ ,err=270) CHG 

if ((CHG .eq. “y’) .or. (CHG .eq. ’Y’)) then 
go to 270 

end if 

call CLR 

CHG = ‘n’ 

if (FLAGRT .eq. 1) then 
write(*,371) 
format(/////,10x,’Enter ALPHA value ==> ’,/,10x,’All ALPHA coef 


*ficients multiplied by ALPHA’ ) 


call PSIT(7,34). 

read(*,*,err=370) ALPHA1 

write(*,373) 

format(////,10x,’Enter BETA value ==> ’ ,/,10x,’All BETA coeffiic 


*ients multiplied by BETA’) 


call PSIT(12,. 33) 
read(*,*,err=372) BETAL 
end if 


**¥* Enter RANGE of frequencies to examine *** 
if (FLAGRT .ne. 1) then 


write(*,297) 
format(1x,//,5x,’Enter the RANGE of frequencies you wish to exam 


*ine’ f/f, 5X, -H a on) 


write(*,298) 

format(1x,//,5x,’Enter minimum frequency (WnMIN) ===> “ ) 
cally Psir.10, 43) 

read(*,*,err=293) WNMIN 

write(*,299) 

format(1x,//,5x,’Enter maximum frequency (WnMAX) ===> ” ) 
call PSIT(13,43) 
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81 
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read(*,*,err=296) WNMAX 
end if 
return 
end 


em wm wc ce cw me me ee ee we ws es ee wes ei we a ee i 
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Subroutine MONPRT -- allows display on monitor or printer 


mm ee em ee eee ee i i i i ee OS 
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Subroutine MONPRT 
integer IOPORT,MODEL, SEL 
common/plottr/ IOPORT,MODEL 


call CLR 

write(*, 80) 

write(*, 81) 

write(*, 82) 

write(*,81) 

write(*, 83) 

write(*, 87) 

write(*, 84) 

write(*,81) 

write(*,93) 

format (1x,///,10x, ° ------------------------------------- ---------- 
* —___—_————-—-—---—---— Ty lOx Ir 21x eERINIER; OUTPUT MENU 22x,°|*) 








format (LOX, ~ FEFFEEATHEEAFEEEFFFEAFEEEEEHEEFF EEE FEE HEF EEE THF + +++ +4 +44 
K++++++4+44+7 ) 























monmmac@lox, | ,2x,° PRINTER NO As 2x lee 20x] CPR ENTER alo x, “) 
mormat¢(l0x,° 1% ,3x,° 1 - , 2x. 4x, (Epson PxX—80-8 Al)” 
eeOOx, |” 

ceeeel@x. «| , 3X,” 2 oW2x,  |° 4x, Epson FX=f00, “er 
fox, |” 

ee elOX, 1° 3X, ~ 3 ‘,2x,°|°,4x,’Epson MX-100, All’ 

eee ox,; |” 

aw 10x |, 3x,” 4 7 2x,’|’,4x,’Epson RX-80, All’ 
Pe2GR | 7 

eanOX, “i . 3X, ~ 5 ‘,2x,’|%,4x, “Epson MX-80 & IBM Printer’ 
welax, |” 

eeelOx, 1° , 3X, ° 6 “2%, | , 4%, "HP 7470f Graphics Plotter’ 
“7a Ey ea Ve 

cow Ora 1 3x, 7 “,2x,°|*,4x,°HP 7475A Graphics Plotter’ 
_ na ale ara 

ee emt neail 3X, * 8 2x0) |e HP a oe xBeoceries pitet ters” 
eS x |< 
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sears ora (Op, Gemma Mller Oe 9 “,2x, |°,4X%, HP 2686A Laser Ver 
* 24x, 7p) 
84 format (Os ange 
ne 10 “,2x,°|°,4x,°’Graphics Monitor (default) ’ 
oat IL (6 re alee 
ef OM os iba ’ 2x,°|1°,4x,° HARDWARE Interface Menu’ 
* TOK he 
«Sf 10x |, ok ‘.2x,’|’,4x,’Input PLOT88 Values for IOPO 
*RT and MODEL’ 
Se xe a 
eo LOX Wea oe 99 *,2x,°|° ,4x, EXIT to™Main Menu’ ,25xo ae 








93 format(1x,///,15x,’Enter integer number for selection ===> “) 
call PSIT(24,56) 
read(*,*) SEL 
read(*,*) SEL 
if (SEL .eq. 1) then 


IOPORT = O 
MODEL = 5 

else if (SEL .eq. 2) then 
ITOPORT = O 
MODEL = 15 

else if (SEL .eq. 3) then 
IOPORT = O 
MODEL = ll 

else if (SEL .eq. 4) then 
IOPORT = O 
MODEL = 1 

else if (SEL .eq. 5) then 
IOPORT = 0 
MODEL = 1 


else if (SEL .eq. 6) then 
IOPORT = 9600 
MODEL = 20 

else if (SEL .eq. 7) then 
ITOPORT = 9600 
MODEL = 30 

else if (SEL .eq. 8) then 
ITOPORT = 9600 
MODEL = 80 

else if (SEL .eq. 9) then 
IOPORT = 9600 


MODEL = 60 

else if (SEL .eq. 11) then 
call PORT 

else if (SEL .eq. 12) then 
call CLR 
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write(*,100) 


100 format(1x,////,5x,°’ Input value to be used by the PLOT88 graphics 
* package’ ,/,5x,°’for IOPORT ==>’ ,//) 
call PSIl(7,21) 
Freaa(*,*) IOPORT 
write(*,110) 
10 format(1x,///,5x,’Input value to be used by the PLOT88 graphics 
package” ,/,5X,° for MODEL ==>‘u//) 
gad) PSIT(12,20) 
read(*,*) MODEL 


else 
IOPORT = 99 
MODEL = 99 
endif 


write(*,98) IOPORT,MODEL 
98 format(//,10x,’PLOT&8 will use these values to output graphics:’ 
* //,10x,”° IOPORT = ’” ,I14,10x,°MODEL = °” ,I4,///) 
PAUSE 
return 
end 


mmm Sm cs cw ww ww we ee we ww we we ww ww ai ia ie ase i i 
———— eo OH ee ee ee SS SS Se SS SS SS SS ES SS SS ES SS SS LS SS SS ES Se Se 


Subroutine PORT -- allows selection of output printer port 


ee ae eae a a i i iS SES 
ee rr 


Subroutine PORT 
integer IOPORT,MODEL, SEL 
common/plottr/ IOPORT,MODEL 


call CLR 
write(*,79) 

79 format (1x,///,15x,’Selection of a printer in the previous menu aut 
*omatically’ ./,10x,’selects the most commonly associated output por 
Bueoe it parallel ’ ,/,10x, 

*“ (LPT) or serial, for that particular device. If the program graph 
mags ,/, 10x, 

*’are not being output correctly, use this menu to properly route t 
*ne ae , 10x, 

*“sraphics data to the output device.’ ,////) 

PAUSE 


call CLR 

write(*, 80) 
write(*, 81) 
write(*, 82) 
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80 


81 


87 


82 
83 


84 
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Write (+ ol) 
write (*, 83) 
write(*, 87) 
write(*, 84) 
write(*, 81) 
write(*,93) 
format(1x,///,10x, ° ----------------------- --- - - - ----- ------ ---- 


Sr AOr Seal 
format (10x, * —---—————-——_-- + == — er 


Format (LOX, ~+FFEFEEFFEFEFEFFEFEFFEFEFEFFFEEFFEFEFEFE EEE E FE 4444444+4444 
F++4++++4++447 ) 


format(10x, ’ fait) 








7, 2x, SELECT NO.” 3x, 1°, 20x) PORE =) 27,0 























format (10x. 4. 3x. 1 * 2x,’|’ ,14x,°LPT1 Printererone 
x LO xe: 

telox,. |, 3x, 2 “2x9 |* 14x, “CPT2" Printer Porm 

+, Lox 

kJ NOR ine BX 3 “,ex, |°, 14x, °“LPT3 Printer Pore 
wreihar.a (ie 

ef Oa, | OX, 4 *,2x,°|* ,14x,°COM1 Serial Pore 

* 16x,7 1" 

tO ano Oe 5 7, 2x, °|° ,14x,°COM2 Serial Porte 

* Lexan) 

forma lO xe ans oe 

AE 99 72%, |, Waxy EXIT to Main) Menu <5.) 
format(1x,///,15x,’Enter integer number for selection ===> ” ) 


call PSIT(21,56) 
read(*,*) SEL 
if (SEL .eq. 1) then 


IOPORT = 0 

else if (SEL .eq. 2) then 
IOPORT = 2 

else if (SEL .eq. 3) then 
IOPORT = 3 

else if (SEL .eq. 4) then 
IOPORT = O 

else if (SEL .eq. 5) then 
IOPORT = 50 

else 
IOPORT = 99 

endif 


if ((SEL .eq. 4) .0or. (SEL .éq. 5)))_then 
if ((SEL .eq. 4) .or. (SEL .eq. 5)) then 


call CLR 
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60 


61 


62 
63 


50 


51 


write (*,60) 
write(*,61) 
write(*,62) 
write(*,61) 
write(*,63) 
write(*,61) 
write(*,93) 


format (1x,///,10X, ° ------~----- 


plex,“ |“) 


homnrenn@lox, | °,2x,° SELECT NO.* ,3x,°|°,1!5x, BAUD RATE”, 22x, “ 


format ( 10x,” |*,3x," il 
oa, | 
mer Ox |’ , 3x,” 2 
sO ax ||| 
eelOox, | 43x," 3 
eeea x, 
mee, *|~ ,3x,° 4 
mmoax, “| ) 
call PSIT(18, 56) 
read(*,*) SEL 
if (SEL .eq. 1) then 
IOPORT = IOPORT + 300 
else if (SEL .eq. 2) then 
ITOPORT = IOPORT + 1200 
else if (SEL .eq. 3) then 
IOPORT = IOPORT + 4800 
else 
IOPORT = IOPORT + 9600 
endif 








call CLR 

write(*,50) 
write(*,51) 
write(*,52) 
write(*,51) 
write(*,53) 
write(*,51) 
weerre( * O03) 


ES SEN (MOS fo LS It 


2) 





Ox, 1, 18x O00" 


28 oo ese, 7 A Our 


SOK | xX 4 8O00- 


any aad 





“ , W8x , AS60000" 


—_— ne *./,10x,’|%,25x, “PARITY MENU’ 


Prax |) 


ae 


o2 


53 


99 


10 


20 


format (10x, ” 





*, 2X, ° SELECT NO.” , 3x, |° , 17x, Pande ee 





“) 








format( lox, |. cx iL *,2x,°|* , 14x," NO@ean it, 
ale 2) 2p ea 

#7 OX pens Xo 2 2x,’ |* , 14x, “EVEN Panaty © 
eles nine 

* oO Xana eo 3 “2%, |° ,14x,° ODD Pawit 

+ 2 ae) 


eal Pst 7.50) 

read(*,*) SEL 

end if 

if (SEL eq. 3). then 
IOPORT = IOPORT + 1 

else if (SEL .eq. 2) then 
IOPORT = IOPORT + 2 

endif 

return 

end 
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SS SS ESS OS OS SS ee ee ee SE OS SS ES ES ES ES ES SS ES SS SS ES ES SS Se ee Oe Se ee ae Ger ee ee ee ee ee ee SS OS Se Se ee ee ae 


Subroutine SAMPLE -- demonstrates input of characteristic equation 


ems me ee a i a a eee i i 
ee 


Subroutine SAMPLE 


Ga lie cer 
write(*,10) 
format(1x,////,1x, ° --<----------- -- -- - -- - - - -- - - - - -- -- -- -- -------- 


*°*** This subroutine demonstrates loading a Characteristic Equatio 
¥*n is fa: 


write(*,20) 

format (1x, 

*“A universal third order system is represented by the characterist 
*ic equation: ,//,20x, sS**3 + As**2 4+ Bs + 1 °="3@0@09 771. 

*“A and B are variables representing user defined compensator param 
*etiers. ~,/7/ , 10x, 

*“The characteristic equation is obtained with a 1/s**3 plant’ ,/, 
*5x, “incorporating both acceleration (s**2) and velocity (Ss)’,/ 
*,5x,° feedback and assigning variables A and B as the respective’ ,/ 
* 5x,° gain terms of the two compensators.’ ,//,10x, 

*“To load this characteristic equation into the computer model firs 
a Oe 

*’select option 2 in the introductory LOAD/INSERT Menu or option 6 
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30 


40 


iniecne , /.oX., 


*°“MAIN Menu. Both options are titled "INPUT Characteristic Equation 
ee //) 

PAUSE 

Gall CLR 

write(*,30) 

fanmnmauclx,//, 10x, 

*’“You will be asked to insert CONSTANT, ALPHA and BETA coefficients 
pee / OX, 

*“In the case of the universal third order system, the appropriate 
xresponse’ ,/,5x,°following each prompt is:’ ,//,20x, 


*“CONSTANT Coefficient of S ** 3 = lige / W20x, 
*“CONSTANT Coefficient of S ** 2 = OArs/ 20 x 
*“CONSTANT Coefficient of S ** 1 = O-, 20x, 
*“CONSTANT Coefficient of S ** 0 = ome, 2ZOX., 
*“ ALPHA Coefficient of S ** 3 = QO”, /M20x,, 
*“ ALPHA Coefficient of S ** 2 = te , CO 
*“ ALPHA Coefficient of S ** 1 = Oy 20x. 
*° ALPHA Coefficient of S ** 0 = O° .//, 20x, 
*° BETA Coefficient of S ** 3 = Oo / , 20x. 
““peTAPCoefficientwof S$ ** 2 = 0” .7,20" 
*°“BETA Coefficient of S ** 1 = law, 20x, 
> BERAMCOGERicient of S ** 0 = O* ,/ , 20x eez/ /) 
PAUSE 

call CLR 


write(*,40) 

format (1x ,//,5x, 

*“When each set of responses iS complete, they will be echoed back 
ear” ,/,5xX, 

es verification.~ ,//,10x, 

*“You will next be asked to enter the minimum and maximum frequenci 
Fes’ ,/,5x, 

*’you wish to examine. This is the natural undamped frequency range 
BeeOver ./,5X, 

*’“which the relative damping coefficient (zeta) will be calculated 
ace , /, 5%, 

*’construct the constant zeta curves.’ ,//,10x, 

*’Generally, the minimum and maximum frequencies will be dictated b 
 veehne ./,5X., 

*“desired system response. For example, we wish to have a system se 
Setlang ,/,5x, 

*“time of less than 2 seconds. Settling time is approximated by the 
* equation’ ,/,5x, 

*“4 divided by the product of relative damping coefficient and unda 
*mped’ ,/, 5x, 


95 


*’“natural frequency. If we wish to look at a range of zeta values f 
* TOM a ey Ose 

*“to 1, we can solve for the applicable min and max frequencies. In 
Pet Saf oe 

*’case they are 2 Hz and 20 Hz.’,////) 


PAUSE 
return 
end 
C 
C sss=s=<<<=2=2=========$2===S5====555=5555S5555555555 5825555555555 2255522=== 
C Subroutine REVCHG -- allows review or change of entered data 
C seseeesee2sse222e222222255S22e5255e5252 2282822285255 eeSSee5See22555== 
C 
Subroutine REVCHG (ORDER, NUMCOF , CONST, ALPHA, BETA,WNMIN,WNMAX, 
SEL, FLAG) 
real CONST(10) ,ALPHA(10) , BETA(10) ,WNMIN, WNMAX 
integer CHG,ORDER , NUMCOF, FLAG, SEL, ROW 
call CLR 
90 write(*,100) 
100 format(1x,///,10x,’***** Characteristic Equation Coefficients **** 
ee) 
write(*,3) WNMIN,WNMAX, ORDER 
3 format (5x, “WNMIN = ~,f12.3,5x%, ° WNMAX = “2 fl2237 5-5 
“ORDER = ’ ,1I2,//) 
do 9 I = NUMCOF,1,-1 
J=I-1 
write(*,4) J,€CONST(1) .J7ALPHA CG) BETA) 
2) continue 
+ format(1X, “CONST(’ ,12,°) = ~ ,£10.3,5x, “ ALPHA(’ ,12,7) =  1tGeeaae 
“BETA(*22—" )) =" fis) 
write(*,5) 
5 format(1x,//,5x,°Do you want to make any CHANGES to the Coefficien 


*t Values? .//7 ,5x, Enter" or "n'3==)343) 


ROW = 16 + NUMCOF 

call PSIT(ROW, 27) 

read(*,°(A)*,err=90) CHG 

if ((CHG .eq. ’y’) .or. (CHG .eq. “Y’)) then 
Spee — 66 
FLAG = 1 

end if 

return 

end 
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CURVE SELECTION MODULE 
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Subroutine CRVSEL -- allows selection of desired constant curves 


cm rr mw 
ee mm re ee ee ee ee SS ee 


Subroutine CRVSEL (CONST, ALPHA, BETA, ZETA,NUMZTA, STOP, NUMCOF , 


* WNMIN , WNMAX, FLAG, SEL , NUMWN , WN, 
STOPO , ORDER, ZWN , NUMZWN , STOPZO) 
real R,WNMIN,U,U1,U2,B1,B2,C1,C2,D1,D2, ZETA(10) , INCZOM, 
> WNMAX , WNMIN1 , CONST (10) , ALPHA(10) , BETA(10) ,A(1002) ,B(1002), 
: INCLOG , MAXLOG ,MINLOG,WN(10) , INC1, ZETA1 , AODMEGA( 1002) , INCZWN, 
BOMEGA (1002) , ALPHA1, BETA1 , ZWN(10) , AZTAOM(1002) , BZTAOM( 1002) 
integer SEL,SEL3,ORDER , NUMCOF ,NUMZTA,I,M,N,K,L,ROW2,COLM, 
ee START , STOP, FLAG, FRQCHG , NUMWN, TABLE , 
* STARTO, STOPO, ZTASEL , FLAGRT , NUMZWN , STOPZO , OUT 


character*l ANS 

character*12 FILENAME 

common A,B,AOMEGA, BOMEGA, AZTAOM, BZTAOM 
common/filenam/ FILENAME 

logical EXIST 

FILENAME = ” . 


call CLR 

write(*,21) 
write(*,32) 
write(*,22) 
write(*,32) 
write(*,23) 
write(*,32) 
write(*,38) 

Poumawelx.,/, 10x, =—--—————————__—-_______ 
*————---------—_— *,./,10x,’|’,16x,”’CURVE DATA POINT DISPLAY MENU’ , 
myx, |") 

mMonmatclox, |||° ,2x,° OPTION NO.” ,2x,°|°,15x,° OUTPUT SELECTION’ 
mealox, ||” ) 


mormauc¢lOx, | ,2x,~ 1 “22x. |", 5x, NO GUL pUL™RDESPLAY or sav 
*e to FILE’ ,9x,’ |’ 

men Ox, |° , 2x; ~ 2 ‘,2x,°’|’,5x,’°DISPLAY every 10th alpha/bet 
~emvaliue” ,7x, ~ |’ 

x io, 1° 2x,” 3 ‘,2x,°|°,5x,°DISPLAY alpha/beta values an 


Oe mGots 7x, | 
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28 


29 
101 


102 


103 


196 


30 


€ 4/10) See 4 *,2x,’|’,5x, “DISPLAY & FILE alpha/berasun 
*nd roots” , 6a, | 

te 1 OX) ae os 5 ‘,2x,°[|[’°,5x,’FILE all alpha/beta values’ 
ep kL Oke 
*nd roots’ ,6x,°|~) 





0 xe 6 ‘,2x,°|°,5x,’°FILE all alpha/beta values a 


format(1x,///,15x,’Enter integer number for selection ===> ” ) 

call PSIT(19, 56) 

read(*,*,err=20) OUT 

call CLR 

if €(OUT .lt. 4)>2or. (OUT ot. Oo eOrvomse 

write(*,101) 

format(1x,////,10x,’What is the file name (fn) and extension (ext) 


*?2 “ //,10x,’Enter in form fn.ext (.ext 18 optional) ===> 7) 


callersll(s +55) 
read(*,102) FILENAME 
format (al2) 
inquire(file=FILENAME, exist=EXIST) 
if (EXIST) then 
write(*, 1030) 
format(///,10x,’That file already exists.’ ,/,10x, 
“Do you want to overwrite? (Y/N)”) 
read(*,*(A1l)’) ANS 
if( (ANS .eq. “Y’) .or. (ANS .eq. “y’)) then 
open(8,file=FILENAME, status=’ OLD’ ) 
rewind 8 
goto 30 
else 
goto 29 
end if 
end if 
open(8&’,file=FILENAME, status=’ NEW’ ) 
write(8,196) WNMIN,WNMAX 
format(1x,///,10x,“WNMIN =’ ,g11.4,9x,”“WNMAX =’ 
se 
call CR 
write(*,31) 
Write. 32) 
Write (* 33) 
write(*,32) 
write(*, 34) 
write(*,37) 
write(*,35) 
write(*,37) 
write(*, 36) 
write(*732) 
write(*, 38) 
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o format(1x,//,10x, ° ------------------------------------ --- --------- 


2 format (10x, ° ---------------------------------~---- ----------------- 


33 format (10x, ’ 





mrex, OPTION NO.” , 2x, ~ 





lox “CURVE “SELECT LON” 
































el6ex, |) 

34 formavirox,~ |~,2x,” li ere ox, Constant ZETA Curves’ 
GEA Ap Cal al 
ee OK OX,” 2 7, 2x,°|°,5x,”’Constant OMEGA Curves’ 
ere xX | 
ee LOX, |, 2x, * 3 ‘,2x,’|°,5x,’Constant ZETA*OMEGA Curves’ 
eeerGxX, © |” 
oy rOxe |” 2x,” 4 ‘,2x,’|°,5x,’°Change Frequency Range’ 
mee2Ux, |) 

35 morma tiiOx,~ |... 2x5 5 ‘,2x,°|°,5x,°PLOT Selected Curves’ 
Seioxk. |) 

36 format (1Ox.° |° .2x, ” 9 2%. i 5x nT he bo Mary Meru 
peo, Il) 





37 format (LOX , ~ FEEEEEEEEEEEEEEA EE EEFHEEEEAEEFEFEEEAFHEA +E FE H+ +4 44444444 
SSP Pim amas 
38 format(1x,///,15x,’Enter integer number for selection ===> ” ) 
call PSIT(21 ,.56) 
read(*,*,err=30) SEL3 
FLAG = 0 
call CLR 


ao SCONSEANT ZETA PLOTS *** 


GGG ©) 


ae) 


330 


40 


41 


if (SEL3 .eq. 
if (FRQCHG 


1) then 
.eq. 1) go to 380 


Constant ZETA curve selection menu 


eal CLR 

write(*,41) 
write(*,42) 
write(*,43) 
write(*,42) 
write(*,44) 
write(*,46) 
write(*,45) 
write(*,42) 
write(* ,47) 


format (1x, /7/,10x,- 


ems mms mh cm ms ym ee ee a se SS Se 
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42 


43 


44 


45 


46 


47 


format(10x,’ |’ ,2x,°’ OPTION NO.’ ,2x,°|’,14x,° ZETA CURVE SELECTION’ 
eS eae |e) 

format (10xeel* 2x. 1 7 ,2x,°|°,5x,°Select particular consta 
*¥nt ZETA Curves 4x. | 




















+f VOX 2s 2 *,2x,°|* ,5x,°ZETA = 0 Curve =27ae 
Xow), LOxG alone 3) ‘,2x,°|°,5x,° ZETA = 0,0.5 and 1.0 Gives 
oah5 x , < || 

say al Op || es 8 + “,2x,°|°,5x,° ZETA = 0,0.25,0.5) 02 foeeeeee 
*urves’ ,9x,’|*) 

format ClOx, |) exe 8) °,2x,°|’,5x,°EXIT to Curve Selection 


*Menu’ ,14x,°|7) 

format (LOX, ~ FAFHEEFEEFHEEHEEEEEEFEFE FEE EEE FEE EE EP HEE A+++ 4444444444 
pera area ieraee | 

format(1x,///,15x,’Enter integer number for Selection ===> ” ) 

Gadi) PSIT (195 56) 

read(*,*,err=30) ZTASEL 


*** CONSTANT ZETA CURVE SELECTION *** 
FRQCHG = O 
if (ZTASEL .eq. 2) then 
NUMZTA = 1 
ZETA(1) = O. 
go to 380 
else if (ZTASEL .eq. 3) then 
NUMZTA = 3 
ZETA(1) 
ZETA(2) 
ZETA(3) 
go to 380 
else if (ZTASEL .eq. 4) then 
NUMZTA = 5 
ZETA(1) = O. 
ZETA(2) = .25 
ZETA(3) = .5 
ZETA(4) = .75 
ZETA(5) = 1. 
go to 380 
else if (ZTASEL .eq. 9) then 
go to 30 
end if 


©; 
.5 
le 


call CLR 
write(*,331) 


100 


331 fOrmac(1x.////7/.10x. This is®@the constant ZETA® curve rout ine’ 


* ,/,10x, ’--------------------- =s==—-------------- ‘) 
C *** Enter Number of Constant ZETA Curves *** 
write(*,341) 
write(*,349) 
341 format(1x,////,10x,°’How many constant ZETA curves do you wish to 
* plot? ”,/) 
349 format(10x,’Enter integer value less than or equal to 10 ==> ”) 
caller siPcl4 60) 
read(*,*) NUMZTA 
Cail CER 
write(*,342) 
342 format(1x,///,5x,’Enter constant ZETA values in the range: 0 <= 


* ZETA <= 1.0’ ,/, 5x, ‘---———------- = == -- -- = - = ----------- ++ 


ROW2 = 6 
do 377 I = 1,NUMZIA,1 
write(*,343) I 
343 forma tielewy7llOm ° ZEERAGeI2p) = *) 
ROW2 = ROW2 + 3 
call PSIT(ROW2, 23) 
read(*,*) ZETA(I) 
377 continue 


380 call CLR 


FRQCHG = 0 
MAXLOG = alog10(WNMAX) 
MINLOG = alog10(WNMIN) 
INCLOG = (MAXLOG - MINLOG)/99. 
START = -99 
STOP = 0 


DO 308 I=1,NUMZTA 
START = START + 100 
STOP = STOP + 100 
if ((START.eq.1).and.((OUT.eq.1).or.(OUT.eq.5).or. 
(OUT.eq.6))) write (*,888) 
888 fommat (////#////07 // 25x * * *ee Computing Data * * *’) 
if ((OUT .ge. 2) .and. (OUT .le. 4)) then 
write(*,326) I,ZETA(I) 
326 monmatc(ly./, \Ome“ ZBMA(° ,12,°) =" ,g11.4,//) 
end Tf 
if ((OUT .ge. 4) .and. (OUT .le. 6)) then 
write(8,326) 1,ZETA(I) 
endif 


10f 


328 R = QO. 
WNMIN2 
DO 306 
Bl = 
B2 = 
Cl = 
C2 = 
Dal 
D2 = O. 
DO 303 N=1,NUMCOF 
K = N-1 
TF (K) 302,301 302 


WNMIN 
= START, STOP 


oO. oC OO Oo 


301 


302 


303 


800 


304 


U = 0.0 
Ul = =120 
U2e = 2. Ot ZETA U_UL 


D1 = (-1.0) **K*CONST(N) *WNMIN2Z**K*UL+D1 
D2 = (-1.0)**K*CONST (N) *WNMIN2**K*U+D2 
C1 = (=1.,0) **K*BETA(N) *WNMEN2Z**K Ul 
C2 = (-1.0) **K*BETA(N) *WNMIN2=*K*tU+CZ 
Bl = (-1.0)**K*ALPHA(N) *WNMIN2**K*U1+B1 
B2 = (-1.0) **K*ALPHA (N) *WNMIN2**K*U+B2 
Ul = U 
U = U2 

CONTINUE 

Zt 0r—o 

R= ti 


check for division by - O 
if (( CABS (B1*C2—-B2*@1) jee le. 


l.e-12) .and. 
(L .ne. 1) oettnen 

A(L) = A(L-1) 

B(L) = B(L-1) 

Write (= s0C ee 

if ((OUT .ge. 4) .and. (OUT .le. 6)) then 
write (8,800) L 

end if 


format (5x,°’Denominator approximately zero at point “,1I3 
,/,5x%,’° This point now assigned previous point value’ ) 


goto 802 
else 
A(L) = (C1*D2—-C2*D1) /(B1*C2-B2*Cl1) 
B(L) = (B2*D1—-B1*D2)7 (817@2-B2*Cl1) 
end if 


check for apparent discontinuity - graphics can’t handle it 


{02 


e) 


801 


802 


305 


325 


307 


309 


306 


381 


308 


i see. ) then 
if ((abs(A(L) -— A(L-1)) .gt. 1.e6) 
.or. (abs(B(L) —- B(L-1)) .gt. 1.e6)) then 
A(L) A(L-1) 
B(L) B(L-1) 
whate (*, 801) L 
if ((OUT .ge. 4) .and. (OUT .le. 6)) then 
write (8,801) L 
end if 
format (5x,’°Apparent Discontinuity at point “,I3, 
/,5x,°This point now assigned previous point value’ ) 
end If 
end if 


if (OUT .eq. 1) goto 309 
it (Coun .eq, 859m.or. (OUT .eque6)o goto 307 
if (amod(R,10.)) 307,305,307 
Wire t™,oco) Gy, ACL) ,Li,B (CE) 
fommact(lix,° A({ ,13,°) = “#e@e4,9ox,°B(@mI3,"°) = ’, 
g11.4) 
if (OUT .eq. 2) goto 309 
if ((OUT .ge. 4) .and. (OUT .le. 6)) then 
wrave(s.3o209) L,A(L),L,8(L) 
endif 
if (OUT .eq. 5) goto 309 
ALPHA] = A(L) 
BETA1L = B(L) 
call CRVOUT (OUT,L,A(L),B(L) ,ALPHA1, BETA1,WN,SEL3, 
MINIEOGy INCLOG, ING], Bl B27! ,C2,D1,D2) 
call ROOTS (ORDER,CONST, ALPHA, BETA,NUMCOF,WNMIN, 
WNMAX , ALPHA1 , BETA1 , FLAGRT , OUT, R) 


continue 
WNMIN2 = 10.**(MINLOG + R*INCLOG) 
continue 
if( (OUT .ge. 2) .and. (OUT .le. 4)) then 
write(*, 381) 
format (1x,///) 
PAUSE 
call CLR 
end if 


continue 


Po eS CONSTANIAOMEGA PLOTS *** 
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else if (SEL3 .eq. 2) then 
write(*,431) 
431 format(1x,////,10x,’This is the constant OMEGA curve routine’ 


*** Enter Number of Constant OMEGA Curves *** 
write(*,441) 
write(*,449) 


441 format(1x,////,10x,”’How many constant OMEGA curves do you wish t 
*o 1p lOtn ee 
449 format(10x,’Enter integer value less than or equal to 10 ==> “”) 


call PSIT(14,60) 
read(*,*) NUMWN 


call CLR 
write(* ,442) 
442 format(1x,///,5x,’Enter constant OMEGA values’ ,/, 
x 5x SS SS SO SSO SSeS SS SOO SS SSS es 
ROW2 = 6 


do 477 I = 1,NUMWN,1 
write(*,443) I 
443 format(1x,//, 10x, OMEGA 42) = —. 
ROW2 = ROW2 + 3 
call PSIT(ROW2, 24) 
read(*,*) WN(I) 
477 continue | 


480 call Cie 
INC1 = 1/99. 
STARTO = -99 
STOPO = 0 
DO 408 I=1,NUMWN 
STARTO = STARTO + 100 
STOPO = STOPO + 100 
if ((STARTO.eq.1).and.((OUT.eq.1).or.(OUT.eq.5).or. 
* (OUT.eq.6))) write (*, 888) 
if ((OUT .ge. 2) .and. (OUT .le. 4)) then 
write(*,426) I,WN(I) 
426 format (1X7 , 10x, OMEGA 212), Qe= sels es 
end if 
if ((OUT .ge. 4) .and. (OUT .le. 6)) then 
write(8,426) I,WN(I) 
endif 
R = OQ. 
ZE LAI 207 
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DO 406 L = STARTO,STOPO 
Bl = 0 
B2 = QO. 
C1 = QO. 
C2 = 0 
Di= 0. 
D2 = C0. 
DO 403 N=1,NUMCOF 
Ke= N-1 
IF (K) 402,401,402 
401 U = 0.0 
Ul = =O 
402 UZ ee .OSZETAI*U-Ul 
Dl = (= 190) * *K* CONST (N) *WN(1I) **K*U14D1 
BZ = (—2.0) + * KR CONST (CN) AWN (1) ** K*U+D2Z 
Cie = 107 BETA (N) *WNCT) **K*U1+Cl 
C2= (=170) 79h? BEPA(N ) *WN( 1) **K*U+C2 
Bie = (=1..0) **KYAPr eA (N)*WN(T) **K*U14+B1 
Be = (S10) * RF ALPHACNY WN (1) * *K*U+B2 
Ul = U 
U = U2 
403 CONTINUE 
Ro R+L 


check for division by - 0 
if ((ABS(B1*C2-B2*C1)) .le. 1.e-12) .and. 
(L .ne. 1)) then 
AOMEGA(L) = AOMEGA (L-1) 
BOMEGA(L) = BOMEGA(L-1) 
write (*,800) L 
mf ((OUZ “ee. 4) .and. (OUT .le. 6)) then 
write (8,800) L 
end if 
goto 803 
else 
404 AOMEGA (L) 
BOMEGA (L) 
end if 


(2102-02 D1) /(Bi*C2-B2er} 
(B2*DiI-B51*D2)7 (B1*C2-B2+Cl1) 


check for apparent discontinuity - graphics can’t handle it 
if ((abs(AOMEGA(L) - AOMEGA(L-1)) .gt. 1.e6) .or. 


“i (abs (BOMEGA(L) - BOMEGA(L-1)) .gt. 1.e6)) then 
AOMEGA(L) = AOMEGA(L-1) 
BOMEGA(L) = BOMEGA(L-1) 
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‘@) 


803 


405 
425 


407 


409 


406 


408 


531 


941 


549 


write (*, 801) er 
if ((OUT .ge. 4) .and. (OUT .le. 6)) then 
write (8,801) L 
end if 
end if 


if (OUT .eq. 1) goto 409 
if ((OUT .eq. 5) sor. (OUT .eqmeG)) Foto raen 
if (amod(R,10.)) 407,405,407 
write(*,425) L,AOMEGA(L) ,L,BOMEGA(L) 
format (11x, “AOMEGA(’ ,I3,°) = *,g11.4,9x, 
“BOMEGA(’ ,I3,’) = °,gl11.4) | 
if (OUT .eq. 2) goto 409 
if ((OUT .ge. 4) .and. (OUT .le. 6)) then 
write(8&,425) L,AOMEGA(L) ,L,BOMEGA(L) 
endif 
if (OUT .eq. 5) goto 409 
ALPHA1L = AOMEGA(L) 
BETA1 = BOMEGA(L) 
call ROOTS (ORDER,CONST, ALPHA, BETA, NUMCOF ,WNMIN, 
WNMAX , ALPHA] , BETA1 , FLAGRT , OUT, R) 


continue 
ZETAL = ZETAI1 + INCI 
continue 
if( (OUT .ge. 2) .and. (OUT .le. 4)) then 
write(*,381) 
PAUSE 
call CLR 
end if 


continue 
*** CONSTANT ZETA*OMEGA PLOTS *** 
else if (SEL3 .eq. 3) then 


write(*,531) 
format(1x,////,10x,’This is the constant ZETA*OMEGA curve routin 


*e / 4 yi 9 10x ’ ‘ oe re re re re re re oe re re re ee ee ee ee SS SS SS SS SS SS SS SS SS a ¢ ) 


*** Enter Number of Constant ZETA*OMEGA Curves *** 
write(*,541) 
write(*,549) 
format (1x,////,10x,”’How many constant ZETA*OMEGA curves do you w 


*ish to plot? ”,/) 


format(10x,’Enter integer value less than or equal to 10 ==> ”) 
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C 


542 


543 


o77 


580 


526 


501 


call PSIT(14,60) 

read(*,*)} NUMZWN 

Call CLR 

write(*,542) 

format(1x,///,5x,’°Enter constant ZETA*OMEGA values’ ,/, 


BX , a a 
ROW2 = 6 


do 577 I = 1,NUMZWN 


write(*,543) I 

fOrmat (LX, 7 /ehOx, ZETA*OMEGA(” , 12, * ) 
ROW2 = ROW2 + 3 

call PSIT(ROW2, 29) 

read(*,*) ZWN(I) 


ll 
. 
— 


continue 
calla CLR 
MAXLOG = alog10(WNMAX) 
MINLOG = alog10O(WNMIN) 
INCLOG = (MAXLOG —- MINLOG)/99. 
START = -99 
STOPZO = 0 


DO 508 I=1,NUMZWN 
START = START + 100 
STOPZO = STOPZO + 100 
if ((START.eq.1).and.((OUT.eq.1).or.(OUT.eq.5).or. 
(OUT.eq.6))) write (*, 888) 
if ((OUT .ge. 2) .and. (OUT .le. 4)) then 
write(*,526) I,ZWN(I) 
format(1X,/,10x, ° ZETA*OMEGA(’ ,I2,°) =" ,g11.4,//) 
end if 
if ((OUT .ge. 4) .and. (OUT .le. 6)) then 
write(8,526) I,ZWN(I) 
endif 
Ra OC. 
WNMIN2 
DO 506 L 
Bl = 0 
B2 = QO. 
Ci = "0: 
0) 
0 


WNMIN 
= START, STOPZO 


C2 = 
Dix= 9. 
D2 = O. 
DO 503 N=1,NUMCOF 
K = N-1 
[heres 2. 501, 502 
U = 0.0 
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Ul = -1.0/WNMIN2**2. 

502 U2 = —-2.0*ZWN(1I)*U — WNMIN2**2*U1 
D1 = CONST(N)*U1+D1 
D2 = CONST(N) *U+D2 
Cl = BETA(N)*U1+C1l 
C2 = BETA(N) *U+C2 
Bl = ALPHA(N)*U1+B1 
B2 = ALPHA(N) *U+B2 
Ul =U 
U = U2 

503 CONTINUE 

R = R+1. 


check for division by ~ 0 
if ((ABS(B1*C2-B2*C1l)) .le. 1.e-12) .and. 
(L .ne. 1)) then 
AZTAOM(L) = AZTAOM(L-1) 
BZTAOM(L) = BZTAOM(L-1) 
write (*,800) L 
if ((OUT .ge. 4) .and. (OUT .le. 6)) then 
write (8,800) L 
eng: if. 
goto 804 
else 
504 AZTAOM(L) (C1*D2-C2*D1) / (B1*C2-B2*C1) 
BZTAOM(L) = (B2*D1—-B1*D2) / (B1*C2-B2*Cl1 ) 
end if 


check for apparent discontinuity - graphics can’t handle it 
if ((abS(AZTAOM(L) - AZTAOM(L-1)) .gt. 1.e6) .or. 
. (abs(BZTAOM(L) - BZTAOM(L-1)) .gt. 1.e6)) then 
AZTAOM(L) AZTAOM(L-1) 
BZTAOM (L) BZTAOM(L-1) 
write 7478015) 
if ((OUT .ge. 4) .and. (OUT .le. 6)) then 
write (8,801) L 
end if 
end if 


804 if (OUT .eq. 1) goto 509 
if ((OUT .eq. 5) vor. {OUT -eq. 6)) goto 507 
if (amod(R, 10.) 5077 505,507 


505 write(*,525) L,AZTAOM(L) ,L,BZTAOM(L) 
525 format (11x, °“AZTAOM(~ ,13,°) = ~ ,g11.4,9x, 
* “BZTAOM(’ ,I13,°) = °g1l1.4) 


5Od if (OUT .eq. 2) goto 509 
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909 


906 


908 


297 


298 


299 


296 


295 


Poet tOuUswee. 4) .and. (OUT .le. 6)) then 
write(8,525) L,AZTAOM(L) ,L,BZTAOM(L) 

endif 

if (OUT .eq. 5) goto 509 

ALPHA1 = AZTAOM(L) 

BETA] = BZTAOM(L) 

call ROOTS (ORDER,CONST, ALPHA, BETA,NUMCOF,WNMIN, 

WNMAX , ALPHA1 , BETA] , FLAGRT , OUT, R) 


continue 
WNMIN2 = 10.**(MINLOG + R*¥INCLOG) 
continue 
if( (OUT .ge. 2) .and. (OUT .le. 4)) then 
write(*,381) 
PAUSE 
call CLR 
end if 


continue 


*** FREQUENCY RANGE CHANGE FOR CONSTANT ZETA PLOT *** 
else if (SEL3 .eq. 4) then 


*** Enter RANGE of frequencies to examine *** 
write(*,296) WNMIN,WNMAX 


write(*,297) 
format(1x,//,5x,’°Enter the RANGE of frequencies you wish to exam 


*ine’ ,/,5x, ° -------------------------- - - - - - ae a 


write(*,298) 
format (1x,/,5x,’°Enter minimum frequency (WnMIN) ===> ”) 
call PSIT(13,43) 
read(*,*) WNMIN 
write(*,299) 
format(1x,//,5x, “Enter maximum frequency (WnMAX) ===> “” ) 
call PSIT(16, 43) 
read(*,*) WNMAX 
write(*,296) WNMIN,WNMAX 
format(1x,////,10x, “WNMIN =’ ,g11.4,9x,”“WNMAX =’ ,g11.4) 
if (FILENAME .ne. “ “”) then 
write(8,295) WNMIN,WNMAX 
format (1x,///,10x,”WNMIN =’ ,g11.4,9x,”“WNMAX =’ 
pears, / / } 
endif 
FRQCHG = 1 
SEue =| | 
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QO 


399 


write(*,381) 
PAUSE 


*** CALL PLOTTING ROUTINE *** 
else if (SEL3 .eq. 5) then 


SEL = 7 

FLAG = 1 
else 

continue 
end if 


Recompute Alpha & Beta points in constant zeta routine when 
frequency range is changed 
if(FRQCHG .eq. 1) go to 330 


if(SEL3 .eq. 5) go to 399 


If EXIT not selected, redisplay curve select menu 
if (SEL3 .ne. 9) go to 30 


if (FILENAME .ne. “ 7”) then 
close(8,STATUS=’ KEEP’ ) 


endif 


return 
end 
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ROOT FINDING MODULE 


KOK OR KK KK OK OK OK OK OK ok KK OK ok OK oR OK OK ok OK OK ok OK ok OK OK OK OK OK OK OK KK OK OK OK OK OK OK OK OK OK OK OK OK OK ok OK OK ok ok OK OK OK ok OK OK OK OK OK OK OK ok OK OK 


mmm rm rm mm re ee me ei i ei 
———— oe ce a cm cr cr ce ec ee ee ee ee ee ee ee ee ee ee ee ee eee oe OO eee eee ee ee 


Subroutine ROOTFD -- calculates the roots of the entered C. E. 


ee ce ems ces rr peers cops copes cep cm cc cece caper corre epee crc cr cr cr er cr cr cr cr cr cr crc cr crs cr cs per rr re cr rr re ee rr rm re re ee ee ee ee ee 
=—— OO oOo ce cs crm cr i ee re ee ee es ee ee ee ee ee ee ee ee eee ee ee ee ee ee SS ee eee ee 


Subroutine ROOTS (ORDER,CONST, ALPHA, BETA ,NUMCOF ,WNMIN ,WNMAX, 


* ALPHAI1 , BETA] , FLAGRT , OUT , R) 
real CONST(10) ,ALPHA(10) ,BETA(10) ,COEF(13),A,ALPHA1,BETA1,R 
real*8 B,C,P,Q@,DENOM, REAL1 , IMAG1,REAL2, IMAG2 , IMAGSQ, 
= REALNO , IMAGNO 
integer ORDER,I,NUMCOF, IMIN1, FLAG, FLAGRT, OUT 
character*1l CHG 
character*12 FILENAME 
common/filenam/ FILENAME 
data FLAG /0/ 


if (FLAGRT .eq. 1) then 
call CHAREQ (ORDER,CONST, ALPHA, BETA , NUMCOF , WNMIN , WNMAX, 
* FLAGRT , ALPHA1 , BETA1 ) 
end if 


call NORMAL (ORDER , NUMCOF , CONST , ALPHA, BETA, COEF ,ALPHA1,BETA1, 


7 FLAGRT ) 
call ROOTS2(ORDER,COEF,FLAG,P,Q,OUT,R) 
return 
end 
Subroutine ROOTS2 -- Reduce polynomial uSing Bairstow’s method 


i si i Se 
ee ee SSS SS 


Subroutine ROOTS2 (ORDER, COEF,FLAG,P,Q@,OUT,R) 
dimension F(13) ,G(13) 

real EPSLON , DELTAP, DELTAQ,COEF(13) ,R,REDUC,PREVDP, 
- PREVDQ 

real*8 DENOM,P,Q,REAL1,IMAG1,REAL2,IMAG2,F,G 

integer J,ORDER,ORD,ORD3,N,ITERAT,MAXIT, FLAG, ORD1 ,OUT 
character*12 FILENAME 

logical PDIFF,QDIFF,PQCHG 

common/filenam/ FILENAME 

EPSLON = .0001 


tit 


'?) 
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TO% 


100 


MAXIT = 500 

ORD = ORDER 

ORD3 = ORDER + 3 
ITERAT = 1 

FLAG 0 

P = 
oa 
Pa 
F(2y 
G(1) 
G(2) = 
DELTAP = 
DELTAQ 


[ee I 
oO OCG 


ofl 
[iS e) 


Check for roots at origin 


l-=70 
if (COEF(ORD3) .eq. 0.) then 
a ae Oe Cae 
ORD3 = ORD3 - 1 
ORD = ORD - 1 
if ((COEF(ORD3) .eq. O.) goto 108 
if € (°C OUT .eqme3s ) .or. ( OUT .eq. 4) ) ame 
( amod(R,10.) .eq. O. ) then 
write (*,107) I 
write (8,*). ~~ 
write (8,107) I 
format(10x,I1,° root(s) at the origin } 
end if 
end if 


if (ORDER .eq. 1) then 
REAL1 —COEF (4) 
IMAG1 Oe 
if ( ¢( ¢€ OUT .eqm3 ) .or. ( OUT 7eq. 4 ) )) and. 
( amod(R,10.) .eq. O. ) then 
write (*,100) REAL1, dabs(IMAG1) 
format(1x,/, 10x, Only Reot)—) Sxne7 4 ee, a , cil ee 
end if 
if ((OUT .eq. 4) .or. (OUT .eq. 6)) then 
write (8,100) REAL1, dabs(IMAG1) 
endif 
FLAG = 1 


Il 


else if (ORDER .eq. 2) then 
P = COEF (4) 
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On eCOnF (o) 
call ROOTS1(P,Q,REALI1, IMAG1,REAL2, IMAG2) 
PamGcCOUL sedues)  .Or., (OUL .eq, 4)) .ands = (amod(R,10.) .eq. 0.) 
*) then 
write (*,105) REAL1, dabs(IMAG1) 
write (*,106) REAL2, dabs(IMAG1) 
WEIte (47) ~ 


105 foarnagaelOoxsarirst RoOt = ~° ,5x,¢g9.3,° + 3° ,g@9.3) 
106 format(10x,’Second Root = ” ,5x,g9.3,° - 3° ,g9.3) 
end if 


mf C(OUT .eq. 4) .or. (OUT .eq. 6)) then 
write (8,105) REAL1, dabs(IMAG1) 

write (8,106) REAL2, dabs(IMAG1) 

wrasee (8, *)m° ” 

endif 

FLAG = 1 


end if 
If all roots computed or max iterations exceeded — finished 
10 if ((FLAG .eq. 1) .or. (ITERAT .gt. MAXIT)) go to 99 
Gonl10.) = 85 ORDS 
F (J) COEF(J) - P*F(J-1) -— Q*F(J-2) 
G(J) = F(J) - P*G(J-1) - Q*G(J-2) 


110 continue 
DENOM = G(ORD+1)**2 — (G(ORD+2)—-F(ORD+2) )*G( ORD) 


if (DENOM .ne. QO.) then 


PREVDP = DELTAP 
PREVDQ = DELTAQ 
DELTAP = (F(ORD+2)*G(ORD+1) — F(ORD+3) *G(ORD) ) *REDUC/DENOM 
DEEPAQ = (GCORDTI)*F(ORD+3) — (G(ORD+2) —F(ORD#2) ) * 
* F (ORD+2) ) *REDUC /DENOM 
Pease et DEE TAP 
Q = Q@ + DELTAQ 


if ((abs(P) .gt. EPSLON) .and. (abs(Q) .gt. EPSLON)) then 
PQCHG = (absS(DELTAP/P) + abS(DELTAQ/Q)) .1t. EPSLON 

else if ((abs(P) .1t. EPSLON) .and. (abs(Q) .gt. EPSLON)) then 
PQCHG = (abs(DELTAP) + abs(DELTAQ/Q)) .1t. EPSLON 

else if ((abs(P) .gt. EPSLON) .and. (abs(Q) .1t. EPSLON)) then 
PQCHG = (abS(DELTAP/P) + abS(DELTAQ)) .1t. EPSLON 

else 
PQCHG = ((abs(DELTAP) .1t. EPSLON) .and. 
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* (abS(DELTAQ) .1t. EPSLON) ) 
end if 
if (PQCHG) then 
call ROOTS1(P,Q,REAL1,IMAG1,REAL2, IMAG2) 
ORD1 = ORD - 1 
if (((OUT .eq. 3) .or. (OUT .eq. 4)) .and. (amod(R,10.) .eq. 


*O,) ) then 
write (*,150) ORD,REAL1,dabs(IMAG1) ,ORD1,REAL2,dabs(IMAG2) 
150 format(1x,° Root( ,12,°) =" 7ell-4) 3 fe ee 
x 5x, Root(~,12,°) =" 7ell 4,” =e) ela 
end if 


if ((OUT .eq. 4) .or. (OUT .eq. 6)) then 

write (8,150) ORD,REAL1,dabs(IMAG1) ,ORD1,REAL2,dabs(IMAG2) 
endif 
ORD = ORD - 2 


Check order of reduced polynomial 
if (ORD .eq. 0) then 


FLAG = 1 
else if (ORD .eq. 1) then 
REAL1 = -F(ORD+3)/F (ORD+2) 
IMAG] = QO. 
if (( (OUT .eq. 3) .or. (OUT .eq. 4)) .and. (amod(R,10.)ieeu 
* O.)) then 
write (*,160) ORD,REAL1, IMAG1 
160 format(1x,“Reot(’ ,12, ) = ,e11.4,° 45 {cli 
end if 


if ((OUT .eq. 4) .or. (OUT .eq. 6)) then 
write (8,160) ORD, REAL1, IMAG1 

endif 

FLAG = 1 

else 

ORD3 = ORD + 3 

qo 18050 = "c,ORDe 
COEF(J) = F(J) 


TPERAT —2e 
180 continue 
end if 
else 
ITERAT = ITERAT + 1 
end if 
else 
P= P+ 1 


Q=aQ +1 
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ITERAT = 1 


end if 
PDIFF = (abs(PREVDP) + abs(DELTAP)) .ne. absS(PREVDP + DELTAP) 
QDIFF = (abs(PREVDQ) + abs(DELTAQ)) .ne. abs(PREVDQ + DELTAQ) 


Pam(eDlL nee anG ~ut@ eelt. .Ol)) .or. (QDIFF .and. (P .2t. .01)) 
* .or. (PDIFF .and. QDIFF)) then 
REDUGC l= REDUC/2. 
end if 
if (FLAG .ne. 1) go to 10 
99 return 
end 


em me eccrine cr cr cr re cr cr cr cr ce cr cs cr cs ce cr cc cs ce cr cs cs cs cs ce cs ee ee ee 2 oe 
a a ee ee 


Subroutine ROOTS1 -- Finds roots of second order polynomial 


me mem sem emcee se cece crm cert etre cr ccc cr i i i rt re tr te er rt i ee ee ee ee ee ee oe 
rr ere er ar 


Subroutine ROOTS1(B,C,REAL1, IMAG1,REAL2, IMAG2) 
real*8 B,C,REAL1,IMAG1,REAL2, IMAG2 ,REALNO, IMAGNO, IMAGSQ 


REALNO = -B/2. 

IMAGNO = B*¥*2 — 4*C 

if (IMAGNO .ge. 0.) then 
IMAGSQ = dsqrt(IMAGNO) /2. 


REALI = REALNO + IMAGSQ 
IMAGI = O. 
REAL2 = REALNO —- IMAGSQ@Q 
IMAG2 = O. 

else 
IMAGSQ = dsqrt(—-IMAGNO) /2. 
REALI = REALNO 
IMAG1 = IMAGSQ 
REAL2 = REALNO 
IMAG2 = -—IMAGSQ 

end if 

return 

end 


ei ee em mw mi i es SS i SS SS 
cr cr cr cr cr cr i cr cc cr cr cc ee cr a cr cr cr cr cr cm ee ee ee 


cm ww ew ww ee a a i ise a 
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Subroutine NORMAL (ORDER , NUMCOF, CONST ,ALPHA, BETA,COEF,ALPHA1,BETAI, 
* FLAGRT ) 
real CONST(10) ,ALPHA(10) ,BETA(10) ,COEF(10) ,ACOEF ,BCOEF,CCOEF, 


iis 


10 


20 


30 
25 


* 


ALPHA! BEAL 
integer J,ORDER, NUMCOF,N,M, FLAGRT 


N= 3 

M ORDER + 3 

do 10 J = NUMCOF,1,-1 
COEF(N) = CONST(J) + ALPHA(J)*ALPHA1L + BETA(J)*BETAI1 
N= Nee 

continue 

if (COEF(NUMCOF) .ne. 1.) then 
do 20 J = 4,M 

COEF (.)} = COEF (J) /COEF (3) 

continue 

COEF(3) = 1. 

end if 

if (FLAGRT .eq. 1) then 

Gal clan 

I = NUMCOF + 1 

do 25 J = 3,M 

l =e 

write (*,730) 1,COEF()) 

format (15x. COEF( 212. )  =s e207) 

continue 

end if 

return 

end 
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PLOTTING MODULE 


ORO OK KOK ROR ORO GK OK IR RK OK OR IK IIR IK I RIK ak ke ok ok 


em mc mm mc ce ce es we ww ec ce we ce we ce ww ce wm cw cme wr we ee we we we we we oe eo oe oe ee 
=———_——— ee cet SS cc ce ccc cr ce ce rm cc cc ce ee ee ee ee ee ee ee ee Se Se ee Se Se Se a See eS ee ee ee ee ee ee 


Subroutine RUNPRO -- executes the parameter plane program 


ee eee 
— Se ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ae eee ae ae a SS SS a ee ee eee 


Subroutine RUNPRO (NUMZTA, STOP, ZETA,CURV,NUMWN ,WN,STOPO, 
* ZWN , NUMZWN , STOPZO) 
real ATEMP(102) ,BTEMP(102) ,ZETA(10) ,A,B,AOMEGA, BOMEGA, 
XMIN, XMAX , YMIN, YMAX, FACT, XCEN , YCEN,WN(10) , ZWN(10), 
EXPAMN , EXPAMX , EXPBMN , EXPBMX , AZTAOM, BZTAOM, 
PeeeaALP¢ 1002) ,EXPBET (1002) ,STITLE, CRVfYP, 
AZFST,AZDLT, BZFST,BZDLT, AOFST , AODLT , BOFST , BODLT, 
AZOFST , AZODLT , BZOFST , BZODLT 
integer I,J, IOPORT,MODEL,NUMZTA, STOP,FRSTZ,DELTZ,INTRVL,CURV, 
* DELT1 , DELT2 ,NUMWN ,L,K,STOPO, FRSTW,DELTW,EXPAND,INTSRT, 
= INTSTP, NC, NUMZWN , STOPZO, SEL, SYMBL 
dimension A(1002) ,B(1002) , AOMEGA (1002) , BOMEGA(1002) , AZTAOM(1002) , 
= BZTAOM(1002) 
character*1l CHG,NMCHAR (30) 
enaracter*30 TITLE 
common A,B, AOMEGA, BOMEGA , AZTAOM, BZTAOM 
common/plottr/ IOPORT ,MODEL 
common/box/ EXPAND, SYMBL 
common/factr/ FACT 
common/symb/ SYMBL 
equivalence (NMCHAR(1),TITLE) 


* %§ «© %8% * 


FACT = O. 
EXPAND = 
SYMBL = 14 


9 
9 


call CLR 

write(*,80) 

write(*, 81) 

write(*, 82) 

write(*, 81) 

write(*, 83) 

write(*, 81) 

write(*,93) 

format(1x,///////,10x, ° -------------------- ---- ---- --------------- 
¥ mn or lOx. ||" , 24seaeP Oring MENU” 25x97 |~) 
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81 format (10x, ° ----—---—-—-——-- ~~ -—-~ = = — — === —_ === - ee 











aaa ) 
82 format (10x. "| a2 er Le NOI" 72x; | “320%. OP TEONe. 21 xem “) 
83 format (10x, | 92x, 1 fox es oie 

*’TITLE output graph and PLOT data’ ,10x,’” |’ 

x /, Ox 2 ee 2 Se a 

*“PLOT data (no title)’ ,22x,’|’ 

xf Oa oe xe 3 ok, en Ok, 

*’ SIZE output graph’ ,25x, ’ |’ 

#7] Ox alee Xe 4 Reo eae eo 

*’ SYMBOL to be plotted at each data point’ ,3x,’|’ 

* f/f , VOxPeelee 2 Xe 9 77,ex%, | ,SX@OBXIT to Main Menu’ , 25x) 
93 format(1x,///,15x,’°Enter integer number for Selection ===> ”) 


call Pslt(237 56) 
read(*,*,err=78) SEL 
if (SEL .eq. 1) then 
goto 4 
else if (SEL .eq. 2) then 
goto 9 
else if (SEL .eq. 3) then 
goto 7 
else if (SEL .eq. 4) then 
goto 5 
else if (SEL .eq. 9) then 
goto 399 
else 
goto 78 
end if 


Cc. . . Adjust size of output plot 
7 Cal lICrER 
write(*, 15) 
eS format (/////,10x, “Output plots are currently sized to fill the con 
*sole screen.’ ,//, 
* 10x,”“Would you like to adjust the plot size?’ ,/,10x, 
* “Enter ity or "n" ==> <5) 
Cal leborrc1o. 31) 
read(*,’ (A)’% ,err=4) CHG 
if ((CHG .eq. ‘’n’) .or. (CHG .eq. ’N’)) then 
go to 78 
end if 
write(*,16) 
16 format(///,10x,’A default plot factor size of 0.9 is used to fill 
*the console screen.’ ,/,10x, 
*“This plot size will be halved by entering a value of 0.45.’ ,/, 
*10x,°It will be doubled with a factor entry of 1.8.’,//,10x, 
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17 


6 


al 
13 


* 


* 


* 


“Enter factor size == “) 
call PSIT(18,32) 

read(*,*) FACT 

goto 78 


Select symbol to be plotted at each data point 
ealt> CLR 
write(*,17) 
format(/////,10x,’Type in an INTEGER number from O through 13 to p 
lace a’,/, 
10x,’symbol at each calculated data point. There are 100 data poi 


*nts* ,/, 


* 


* 


10x,’for each curve. Examples of symbols with associated numbers: 
eS 

SX, 2 Triangle’ ,/, 

ox Ss te oe 

15x, 4 ) a a 

15x,°8 Bf 

5X eo b BO ya 

he Ne 

5x, re Vertical Line’ ,///, 

10x,’“Enter INTEGER number (0 - 13) ==> ”) 
call PSIT(20, 45) 
read(*,*,err=78) SYMBL 

go to 78 


Enter TITLE of graph 
call CER 
write(*,6) 
format (/////,10x,’Enter plot title == een 


*10x,°(30 characters max)’ ) 


call’ Psil(7, ol) 
read(*,8) TITLE 
format (A30) 
do J1°P = 30,1,-1 
if (NMCHAR(I) .ne. ” ”) go to 13 
continue 
NC = I 
STITLE = 4.3 -— .12*NC 


call CLR 

INTRVL = O 

Seee— 3 

CURV = 0 

call PLOTS (0, IOPORT,MODEL) 
cal] NEWPEN (1) 
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Cs hee Calculate Constant ZETA Plot 


if (NUMZTA .ge. 1) then 
CRVTYP = .1l 
write (*,119) 
119 format (///,25x, “CONSTANT ZETA CURVES’ ) 
call PLTCRV (STOP,A,B,NUMZTA ,CRVTYP, ZETA, 


- TITLE, STITLE,NC,AZFST , AZDLT , BZFST , BZDLT) 
call PLOT (0.0,0.0,-999) 
end if 
eee Calculate Constant OMEGA Plot 
if (NUMWN .ge. 1) then 
call CLR 
CRVTYP = -.1 
write (*,219) 
219 format(///,25x,° CONSTANT OMEGA CURVES’ ) 
call PLTCRV (STOPO, AOMEGA , BOMEGA , NUMWN ,CRVTYP, WN, 
* TITLE, STITLE,NC,AOFST,AODLT, BOFST , BODLT ) 
call PLOT (0.0,0.0,-999) 
end if 


eee Calculate Constant ZETA*OMEGA Plot 


if (NUMZWN .ge. 1) then 


Cae Cin 
GRVIYE =-—-2 
write (*,319) 
319 format(///,25x, “CONSTANT ZETA*OMEGA CURVES’ ) 

call PLTCRV (STOPZO, AZTAOM, BZTAOM, NUMZWN ,CRVTYP, ZWN, 

~ TITLE, STITLE,NC,AZOFST , AZODLT , BZOFST , BZODLT ) 
call PLOT (0.0,0.0,-999) 

end if 


Calculate Zeta & Omega Plot 


if ((NUMZTA .ge. 1) .and. (NUMWN .ge. 1)) then 


eal PeCinn 

ATEMP(101) = AZFST 
ATEMP (102) = AZDLT 
BTEMP(101) = BZFST 
BTEMP(102) = BZDLT 


write (*,492) 
492 format(///,25x,”’CONSTANT COMBINED CURVES’ ) 
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write (*,493) ATEMP(101) ,ATEMP(102) 
493 format (///,5x,°X AxiS Min Value = ” ,g11.4,5x,’X Axis DELTA = ’”, 
‘ pi a /) 
write (*,494) BTEMP(101) , BTEMP(102) 
A204 fonmat (5x,°Y Axis Min Value = “ ,@11.3,5x)°Y Axis DELTA = ’ ,gl1l1.4, 
s Vi TT) 
write (*,495) 
495 format (25x,’*** Calculating Data ***’ ,////) 


INTRVL = O 

call FACTOR (FACT) 

call ASPECT (1.2) 

Colle SOTAXLS (.13,.20,.15,0.1,2) 

€atlt SeMBOL (STITLE,6.0,.25,TITLE,0O. ,NC) 

call AXIS (.8,.8,’Alpha’ ,-5,-7.0,0. ,AZFST,AZDLT) 
call AXIS (.8,.8,’Beta’ ,4,-5.0,90. ,BZFST,BZDLT) 


C 
C .. . Draw dashed lines at axis increments and surrounding box 
call GRID 
C 
do 470 I = 1,NUMWN 
do 480 J = 1,100 
ATEMP(J) = AOMEGA(J + INTRVL) 
BTEMP(J) = BOMEGA(J + INTRVL) 
480 continue 


call CURVE (ATEMP,BTEMP,100,-.1) 
call WHERE (X,Y,FACT) 
call NUMBER (X,Y,0.15,WN(I),0. ,2) 
INTRVL = INTRVL + 100 
470 continue 
c 
INTRVL = 0 
do 440 I = 1,NUMZTA 
do 450 <i2="17100 


ATEMP(J) = A(J + INTRVL) 
BTEMP(J) = B(J + INTRVL) 
450 continue 


call CURVE (ATEMP,BTEMP,100, .1) 
Call WHERE (X,Y,FACT) 
call NUMBER (X,Y,0.15,ZETA(I),0O. ,2) 
INTRVL = INTRVL + 100 

440 continue 


if (NUMZWN .ge. 1) then 


INTRVL = 0 
do 540 I = 1,NUMZWN 


be 


OQ 


550 


940 


30 


31 


32 


33 


34 


35 


36 


37 
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pa esa nee aa “,/;10x,° |", 19x, EXPAND PLOT MENU aia 


do. 550 Je 400 


ATEMP(J) = AZTAOM(J + INTRVL) 
BTEMP(J) = BZTAOM(J + INTRVL) 
continue 

call CURVE (ATEMP,BTEMP,100,-.2) 
call WHERE (X,Y,FACT) 
call NUMBER (X,Y,0.15,ZWN(I),0.,2) 
INTRVL = ENTRVE + 100 

continue 

end if 

end if 


call PLOT (0.0,0.0,999) 


Expand plot about a selected area or point 


cal cer 

write(*,31) 

write(*, 32) 

write(*,33) 

WEHlUte(. 32) 

write(*,34) 

write(*,36) 

write(*,35) 

write(*, 32) 

write(*,37) 

format(1x,// ,10%, =———_ = ee 
“*) 


a ee el 











formacvwo< 
* 14x eS) 

format(10x, * 
*xes values’ , 8x,’ 


“7, 2x, “OPTION NO.’ ,2x,°|° ,14x, “EXPANSION SELECEIONs 














" Veen al Cooked 


4 


°,5x,°Expand area defined by a 

















+ LO xen OX, 2 7 ,2x,°|°’,5x,°’Expand around a selected poi 
bap @ Gales) > cna) Wee 

format tlon, || .2x,. 4 °,2x,°|°,5x%, EXIT plotting rowtmes 
+2 eee) 


format (10x, 
*+++++++4+ | ~ ) 
format(1x,///,15x,’Enter integer number for selection ===> ”) 
call PSIT(17,56) 

read(*,*,err=30) EXPAND 





FEFEFFFELFFE FALE FAEFEFEAFEFEFEFE FFF FEE FFE EAFEF EEF EF EF EF 


if (EXPAND .eq. 2) then 


22 
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297 


298 


300 


Expand plot around selected point 


Carl CER 

write(*,295) 

format(1x,//,5x,”’Enter the POINT about which you wish to expand’ , 
*/ 5x, ° -------- - - - - - - - - - esx. (Tnis 
*point must be entered in # of divisions on alpha & beta’ ,/,5x,’axe 
*s. An expansion about the center of the diagram would ”’,/,5x,’requ 
*ire an entry of 3.5 for the alpha axis and 2.5 for the ’,/,5x,’bet 
Pa axis.)° ,//) 

write(*,297) 

format(1x,//,5x,’Enter CENTER value for ALPHA axis (range is 0 to 
*7) ===> “) 

call PSIT(14,63) 

read(*,*) XCEN 

write(*,298) 

format(1x,//,5x,’Enter CENTER value for BETA axis (range is 0 to 5 
*) ===> “) 

cals BSIT (17,62) 

read(*,*) YCEN 

write(*,300) 

format(1x,//,5x,’Enter value for FACTOR (enlargement) ===> ”) 

call PSIT(20,50) 

read(*,*) FACT 


Compute axis scaling values 


XCEN ACEN * FACT 
YCEN = YCEN * FACT 


XMIN = XCEN —- 3.5 
XMAX = XCEN + 3.5 
YMIN = YCEN - 2.5 
YMAX = YCEN + 2.5 


if (NUMZTA .eq. 0) then 
ATEMP(101) = AOFST 
ATEMP(102) = AODLT 
BTEMP(101) = BOFST 
BTEMP(102) = BODLT 

else 
ATEMP(101) = AZFST 
ATEMP(102) = AZDLT 
BTEMP(101) = BZFST 
BTEMP(102) = BZDLT 

end if 
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call CLR 
write (*,392) ATEMP(101) ,ATEMP(102) 

392 format (/////,5%,°X Axis Min Value = “ ,g11.4,5x,°X Axis DELTA 
‘i gl .4,/77) 
write (*,393) BTEMP(101) , BTEMP(102) 

393 format (5x,°Y Axis Min Value = ” ,g11.3,5x,°Y Axis DELTA =) (2m 
= IIL), 


4 


call »sPLOTS (0, IOPORT,MODEL) 
call WINDOW (XMIN, YMIN, XMAX, YMAX) 
call FACTOR (FACT) 
call ASPECia(!). 2) 
cally STAXIS (137.20. .15.043,2) 
call “SYMBOL (STITLE, 6.05425, 711LE. 0. NG) 
call AXIS (.8,.8,’Alpha’ ,-5,-7.0,0.,ATEMP(101) ,ATEMP(102) ) 
call AXIS (.8,.8,’Beta’ ,4,-5.0,90. ,BTEMP(101) ,BTEMP(102) ) 
INTRVL = O 
Wiebe (Soo) 
331 format (/7//7,25x, *7** Calculating Data, ***... 


Draw dashed lines at axisS increments and surrounding box 
call GRID 


do 340 I = 1,NUMZTA 
do 350 J = 1,100 


ATEMP(J) = A(J + INTRVL) 
BTEMP(J) = B(J + INTRVL) 
350 continue 


call CURVE (ATEMP,BTEMP,100,.1) 

call WHERE (X,Y,FACT) 

call NUMBER (X,Y,0.15,ZETAQ) 0.2) 

if ((SYMBL .le. 13) .and. (SYMBL .ge. 0O)) then 
call LINE (ATEMP,BTEMP,100,1,-1,SYMBL) 


end if 
INTRVL = INTRVL + 100 
340 continue 
INTRVL = O 
do 370 I = 1,NUMWN 
do 380 J = 1,100 
ATEMP(J) = AOMEGA(J + INTRVL) 


BTEMP(J) = BOMEGA(J + INTRVL) 
380 continue 
call CURVE (ATEMP,BTEMP,100,-.1) 
call WHERE (X,Y,FACT) 
call NUMBER (X,Y,0.15,WN(1)-032) 
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@rm 1 


381 


370 


700 
Ol 


702 
703 


704 
705 


706 
107 


e e s 


710 


if ((SYMBL .le. 13) .and. (SYMBL .ge. 0)) then 
call LINE (ATEMP,BTEMP,100,1,-1,SYMBL) 
end if 
continue 
INTRVL = INTRVL + 100 
continue 
call PLOT (0.0,0.0,999) 
FACT = 0.9 


Enter RANGE of ALPHA & BETA values to examine 


else if (EXPAND .eq. 1) then 
ea lily CER 
write(*,701) 
format (////,10x,’Enter Minimum ALPHA axis value ==> ‘”) 
call PSIT(6,48) 
read(*,*,err=700) EXPAMN 
write(*,703) 
format(/,10x,”’Enter Maximum ALPHA axis value ==> 7”) 
call PSIT(8,48) 
read(*,*,err=702) EXPAMX 
write(*,705) 
format(/,10x,’Enter Minimum BETA axis value ==> 7”) 
call PSIT(10,47) 
read(*,*,err=700) EXPBMN 
write(*,707) 
format(/,10x,”’Enter Maximum BETA axisS value ==> ‘”) 
Gall PSIT¢l2 .47) 
read(*,*,err=702) EXPBMX 


.. Calculate Constant ZETA Plot 


if (NUMZTA .ge. 1) then 
J = 0 
do 710 I = 1,STOP 
if((A(I) .ge. EXPAMN) .and. (A(I) .le. EXPAMX)) then 


if((B(I) .ge. EXPBMN) .and. (B(I) .le. EXPBMX)) then 


J=J+1 
EXPALP (J) 
EXPBET (J) 
end if 
end if 
continue 
FRSTZ = J +1 
DELTZ = J + 2 


A(T) 
B(I) 
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C 


call SCALE (EXPALP,7.0,J,1) 
eal] SCALE (EXPERI. 5.000. 


write 7 37) 
9falal format (///////,30x, “WINDOW EXPANSION’ ) 
write (*,712) EXPALP(FRSTZ) , EXPALP(DELTZ) 
712 format (/////,5x,°X Axis Min Value = “ ,#11.4,5x, X Axis DELTAg =e 
: oe) 
write (*,714) EXPBET(FRSTZ) , EXPBET(DELTZ) 
714 format (5x,°Y Axis Min Valuej= “ ,g11.3,5x,° Y Axis DELTA = {gine 
. (ie a 
write (*,716) 
716 format (77/7 ,25x, *** Calculatingsbata “** —=7o 


call PLOTS (0, IOPORT,MODEL) 

call FACTOR (FACT) 

call NEWPEN (1) 

call ASPECT (1.2) 

call STAXES: (G32. 20:35 70 see) 

call SYMBOL (STITLE, 6-0-2425. EL tLe, O.- NC) 

call AXIS (.8,.8,’Alpha’ ,-5,-7.0,0. ,EXPALP(FRSTZ) , EXPALP(DELTZ) ) 
call AXIS (.8,.8,’Beta’ ,4,-5.0,90. ,EXPBET(FRSTZ) , EXPBET(DELTZ) ) 


Draw dashed lines at axis increments and surrounding box 
call GRID 


INTSRT 1 
INTS EP 100 
do 7260 K = I NOMZIA 
J=0 
Jl = 0 
dov7g60 1 = INSERT. INioae 
if((A(1I) .ge. EXPAMN) .and. (A(I) .le. EXPAMX)) then 
if((B(I) .ge. EXPBMN) .and. (B(I) .le. EXPBMX)) then 
J=J+1 
ATEMP (J ) 
BTEMP (J) 
end if 
end if 
1616) continue 
ATEMP (J+1) 


A(I) 
B(I) 


EXPALP (FRSTZ) 

ATEMP(J+2) EXPALP (DELTZ) 

BTEMP(J+1) EXPBET (FRSTZ) 

BTEMP(J+2) = EXPBET(DELTZ) 

call CURVE (ATEMP,BTEMP,J,.1) 

if «((SYMBL .le. 13) .and. (SYMBL .ge. 0O)) then 
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call LINE (ATEMP,BTEMP,J,1,—-1,SYMBL) 
end if 
Jl = J-l 
call WHERE(X,Y,FACT) 
camel NUMBER @(X;Y,0.15,ZETA(K) ,O. v2) 
INTSRT = INTSRT + 100 
INTSTP = INTSTP + 100 
720 continue 
end if 


QO 


Constant Wn curves 


if (NUMWN .ge. 1) then 
if (NUMZTA .ge. 1) go to 819 
J =0 
do 810 I = 1,STOPO 
if ((AOMEGA(I) .ge. EXPAMN) .and. (AOMEGA(I) .le. EXPAMX) ) 
‘ then 
if ((BOMEGA(I) .ge. EXPBMN) .and. (BOMEGA(I) .le. EXPBMX)) 
* then 
J=J+1 
EXPALP (J) 
EXPBET (J) 
end if 
end if 
810 continue 
FRSTZ = J +1 
DELTZ = J + 2 


AOMEGA (T) 
BOMEGA (T) 


eal SCALE” (EXPAEP,7.0,J,1) 
call SCALE (EXPBET,5.0,J,1) 


write (*,811) 
811 format (///////,30x, “WINDOW EXPANSION’ ) 
write (*,812) EXPALP(FRSTZ) , EXPALP(DELTZ) 
bee tormae (/////,5x, X Axis®Min Value = “ ,g11.4,5x,”°X Axis DELTA = ” 
. Pai) ) 
write (*,814) EXPBET(FRSTZ) , EXPBET(DELTZ) 
814 format (5x,°Y AxiS Min Value = ” ,g11.3,5x,’Y Axis DELTA = ” ,gll1.4, 
ss eS 
write (*,816) 


oo, fOmmat (///,25x, *** Calculating Data ***’ ,//) 
c 


¢ 


call PLOTS (0, IOPORT ,MODEL) 
call FACTOR (0.9) 
call NEWPEN (1) 
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© 


819 


830 


820 


* 


call 
call 
eal 
call 
call 


Draw 
call 


INTSR 
INTST 
do 82 
J=0 
ili= 
do 

i 

£ 


ASPECT (132 

STAXIS 1, 1 otercO. 157/05 2) 

SYMBOL (STITLE,6.0, .25, TITLE,O.7, NC) 

AXIS (.8,.8,’Alpha’ ,-5,-7.0,0. ,EXPALP(FRSTZ) , EXPALP(DELTZ) ) 
AXIS (.8,.8,° Beta’ ,4,-5.0,90. , EXPBET(FRSTZ) , EXPBET(DELIZg® 


dashed lines at axiS increments and surrounding box 
GRID 


Hl 
— 


Ae 
P = 100 
O K = 1,NUMWN 


0 
830 I = INTSRT,INTSTP 
f ((AOMEGA(I) .ge. EXPAMN) .and. (AOMEGA(I) .le. EXPAMX)) 
hen 
if ((BOMEGA(I) .ge. EXPBMN) .and. (BOMEGA(I) .le. EXPBMxX) ) 
then 
J=J++1 
ATEMP (J) 
BTEMP (J) 
end. if 


AOMEGA (T) 
BOMEGA (T) 


end if 
continue 


ATEMP(J+1) 
ATEMP(J+2) 
BTEMP(J+1) 


EXPALP(FRSTZ) 
EXPALP (DELTZ) 
EXPBET (FRSTZ) 


BIEMP(J+2) = EXPBET (DELTZ) 

call CURVE (ATEMP,BTEMP,J,-.1) 

if ((SYMBL .le. 13) .and. (SYMBL .ge. 0O)) then 
call LINE (ATEMP,BTEMP,J,1,-1,SYMBL) 


end 
call 


tee 
WHERE (X,Y,FACT) 


call NUMBER (X,Y,0.15,WN(K),O. ,2) 


J1 


= J-1 


INTSRT = INTSRT + 100 


INTSTP 


INTSIE st aoe 


continue 
end if 


call 


else 


PLOT (0.0,0.0,999) 


go to 399 


end if 
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239 


399 


238 


—|———_ 
—— = 


call CLR 
write(*,239) 
format(1x,///////,9x,°Do you want to make another WINDOW selection 


~)° Weeox , Senter wyt) ope "n" ==> “) 


Cause PSIT (11,31) 

read(*,° (A)’ ,err=399) CHG 

if ((CHG .eq. ’y’) .or. (CHG .eq. “Y’)) then 
go to 30 

end if 

ai! CER 

write(*,238) 


format(1x,///////,9x,°Do you want to SAVE your constant CURVE sele 


Serrons? .//,0x., Enter "y" or "n" ==> 7 ) 


Galt PSIT(1i1,3l) 

read(*,° (A)* ,err=399) CHG 

if ((CHG .eq. ’n’) .or. (CHG .eq. “N’)) then 
NUMZTA = O 
NUMWN = O 
NUMZWN = 0 

end if 

return 

end 


Subroutine GRID -- draws dashed lines at axis increments 


Subroutine GRID 
integer EXPAND, SYMBL 
common/box/ EXPAND, SYMBL 


if (EXPAND .eq. 2) then 
Caml PLOT (0.8.,0.8;=3) 
else 
call PLOT (0.8,0.8,-3) 
endif 
call PLOT (0.0,5.0,3) 
@att PLOT (7.0,5.0;, 2) 
eail PLOT (7.0,0.0, 2) 


call STDASH (.01,.10) 

Cart PEOTD: (1.0,0.0,3) 
call PLOTD (1.0,5.0,2) 
Cail, FEOTID (2.0.5.0, 3.) 
call PLOTD (2.0,0.0,2) 
CalleerLOITD ($.0,0.0, 3) 
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call PLOIDV(3.-0.5.0.,c) 
eall PLOUTD (4.075.043) 
caleiry PLOTD (4.0,0.0,2) 
call PLOID (5.0,0-0,3) 
call PLOIFD (5..0,-57072) 
call PLOID (6.0,5.0,3) 
call PLOTD (620;0.0,2) 
call PEGID (7.0 ,1..0-3) 
call PLOITD (OFC. 1.0, 2) 
call PLOTD (0.0,2.0,3) 
call PLOTD (7.0,2.0,2) 
call PEOID (720,320, 3) 
call PEOTD: (070.320. 2) 
call PLOTD (0.0,4.0,3) 
call PLOTD (7.0,4.0,2) 

return 

end 

Subroutine PLTCRV -- calculates curves based on input arrays 


mm me mc cs em we se ws we ce we ww ms we ws we we es we ee we we ee ss ee me ee ee se ee ee es ee es ee aoe oe oe oe =. —— oe 
===. SS EE LEE Lee ee ee ee ee ee ee ee ee ee ee ae ee ee ee ee eae ae ae ae ae ae ae is i SS i SS i SS SS SS ee Se ee Se ee ere nel a ee 


Subroutine PLTCRV (STP,APTS,BPTS,NUMPTS,CRVTYP,CRVNUM, 


= TITLE, STITLE,NC,AFIRST,ADELTA, BFIRST, BDELMA® 
real ATEMP(102) ,BTEMP(102) ,APTS(1002) , BPTS(1002) ,CRVNUM(10) , 
ys STRITLE, CRVIYE JAP IRST PADELIAWVer IRST , BDELTA 


integer STP,FRST,DELT, INTRVL,NC,NUMPTS , EXPAND, SYMBL 
character*30 TITLE 

common/box/ EXPAND, SYMBL 

common/factr/ FACT 

common/symb/ SYMBL 


call FACTOR (FACT) 
call ASPECIMG. 2) 
call @STAKis= (213.20, .15 Oe 2a 
INTRVEE— © 
FRST = STP + 1 
DELT = STP + 2 
call SCALE (APTS,7.0,STP,1) 
ATEMP(101) = APTS(FRST) 
ATEMP(102) = APTS(DELT) 
AFIRST = APTS(FRST) 
ADELTA = APTS(DELT) 
write (*,120) ATEMP(101) ,ATEMP(102) 
120 format (///,5x,°X Axis Min Value = “ ,g11.4,5x,’X Axis DELTA = 7’, 


130 


: g11.4,//) 
call SCALE (BPTS,5.09STP;1) 
BTEMP(101) = BPTS(FRST) 
BTEMP (102) BPTS (DELT ) 
BFIRST BPTS(FRST) 
BDELTA = BPTS(DELT) 
write (*,130) BTEMP(101) , BTEMP(102) 


130 foemmat, (5x, Y AxiS Mim Value = ~ ,@11.3,5x,°Y Axis DELTA = ’, 
‘i oa / // 77) 
write (*,131) 
PS format (25x, *** Calculating Data ***° ,////) 


call AXIS (.8,.8,’Alpha’ ,-5,-7.0,0. ,APTS(FRST) ,APTS(DELT) ) 
coll AXIS (.8,.8, Beta ,4,-5.0,90.,BPIS(FRST) ,BPTS(DELT) ) 
Gabtl SYMBOL (STITLE,6.0,.25,TITLE,O. ,NC) 


O 

C .. . Draw dashed lines at axis increments and surrounding box 
call GRID 

c 


do 140 I = 1,NUMPTS 
Gowr1 5060" = 100 
ATEMP(J) = APTS(J + INTRVL) 
BTEMP(J) = BPTS(J + INTRVL) 
150 continue 
call CURVE (ATEMP,BTEMP,100,CRVTYP) 
call WHERE (X,Y,FACT) 
call NUMBER (X,Y,0.18,CRVNUM(T) ,0O. ,2) 
if ((SYMBL .le. 13) .and. (SYMBL .ge. 0O)) then 
call LINE (ATEMP,BTEMP,100,1,-1,SYMBL) 
end if 
do 161 K = 1,100,99 
ATEMP(K) = ATEMP(K) / APTS(DELT) 
BTEMP(K) = BTEMP(K) / BPTS(DELT) 
call NUMBER (ATEMP(K) , BTEMP(K) ,0.18,CRVNUM(I) ,O. ,2) 
cl6l continue 
INTRVL = INTRVL + 100 
140 continue 
C 


QO GOOG 


return 
end 
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KK KKK KOK OK KOK OK KKK KOK KOK KKK OK OK KOK KOK KOK KOK KKK KOK KOK KKK KKK KOK KOK KOK KKK KOK KOK KOK KK KK KK KKK K 


ASCII MODULE 


MK KK kK KK OK OK KKK OK KK OK ok OK OK OK OK OK OK OK OK OK KK OK OK KOK OK OK OK OK OK OK OK KOK OK OK OK KOK KOK OK KK OK OK OK OK OK OK KOK OK KOK KOK KK KK 


OLS SS a @, 


C ssssssssssssssseeee2222e2es222ee=2222252525252255SSSs22e22SSeee=s=555S5 
C Subroutine CLR -- clears the screen 
C sensansn2enseeeee eee sees ee esses esse 282SSSe58882 855888252 2555222252== 
C 

Subroutine CLR 

character*i Ci. C2.) C3, C4 

integer IC(4) 

equivalence (C1,IC(1)),(C2,1C(2)),(C3,1C(3)), (C4 , 2ea 458 

data IC/16#41B,1645B,16432,16#4A/ 

write(*,1) C1,C2,C3,C4 
al format (1xX,4Al1) 

return 

end 
C 
SS SS SS SS SS SS SS SS SS SSS SS SS 
C Subroutine PSIT -- positions cursor by row and column 
C s===$==$=$=$-22=2s2=<222=s2222=2222Ss22S2Se22eee5S=e82228222252222222===== 
C 

Subroutine PSIT( ROW, COLUMN) 

integer IC(4),ROW,COLUMN,L 

eharacter*t| Cl. C23@57Cs-Le ro 

cnaracter*s CBUFF 

equivalence (C1,1C(1)),(C2,1C(2)),(C5,1C(3) )\9ee@er1C(4)), 

~ (CBUFF ,LC(1) ) 

data IC/16#1B,16#5B,16#43B,16#66/ 
C 

L=10000+100* ROW+COLUMN 
C 


C +++ Write EScape Codes to a Character Buffer +++ 
write(CBUFF,2) L 
2 format(I15) 


C +++ Write Escape Codes to Display +++ 
Write(*.5)) Cl,62,LC(2).L6 (see, be(4) Ree Gmc 
3 format(1X, 8Al1,\) 
return 
end 
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